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The  final  research  report  containes  studies  of  rare  earth  (RE)  centers  in  semiconductors, 
namely  Yb,  Nd,  and  Tm  in  GaAs,  InP,  ZnS  and  CdS.  Special  emphasis  has  been  put  on  the 
optical  properties  of  the  RE  centers  in  semiconductors  systems  by  studing  the  luminescence  rise 
time  of  the  rare  earth  doped  semiconductors,  excited  indirectly  above  the  band-gap,  contain 
information  about  the  energy  transfer  processes  from  the  host  to  the  4P  electron  system. 

The  results  are  collected  in  the  following  paragraphs  and  papers: 

I  Introduction 

II  Experimental  Techniques 

A.  Sample  Preparstion  Ion  Implantation 

B.  Annealing 

C.  Spectropolarimeter 

III  Photoluminescence  Study  of  GaAs:Er  diffused  with  Li 

A.  Introduction 

B.  Sample  Preparation 

C.  Experimental  Techniques  and  Results. 

D.  Summary  and  Conclusion 

IV  Electroluminescence  of  Tm  Doped  ZnS 

V  Kinetics  of  luminescence  of  isoelectronic  rare  earth  ions  in  III-V  semiconductors 

VI  Kinetics  of  the  Luminescence  of  Isoelectronic  Rare-Earth  Ions  in  III-V  Semiconductors 

VII  Luminescence  Properties  of  Yb  Doped  InP 

VIII  Luminescence  Properties  of  Nd  and  Yb  doped  CdS 

IX  Digital  spectropolarimeter  for  the  measurement  of  optical  polarization 

X  Optical  properties  of  CdTe/Cd^i^Te  strained-layer  single  quantum  wells 
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I  Introduction 

This  final  research  report  containes  studies  of  a  rare  earth  (RE)  centers  in  a  semiconductors, 
namely  Yb,  Nd,  and  Tm  in  GaAs,  InP,  ZnS  and  CdS.  Special  emphasis  has  been  put  on  the 
optical  properties  of  the  RE  centers  in  semiconductors  systems  by  studing  the  luminescence  rise 
time  of  the  rare  earth  doped  semiconductors,  excited  indirectly  above  the  band-gap,  contain 
information  about  the  energy  transfer  processes  from  the  host  to  the  4f”  electron  system.  It  is 
shown  that  the  study  of  the  rise  time  at  different  temperatures,  excitation  intensities,  and  excitation 
pulse  durations  can  provide  important  information  about  the  energy  transfer,  the  radiative,  and 
non-radiative  processes  respectively.  We  discuss  only  the  structured  isoelectronic  traps  in  III-V 
semiconductors  introduced  by  triply  charged  rare  earth  ions  replacing  the  element  from  column 
in  (or  athers  more  complex  RE  isoelectronic  traps).  Furthermore,  we  develop  the  luminescence 
kinetic  models  (for  InP:  Yb  which  can  be  applied  toather  semiconductors)  that  describe  the  energy 
transfer  and  recombination  processes.  The  numerically  simulated  luminescence  rise  and  decay 
measurements  for  InP:  Yb  show  a  good  quantitative  agreement  with  experiment  over  a  wide  range 
of  generation  rate.  The  presence  of  low  lying  empty  core  orbitals  in  rare  earth  impurities 
introduces  new  excitation  and  recombination  phenomena,  which  we  discussed  in  detail.  Finally 
the  possible  quenching  mechanisms  and  the  temperature  dependence  of  rare  earth  luminescence 
are  discussed. 

The  study  of  photoluminescence  (and  electroluminescence)  of  some  II-VI  compounds  include: 
ZnS:  Tm3+,  CdS:  Yb3+,  CdS:  Nd3+.  The  deep  understanding  of  excitation  mechanism  in  II-VI 
compounds  system  can  lead  to  better  understanding  of  excitation  mechanism  in  other  more 
complicated  RE3+  systems  in  semiconductors.  Figure  la  shows  the  energy  levels  scheme  of  the 
Nd3+  ion  in  cubic  (Td)  and  noncubic  symmetry  (Cv3,  Cv2,  CJ  environment.  In  cubic  symmetry 
4p3/2  level  is  not  split  (rg  fourfold  degenerated  crystal  field  state).  Therefore,  in  the  cubic  hosts, 
any  splitting  of  the  ^3/2  level  is  indicative  for  association  of  the  Nd3+  ion  with  other  defects.  The 
same  argument  holds  for  the  multiplicities  of  the  ground  state.  The  ground  state  4  l9n  in  cubic 
symmetry  is  split  into  three  states  (Tg,  rg,  r7)  and  the  first  excited  state  \m  into  four  crystal  f-leld 
states  as  in  the  left  part  of  Fig.  la.  In  noncubic  symmetry,  the  fourfold  degenerate  rg  states  of 
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Nd3+  ion  are  split  into  two  doublet  each  and  the  energy  level  structure  is  shown  in  the  right  part 
of  Fig.  la.  The  emission  spectra  for  the 4  F3/2  - 4 transitions  of  Nd3+  in  noncubic  symmetry 
should  exhibit  all  ten  transitions.  At  low  temperature,  when  only  the  lower  level  of  4  P3/2 
occupied,  only  5  lines  should  be  observed  in  luminescence  emission.  In  absorption  at  low 
temperature  such  that  only  the  lowest  component  of 4  4/2  state  is  occupied,  only  two  lines  should 
be  observed.  In  Fig.  lb  we  show  the  possible  absorption  and  emission  transition  with  wavelength 
(in  urn)  for  Nd 3+  in  cubic  symmetry  in  semiconductors  with  band  gap  E  g  ^  2.14  eV.  Figure  1 
c  shows  energy  level  schemes  for  Yb3+  ion  in  cubic  (Td)  and  noncubic  symmetry,  and  the 
luminescence  transition  for  Tm3+  ion  in  ZnS  is  shown  in  Fig. I  d. 

During  the  period  of  Fall  1991  -  Fall  1992,  I  was  on  my  Faculty  Fellowship  Leave 
sponsored  partially  by  Ohio  University  Faculty  Fellowship,  the  Air  Force  Office  of  Scientific 
Research,  and  Nippon  Telegraph  and  Telephone  Basic  Research  Laboratory  in  Tokyo  Japan. 

The  first  part  (4  months)  of  my  sabbatical  leave  I  spent  at  Wright  Laboratory  Solid  State 
Electronics  Directorate  (WL  /  ELRA),  Wright-Patterson  AFB  in  Dayton  Ohio.  I  collaborated  with 
Dr.  K.  Evans  (and  his  co-workers),  and  Dr.  Reynolds.  Part  of  published  paper  (4,  see  Ch.VHI) 
namely  the  streak  camera  measurement  was  done  in  Dr.  Reynolds  laboratory  at  AFB  in  Dayton. 
During  this  period  I  also  worked  on  the  theory  of  indirect  photoluminescence  excitation  of  rare 
earth  ions  in  ni-V  semiconductors.  I  spent  the  second  part  of  my  Faculty  Fellowship  Leave  (four 
and  1/2  months)  at  the  Nippon  Telegraph  and  Telephone  Basic  Research  Laboratory  in  Tokyo, 
Japan.  I  collaborated  with  Dr.  K.  Takahei,  M.  Taniguchi,  and  A.  Taguchi,  continuing  the 
research  effort  initiated  at  Dayton,  on  EL  structures.  The  next  two  projects  are:  1)  Study  the 
energy  transfer  mechanism  from  the  host  to  rare  earth  isoelectronic  traps  in  MOCVD  grown 
GaAs:  Nd,  and  InP:  Yb  by  investigating  the  PL  rise  kinetics  as  the  function  of  excitation  intensity 
(in  preparation).  2)  Study  of  the  lithium  influences  on  the  photoluminescence  spectra  of  MOCVD 
grown  GaAs:  Er  sample  (see  Chapter  III). 

During  my  sabbatical  leave  I  presented  several  seminars  on  photoluminescence  and  its 
kinetics  of  rare  earth  doped  semiconductors,  at  Wrigth  Laboratory  and  at  NTT  Tokyo. 
Currently,  Dr.  T.  Li  and  I  collaborate  with  Drs.  D.  Reynolds,  and  K.  Evans  on  polarization 
spectroscopy  of  coupled  double  quantum  wells  AlGaAs  /  GaAs  (research  in  progress). 
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To  study  the  electroluminescence,  and  especially  the  direct  impact  excitation  mechanism 
of  rare  earth  in  III-V  semiconductors  I  designed  the  four  EL  structures  and  Dr.  K.  Evans  with 
coworkers  grew  them  using  MBE  system.  (The  grown  structures  unfortunately  were  not  doped 
with  Er  caused  by  depleted  Er  source  in  MBE  machine).  The  reason  why  we  are  trying  to  excite 
the  RE  ions  by  direct  impact  is  that,  the  RE3+  ions  cannot  be  effectively  excited  by  energy 
transfer  from  host  to  RE  core  states.  We  believe  that  this  is  a  reason  the  efficiencies  of  the  III-V: 
RE3+  crystals  are  low.  Some  materials  such  as  GaAs:  Yb  and  InP:  Nd  does  not  show 
luminescence  under  band  to  band  host  excitation.  The  direct  excitation,  either  optically  or 
electrically  is  me  only  way  to  excite  these  ions.  We  would  also  like  to  prove  it  experimentally. 
Some  of  the  hetrostructures  proposed  below  are  the  devices  where  the  carriers  can  gain  energy 
from  the  bandgap  discontinuity  (as  was  proposed  in  my  publications  [1,2]).  Also,  by  using 
superlattice  structure  AlGaAs  /  GaAs:  Er  /  AlGaAs,  we  could  manipulate  the  band  gap  energy. 
In  each  case,  more  upper  states  of  rare  earth  RE3+  (Er)  could  be  used  as  the  pumping  levels  that 
might  result  in  a  higher  luminescence  (EL,  PL)  efficiency.  It  will  be  important  to  compare  this 
direct  impact  excited  devices  with  the  conventional  LEDs  to  see  the  difference  in  excitation 
mechanisms  and  to  see  which  of  them  are  more  efficient.  It  is  interesting  to  understand  what 
happens  physically  when  the  rare  earth  ions  collide  with  hot  carriers.  The  rare  earth  ions  could 
be  excited  by  the  following  mechanisms:  (a)  A  hot  electron,  with  energy  equal  to  the  RE3+ 
4f-shell  transition  energy,  collides  with  the  RE3+  ions  and  directly  transfers  its  energy  to  the  RE 
4f-shell  electron,  (b)  If  an  electron  gains  enough  energy,  it  might  ionize  a  rare  earth  ion.  The 
ionized  RE3+  ions  can  then  recapture  another  electron  with  RE  4f-shell  in  the  excited  state,  and 
(c)  The  high  electric  field  creates  hot  electrons  which  by  collisions  with  the  lattice  can  generate 
electron  hole  pairs.  In  such  a  case,  the  RE3+  ions  could  be  excited  through  host  energy  transfer 
mechanism,  similar  to  above  bandgap  photoluminescence  (PL)  excitation. 

With  a  careful  study,  we  can  determine  which  excitation  mechanism  is  dominant  one,  and 
can  probably  obtain  information  about  the  hot  electron  distribution.  The  most  likely  excitation 
mechanism  is  the  (a)  mechanism  that  uses  the  hot  carriers  (electrons)  with  energy  equals  or  higher 
than  the  RE  4f-shell  internal  transition  energy.  The  direct  inelastic  collision  of  such  hot  electrons 
with  the  RE3+  ions  transfers  its  energy  to  the  4f  electron,  and  creates  the  excited  RE3+  ions.  The 
designed  electroluminescence  structures  are: 
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1.  Figure  2a  shows  the  metal  -  semiconductor  (doped  with  RE)  -  insulator  -  semiconductor  (doped 
with  RE)  -  metal  (MSISM)  structure.  The  device  operates  based  on  tunneling  of  electrons  through 
a  thin  (  ~150A-200A )  dielectric  layer.  In  the  designed  device  based  on  GaAs  instead  of  insulators 
with  very  wide  bandgap,  we  used  undoped  Al^Gao  SJAs  -  200A  thick  layer.  On  both  sides  of  the 
AlGaAs  layer,  we  have  50A  GaAs  layers  and  0.4  ftm  thick  GaAs  layer  doped  with  Er  and  Si 
(lightly  n-type)  covered  by  heavily  doped  (n+)  GaAs:  Si  layer.  The  ohmic  contacts  were  applied 
to  both  n+  GaAs  sides  completed  the  structure.  The  MSISM  structure  is  symmetric  and  can 
operate  by  applying  AC  voltage  as  shown  in  Figure  2b.  The  hot  electrons  tunneled  through  the 
dielectric  layer,  (where  most  of  the  applied  voltage  dropped)  should  have  sufficient  energy  to 
excite  by  impact  the  RE  centers.  The  structure,  should  operate  at  low  voltage. 

Figure  3  shows  the  asymmetric  structure  in  which  the  left  side  is  the  same  as  in  Figure  2  while 
the  right  side  contains  the  GaAs  p-n  junction  with  very  thin  (  *  150A)  p  region.  When  the  p-n 
junction  is  forward  bias,  an  electron  will  be  injected  in  the  conduction  band  of  the  p-layer.  The 
hot  electrons  will  be  tunneled  across  the  AlGaAs  layer  into  GaAs:  Er,  Si  layer  ,  and  by  direct 
impact  will  excite  the  RE3+  ions. 

2.  The  p+-i-n+  structure  grown  by  MBE  on  an  n+  GaAs  substrate  is  shown  in  Figure  4.  The 
p+-i-n+  diode  consists  of  three  semiconductor  regions:  the  2500  A  thick  p+  Al0  45Gao55As:  Be  (* 
5xl019  cm'3),  the  n+  GaAs:  Si  region  and  the  semi-insulating  linearly  graded  GaAs- AlGaAs  *i’ 
region  of  width  d  doped  with  Er  (  *  6xl018  cm'3).  The  three  regions  are  capped  with  a  50A  thick 
p+  GaAs:  Be  layer.  By  applying  the  reverse  bias,  the  depletion  region  and  internal  field  is 
forming  in  the  graded  'i'  region.  With  a  sufficiently  large  voltage,  the  depletion  region  can  punch 
through  from  p+  to  n+  layer  as  shown  in  Figure  4.  The  injected  electrons  (by  light)  through  the 
thin  p+  region,  are  accelerated  in  a  strong  and  well-defined  electric  field,  and  can  gain  energy 
equal  to  or  larger  than  the  RE  4f-shell  internal  transition  energy.  If  such  a  hot  electron  collides 
with  the  RE3+  ion,  it  can  directly  transfer  its  energy  to  the  RE  4f-shell  electrons  and  therefore 
creates  an  excited  RE3+  centers. 


n*-4oaA 

G*As 


n+  -400A  A1(.45)G*  (.55)  As 

GaAs  or 

Insulator 


Figure  3 


Figure  4  The  p-i-n  impact  light  emitting  device 
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3.  The  AlGaAs/GaAs:  Er/AlGaAs  superlattice  structure.  For  the  first  time  we  designed  the 
unipolar  n+  -  superlattice  -  n+  structure  as  shown  in  Figure  5.  The  GaAs:  Er/Al^jGaosjAs 
superlattice  was  grown  by  MBE  on  an  n+  GaAs:  Si  substrate.  It  consisted  of  60  periods  of 
nominally  350A  wide  GaAs:  Er  (  "  2xl0‘®cm'3 )  quantum  wells  and  100  A  thick  Al^Ga^As 
barriers,  sandwiched  between  450  A  thick  GaAs:  Er,  Si  and  0.5  ptm  thick  n+  GaAs:  Si  cap  layers. 
The  calculation  made  for  the  above  super  lattice  shows  that  because  of  thick  barriers  (100  A  ),  the 
confined  levels  are  localized  states  rather  than  extended  subbands.  In  Figure  6,  we  show  the 
conduction  and  the  valence  band  edge  of  the  superlattice  doped  with  Er3+  (  2xl018  cm'3 ),  at  a  bias 
voltage  above  the  impact  excitation  threshold  of  4f  electrons.  After  further  increase  of  the  external 
bias,  some  electrons  are  accelerated  until  they  will  reach  the  threshold  energy  for  impact 
generation  of  electrons-hole  pairs,  and  electroluminescence  of  the  superlattice  should  be  observed. 

The  light  emitting  structure  described  in  Figure  4  required  for  operation  injection  of 
electrons  to  region  "i."  It  can  be  done  by  light  generation  of  electron  hole  pairs  in  p+  region.  To 
overcome  this  drawback,  I  designed  during  my  stay  at  NTT  Tokyo  new  structure,  the  light 
emitting  transistor  (LET).  The  new  LET  structure  is  shown  in  Figure  7  (with  and  without  bias). 
If  the  emitter  base  heterojunction  is  forward  biased,  the  electrons  are  injecting  through  the  base 
region  to  "i"  region  of  the  reverse  bias  p-i-n  structure.  The  electrons  are  accelerated  in  a  strong 
and  well-controlled  electric  field  existing  in  Mi"  region.  The  hot  electron  by  direct  impact  with 
RE3+  ions,  excite  the  core  states  of  rare  earth  ions. 

1.  H.  J.  Lozykowski,  Solid  State  Comm.  66,  755  (1988) 

2..  H.  J.  Lozykowski,  in  "Electroluminescence"  Edited  by  S.  Shionoya  and  H.  Kobayashi, 

Springer  Proceedings  in  Physics,  38,  60  (1989) 

II  Experimental  Techniques 

A.  Sample  Preparation  Ion  Implantation 

The  ion  implantation  of  rare  earth  was  done  by  Dr.  I.  G.  Brown  from  Lawrence  Berkeley 
Laboratory  at  the  University  of  California,  Berkeley.  Implantation  of  Nd  '  Yb,  and  Tm  was 
performed  with  an  unconventional  implantation  technique  using  a  new  kind  of  high  current  metal 
ion  source  developed  at  Berkeley  by  Dr.  Brown  et  al.  This  Metal  Vapor  Vacuum  Arc  (MEW A) 


Supcriattice  (60  periods) 


Figure  5  AIGaAs/GaAs:Er/AIGaAs,  60  periods  superlattice  doped  with  Er 


Figure  6  Conduction-  and  valence-band  edge  of  the  superlattice  doped  with  rare-earth  at  a  bias 
voltage  above  the  impact  excitation  threshold  of  rare  earth  ions.  The  arrows  indicate  electron 
transport  the  impact  excitation  and  recombination  process. 


Figure  7  Light  Emitting  Transistor 
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ion  source  makes  use  of  a  dense  metal  plasma  generated  directly  from  the  metal  cathode  material 
as  the  medium  from  which  ions  are  extracted.  The  source  is  repetitively  pulsed  and  the 
implantation  is  done  on  a  broad  beam  mode,  with  a  direct  line  -of-sight  from  the  ion  source  to  the 
target.  High  current  beams  of  a  wide  range  of  multiply-charged  metal  ions  can  be  produced. 
Beam  extraction  voltage  is  up  to  100  kV,  for  a  mean  ion  energy  of  up  to  several  hundreds  of  key 
by  virtue  of  a  mean  ion  charge  state  that  is  typically  2  to  3.  The  beam  current  can  be  up  to  several 
amperes  peak  and  around  10  mA  time  averaged  delivered  onto  target.  The  source  has  been 
operated  with  virtually  all  of  the  solid  metals  of  the  periodic  table.  The  high  energy  implants  of 
1  MeV  were  implanted  at  Universal  Energy  Systems  (LYES),  Dayton,  OH. 

B.  Annealing 

The  samples  were  cut  to  small  pieces  (4  mm  x  4  mm  squares)  and  then  cleaned  by  standard 
methods.  Thermal  heat-pulse  using  a  strip  heater  was  used  to  anneal  the  samples  up  to  900°  C  for 
10- 15  sec.  The  annealing  process  was  performed  using  the  apparatus  shown  in  Fig.  8  which  we 
designed  and  built  in  our  laboratory.  A  flat  graphite  plate  of  3  mm  thick  is  placed  on  the 
tantalum  strip  heater.  The  samples  are  placed  with  the  implanted  face  down  on  the  graphite  plate 
and  covered  with  a  graphite  lid  of  1  mm  thickness.  The  temperature  is  measured  by 
microprocessor  controlled  temperature  indicator  OMEGA  model  650  through  a  thermocouple 
placed  inside  the  bottom  graphite  plate.  The  variation  of  annealing  temperature  with  time  is 
recorded  by  an  y-x  recorder  as  shown  in  the  insert  of  Fig.  1.  The  rise  time  of  the  temperature  was 
controlled  by  the  AC  voltage  applied  to  the  strip  heater.  For  temperatures  up  to  900  °C,  our 
system  maintains  a  rise  time  of  6  sec  and  a  decay  time  from  900  °C  to  500  °C  of  22  sec. 

C.  Spectropolarimeter 

For  photoluminescence  (electroluminescence)  polarization  measurements  we  modified  our 
digital  spectropolarimeter  to  operate  at  the  spectral  range  from  450  nm  to  2000  nm  (in  digital 
mode  the  range  from  450  nm  to  1200  nm  is  limited  by  the  photomultiplier  spectral  response.)  The 
modified  optical  range  was  accomplished  utilizing  Oriel  Fresnel  rhomb  achromatic  retarders.  The 
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Figure  1  Thermal-Heat-Pulse  Annealing  System 


quarter  wave  or  the  half  wave  achromatic  retarders  are  in  series  with  a  Glan-Thompson  crystal 
polarizer  in  the  excitation  channel  to  provide  90°  ±  2°  (o+  or  o')  or  180°  ±  2°  (linear 
polarization)  phase  retardation.  The  heart  of  the  setup  consists  of  a  photo-elastic  modulator 
(PEM-80)  and  a  dual  channel  gated  photon  counter  which  has  the  capability  to  compensate  the 
variation  of  excitation  intensity  with  time  and  wavelength.The  second  important  modification  is 
automatic  wavelength  tracking  of  the  modulator  to  insure  the  constant  retardation  phase  (it/2  or 
x/4)  when  the  wavelength  of  the  analyzed  light  is  scanned  in  wide  range.  The  angular  aperture 
of  PEM-80  modulator  comprises  a  50  degree  full-cone  angle.  Unlike  the  conventional  lock-in  or 
other  analog  systems,  where  the  difference  of  two  oppositely  polarized  light  components  (AC)  and 
the  average  of  these  components  (DC)  are  extracted  electronically,  and  the  ratio  AC/DC  is 
recorded  as  degree  of  polarization,  our  setup  directly  records  theo+  and  o'  components  (or,  in  the 
case  of  linear  polarization,  I„  and  Ll  components)  of  the  light  signal  in  the  same  run  (of 
wavelength,  temperature,  etc.)  A  detailed  description  of  the  spectropolarimeter  system  is  in  paper 
Chapter  VIII. 
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III  Photoluminescence  Study  of  GaAs:Er  diffused  with  Li 
A.  Introduction 

II-VI  semiconductors  containing  rare  earth  activators  have  been  of  interest  for  many  years. 
Such  materials  have  a  discrete  ion  as  the  principle  source  of  luminescence  emission,  the  4f-4f 
transitions  yielding  narrow  emission  band  or  lines.  For  all  earliest  investigations  of  rare  earth  in 
II-VI  compounds  see  review[l].  The  4f-4f  line  emissions  from  many  rare  earths  co-doped  with 
Cu  and  Ag  are  accompanied  by  characteristic  broad  band  emission  of  the  host  material.  In  1965, 
Yocom  and  Larach  [2]  showed  that  the  effect  of  charge  compensation  particularly  by  alkali  metals 
reduces  the  band  emission  of  ZnS  and  ZnCdS  doped  with  rare  earth  and  enhances  the  rare  earth 
emission  at  room  temperature.  Shrader  [3],  extending  the  work  of  Yocom  and  Larach,  reported 
on  the  details  of  charge-compensation  effect  in  ZnS  doped  with  dysprosium. 

The  proposed  model  [1-3]  of  rare  earth  center  in  II-VI  compounds  assumed  that  the  trivalent 
rare  earth  was  in  a  substitute  lattice  site  for  Zn  with  a  near-neighbor  vacancy  being  present  for 
charge  compensation  purposes.  They  assumed  that  such  a  rare  earth-vacancy  system  resulted  in 
broad  band  emission,  and  that  filling  the  vacancy  by  simplest  type  of  cationic  charge  compensating 
species  chosen  from  the  group  IS  ions  (the  alkali  metal  ions)  in  particular  lithium,  would  result 
in  strong  narrow  line  emission  normally  expected  from  4f-4f  transition.  The  experimental  results 
indeed  have  shown  that  the  addition  of  alkali  metals  ions  to  ZnS  doped  with  rare  earths  decreases 
or  eliminates  the  band  emission  of  the  host  crystals  while  selectively  intensifying  the  sharp 
emission  line  of  4f-4f  transitions. 

The  fluoride  ion  also  reduced  the  broad  band  emission  and  enhanced  strongly  the  line  emission 
of  the  rare  earth.  In  this  case  (fluoride),  it  is  assumed  that  interstitial  fluorides  compensate  the  rare 
earth  eliminating  the  need  for  a  vacancy.  Despite  the  model  for  compensation  mechanism,  the 
several  crystallographically  nonequivalent  positions  exist  for  all  the  possible  compensators  to  be 
found.  The  luminescence  emissions  from  such  nonequivalent  sites  cause  the  line  broadening  and 
the  appearing  of  extra  emission  lines.  The  rare  earth  alkali  metal  ions'  complex  centers  are 
responsible  for  the  drastic  intensification  of  certain  4f-4f  transitions,  as  opposed  to  the  non-alkali 
containing  materials. 

Studies  of  Li  in  GaAs  were  quite  extensive  in  the  early  seventies,  and  they  have  provided 
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information  about  solubility,  diffusion  data  and  local  mode  energies  from  IR  spectroscopy  [4]. 
Recently  photoluminescence  studies  of  Cu^-Lij  complexes  and  double  ionize  Lio,  acceptor  [5], 
and  influence  of  Li  incorporation  on  luminescence  of  n-type  Si,  Se  or  Sn  doped,  and  p-type  Zn 
doped,  or  SI  GaAs  ware  reported  [6a, b].  The  interactions  of  Li  and  Cu  in  GaP  create  the  neutral 
complexes  of  isoelectronic  character  [7,8].  Both  these  group-I  species  are  fast  interstitial  diffusers, 
and  in  complexes  both  can  occupy  Ga  sites.  At  least  five  different  Cu-Li  complex's  ware 
observed  in  GaP  after  a  subsequent  Li  diffusion  into  Cu  doped  material.  These  centers  are 
characterized  by  their  different  bond  exciton  luminescence  spectra,  all  of  them  being  neutral 
associates  of  isoelectronic  molecular  type.  In  addition  Li  form  neutral  linear  Lij-Li^-Op  <  1 1 1  >  - 
oriented  defect  in  GaP,  which  has  been  studied  via  its  bound  exciton  spectra  [9]. 

In  view  of  the  above  strong  enhancement  role  of  Li  on  luminescence  properties  of  rare  earth 
in  II-VI  compounds,  and  involve  in  creation  "isoelectronic"  defects,  we  focus  our  attention  on  the 
behavior  of  Li  in  GaAs  doped  with  Er.  To  our  knowledge,  no  study  on  influence  of  Li  on  Er 
doped  GaAs  has  been  reported  in  the  literature.  Below,  we  present  for  the  first  time 
photoluminescence  spectra  for  GaAs:  Er,  Li.  In  section  B  we  describe  the  Li  doping  procedure 
used  in  this  work  emphasizing  the  different  types  of  MOCVD  samples  used  to  avoid  systematic 
errors.  We  also  briefly  describe  the  experimental  techniques  employed  (section  C).  Section  D 
contains  the  experimental  results  and  discussion. 

B.  Sample  Preparation 

The  GaAs  :  Er  sample  used  in  this  study  was  grown  by  MOCVD  at  two  different  temperatures 
on  a  non  doped  GaAs  substrate  (orientation  (100)).  The  epitaxial  layers  of  about  1.3  pm  and  1.1 
pm  with  Er  concentration  1.5  x  1018  /cm3  were  grown  at  temperature  500°C  and  600°C 
respectively.  Lithium  diffusions  were  done  into  GaAs:  Er,  and  into  n-type  and  p-type  MOCVD 
epitaxial  layer  at  different  temperatures.  Table  I,  II  and  in  summarize  the  properties  of  the 
starting  materials  and  Li  diffusion  conditions.  The  Li  (purity  99.9  %)  diffusion  was  done  in  sealed 
quartz  ampules,  sample  was  in  small  crucible,  the  Li  metal  foil  and  As  (purity  99.9999  %)  were 
in  additional  quartz  liner  with  one  end  closed  in  order  to  protect  the  ampule  from  chemical 
reaction  with  Li.  The  ampules  and  Li  foil  were  cleaned  in  toluene  (ultrasonic  base),  alcohol  and 
deionized  water  and  evacuated  to  4  x  1C7  mbar.  The  As  metal  was  added  for  protection  of  GaAs 
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from  changing  stoichiometry.  The  samples  were  rapidly  cooled  to  room  temperature  after  heat 
treatment  by  immersing  the  quartz  ampule  into  water.  The  reference  samples  were  given  identical 
heat  treatment  as  the  Li-diffused  one  but  without  the  Li  source  and  only  with  As. 


Table  I 


Sample 

[substrate 

GaAs  SI 

(100)] 

Dopant:  Er 

concentration 

Epitaxial 

layer 

thickness 

Growth  temperature 

V  /  III  ratio  "30 

P52 

(1, 2,3,4) 

1.5  1018/cm3 

1.3  pm 

500  °C 

P4X 

(1,2, 3,4) 

1.5  018/  cm3 

1.1  pm 

600  °C 

Sample  P52  &  P4X 


6 

5 

1 

2 

3 

4 

Flow  direction  =» 


thicker 


thinner 


Table  n 


Sample 
GaAs  (SI) 
(100) 


P52(l) 


P52(2)a 


P52(2)b 


P52(3) 


P52(4) 


P52(5) 


P52(6) 


Dopant :  Er 
concentration 


1.5  0"/  cm3 


1.5  0“  /  cm3 


1.5  018  /  cm3 


1.5  018  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


1.5  0"  /  cm3 


Epitaxial 

layer 

thickness 


1.3  pm 


1.3  pm 


1.3  pm 


1.3  pm 


1.3  pm 


1.3  pm 


1.3  pm 


Diffusion  conditions, 
(all  sample  cooled  in 
water) 


Li + As, 570  °C,  6h 


Li+As,500°C,  6h 


Li+As,500  °C,  6h  + 
LH-As,802°C,lh 


Li+As,646  °C,  6h  + 
Ii+As,802  °C,  lh 


As,  802  °C,  lh 


P4X(1) 


P4X(2) 


P4X(3) 


P4X(4) 


P4X(5) 


P4X(6) 


1.1 

pm 

1.1 

pm 

1.1 

pm 

1.1 

pm 

1.1 

pm 

1.1 

pm 

Li+As,802  °C,  6h 


Li + As, 500  °C,  6h 


Li+As,646  °C,  6h 


As,  802  °C,  lh 


Table  III 


Sample 

GaAs  (100) 

Dopant  and 

concentration 

Epitaxial 

layer 

thickness 

Diffusion  conditions, 

(all  sample  cooled  in 

water) 

P5E 

(sub.:Si) 

Se  7.7  1  018/cm3 

1.2  pm 

Li+As,  800  °C,  lh 

P5F  (sub. 

undoped) 

Se  7.7  1  0“/cm3 

1.2  pm 

Li+As,  802  °C,  80 

min. 

P5X  (sub. 

undoped) 

Zn  1.8  1019/cm3 

2.3  pm 

Li+As,  802  °C,  lh 

P5Y  (sub. 

undoped) 

Zn  1.4  10,7/cm3 

2.3  pm 

Li+As,  802  °C,  lh 

3.5  pm 

Li,  500  °C,  7h 

3.5  pm 

Li,  570  °C,  6h 

3.5  pm 

Li,  409  °C,  72h 

3.5  pm 

As,  570  °C,  6h  (ref.) 

3.5  pm 

Li+As,  800  °C,  80 

min. 

P2A(1) 

(sub. 

undoped) 


P2A(2) 


P2A(3) 


P2A(4) 


P2A(5) 


n  =  1.5  0,5/cm3 
p=4.1103[cm2/Vs] 
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C.  Experimental  Techniques  and  Results. 

The  photoluminescence  measurements  were  done  at  2  K  and  at  77  K.  For  excitation,  the 
632.8  nm  line  of  He-Ne  laser  (lOmW)  and  the  5 14.5  nm  line  of  Ar+  laser  (20  mW)  were  used. 
The  photoluminescence  signal  was  analyzed  with  a  1.25-m  monochromator  and  was  detected 
with  a  Ge  photodetector  cooled  to  77  K.  Figure  1  shows  an  example  of  SIMS  bulk  analysis  of 
two  epilayers  grown  by  MOCVD  and  diffused  with  Li.  The  concentrations  of  Li  are  3  -  4  x 
101®  cm'3  to  4  x  10'®  cm’3  dependent  on  diffusion  temperature  and  time.  Figure  2  a,  b,  c  shows 
photoluminescence  spectra  from  sample  P4X  in  the  wavelength  range  of  800  nm  to  1700  nm 
grown  under  condition  shown  in  Table  I.  The  reference  spectrum  for  sample  as  grown  is  shown 
in  figure  2a.  The  dominant  is  the  RE  emission,  the  near  band  edge  bound  exciton  peak  is  also 
seen.  The  PL  spectra  after  Li  diffusion  at  temperature  500°C,  and  646°C  during  6  hours,  are 
shown  in  Fig.  2  b,  and  c,  respectively.  Sample  diffused  at  500°C  (Fig.2b)  shows,  beside  the  RE 
emission,  new  broad  band  with  peak  at  1.02  eV,  and  much  weaker  band  edge  emission. 
Figure  2c  shows  the  drastic  change  in  spectrum  after  Li  diffusion  at  646°C.  Strong  DA  pan- 
emission  appears  at  1.465  eV,  and  the  Cu-Li  peak  at  1.41  eV  [5]  shows  up.  The  strong  band 
at  1.36  eV  with  the  shoulder  at  1.34  eV,  and  1.02  eV  band  grows  strongly.  In  Fig.  3  a,  b,  c  we 
show  the  detail  of  rare  earth  spectra  for  above  sample  recorded  in  the  spectral  range  of  1 .5  pm 
to  1.6  pm.  Spectra  on  figure  3a,  b,  for  the  sample  as  grown  and  Li  diffused  at  500°C  are  almost 
identical.  The  erbium  spectrum  of  sample  doped  with  Li  at  646°C  (Fig. 3c)  shows  the  different 
relative  intensities  of  the  emission  lines  in  comparison  with  the  first  two  samples.  Figure  4  and 
Fig.  5  show  the  erbium  spectra  from  samples  P52  with  grown  condition  described  in  Table  I. 
In  Fig.4a  the  PL  spectrum  of  a  reference  sample  (as  grown)  is  shown,  b  and  c  are  the  spectra 
after  diffusion  of  Li  at  500  °C  and  570  °C  during  6  hours,  respectively.  The  only  differences 
between  these  samples  are  different  relative  intensities  of  the  emission  peaks,  as  seen  on  Fig. 
4d  where  we  overlap  the  spectrum  (a)  and  (c).  The  difference  reflects  the  Li  doping  of  the 
GaAs:Er  epilayers  at  different  temperatures.  PL  spectra  shown  on  Fig.5  a,  b  are  recorded  for 
samples  P52(4)  (a)  and  P52  (3).  Sample  P52  (4)  was  diffused  with  Li  twice  at  646°C  for  6h 
and  additionally  at  802°C  for  1  hour  (Table  II),  sample  P52  (3)  as  grown.  The  diffusion  at 
802°C  reduced  the  intensity  of  Er  emission  and  broadened  the  spectrum. 
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Fig.  1  Examples  of  SIMS  bulk  analysis  of  two  epilayers  grown  by  MOC  VD  and  diffused  with  Li  (concentration  of 
4  x  101*  cm*3),  (a)  sample  P5F  and  (b)  P52(2)  ( see  tables] 
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Fig.  2  The  photoluminescence  spectra  of  P4X  samples  recorded  at  2  K  in  the  range  800  -  1700  nm  using  He-Ne 
laser  (10  raW).  (a)  is  the  reference  spectrum  of  GaAs:Er  as  grown.  In  (b)  and  (c)  the  samples  are  diffused  with  Li 
at  temperature  500  °C  and  646  °C  for  6  hours,  respectively. 
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Fig.  3  The  photoluminescence  spectra  of  P4X  samples  recorded  at  2  K  in  the  range  1500  nm  -  1 600  nm  using  He- 
Ne  laser  (10  mW).  (a)  The  reference  spectrum  of  GaAs:Er  as  grown.  In  (b)  and  (c)  the  samples  are  diffused  with  Li 
at  temperature  500  °C  and  646  °C  for  6  hours,  respectively. 


is  the  overlapping  of  spectra  (a)  and  (c). 
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Fig.  5  The  photoluminescence  spectra  of  P52  samples  recorded  at  2  K  in  the  range  1500  nm  -  1600  nm  using  Ar+ 
laser  at  514.5  nm  (20  mW).  In  (a)  sample  P52(4)  was  diffused  with  Li  twice  at  646  oC  for  6  hours  and  additionally 
at  802  oC  for  1  hour,  and  (b)  sample  P52(3)  as  grown. 
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D.  Summary  and  Conclusion 

Photoluminescence  spectra  are  reported  in  GaAs  doped  with  Er  and  Li.  A  PL  bands  peaking  at 
1.41  eV  involved  a  0. 1  leV  Cuq,  -Li;  acceptor  associate  [5].  Similarly,  a  PL  band  peaking  at 
1.02  eV  is  ascribed  to  the  0.45  eV  acceptor  [10].  The  strong  band  at  1.36eV  is  related  to  the 
0.15  eV  Cu  acceptor  [5],  the  shoulder  at  1.34  eV  is  probably  related  to  a  double  Li^  acceptor 
[6].  In  the  samples  used  in  the  study  the  Cu  concentration  is  limited  to  accidental  contamination 
during  the  heat  treatment.  The  photoluminescence  of  Er  co-doped  with  Li  showed  different 
changes  in  erbium  spectrum  for  sample  P4X  and  P52.  The  erbium  spectrum  of  sample  P4X 
doped  with  Li  at  646^  in  Fig.3c,  shows  the  different  relative  intensities  of  the  emission  lines 
in  comparison  with  the  samples  as  grown  and  doped  with  Li  at  500°C  Fig. 3a  and  b  respectively. 
The  line  at  0.806  eV  (1537.5  nm)  increase  twice  in  intensity  comparing  with  as  grown  sample 
the  other  lines  showed  smaller  increase  in  intensity.  Spectra  in  Fig.4  a,-d,  are  for  sample  P52 
grown  at  500°C  and  doped  with  Li  at  different  temperature.  In  figure  4d  we  overlap  the 
spectrum  (a)  and  (c)  to  shown  the  differences  between  these  samples  as  grown  and  doped  with 
Li  at  570°C  during  6  hours.  PL  spectra  shown  on  Fig.5  a,  b  are  recorded  for  samples  P52(4) 
(a)  and  P52  (3).  Sample  P52  (4)  was  diffused  with  Li  twice  at  646°C  for  6  hours  and 
additionally  at  802°C  for  1  hour  (Table  II),  sample  P52  (3)  as  grown.  The  diffusion  at  802°C 
reduced  the  intensity  of  Er  emission  and  broadened  the  spectrum.  In  conclusion,  we  can  say  that 
the  Li  co-doping  with  Er  in  GaAs  causes  the  difference  in  relative  intensities  of  the  emission 
peaks  of  Er.  The  biggest  changes  in  spectra  are  observed  for  sample  doped  with  Li  at 
temperatures  range  570°C  to  650°C.  More  investigations  must  be  done  to  understand  the  Li 
involvement  in  Er  luminescent  center. 
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College  of  Engineering  and  Technology,  Ohio  University,  Athens,  Ohio  45701-2979 

Abstract 

The  ac  electroluminescence  of  thulium  doped  ZnS  embedded  in  boric  acid  matrix  was 
investigated.  The  El  emission  spectra  were  recorded  as  a  function  of  voltage,  frequency  and 
temperature.  The  emission  spectra  at  all  temperatures  show  only  strong  sharp  emission  lines 
which  were  assigned  to  transitions  within  the  4f  shell  of  Tm3+  ion.  The  emission  spectra 
consisted  of  five  groups  of  lines  centered  around  483  nm  ( ‘G4  -  3H$),  661.7  nm  (‘G4  -  3H4), 
710.5  nm  (^j  -  3H6),  785  nm  (*G4  -  3H5)  and  805.5  nm  (3F4  - 3  HJ.  The  emission  intensities 
(  I  )  of  different  peaks  where  recorded  as  a  function  of  the  applied  voltage, V.  The  plot  of 
natural  logarithm  of  /emission  intensity  (I)  versus  V'm  shows  a  straight  line  characteristic  over 
many  orders  of  emission  intensity  which  indicates  that  the  direct  impact  excitation  mechanism 
is  a  dominant  process.  The  experimental  values  of  the  483  emission  line  (transition  *G4  -  3H6) 
intensity  versus  the  applied  voltage  were  reproduced  theoretically  on  the  basis  of  the  direct 
impact  excitation  mechanism. 
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Introduction 

In  recent  years,  there  has  been  growing  interest  in  rare  earth  doped  ZnS  due  to  their 
promising  potential  in  obtaining  thin  film  electroluminescent  display s[l].  In  particular 
several  studies  have  been  reported  on  the  luminescence  of  Tm -doped  ZnS.  Thulium  is  of 
particular  interest  among  all  rare  earth  ions  because  it  shows  the  most  efficient 
cathodoluminescence  (0.216  W/W)  in  the  blue  emission  region  around  478  nm,  and  the 
line  centered  around  800  nm  has  a  power  efficiency  of  0.59  W/W  as  reported  by  Shrader, 
et  al  [2].  This  means  that  if  the  electrons  can  be  accelerated  to  ballistic  energies,  more 
efficient  devices  can  be  fabricated  which  use  direct  impact  excitation  of  rare  earth  ions[3]. 
Tm3+  emission  was  first  observed  by  Rothschild  [4]  who  found  emission  lines  centered  at 
about  478  nm.  Ibuki  and  Langer  [5,6]  reported  four  groups  of  emission  lines  around  365 
nm,  480  nm,  650  nm  and  800  nm.  They  investigated  the  optical  absorption  and  emission 
spectra  of  the  ZnS:Tm  powder  between  300  and  900  nm  at  different  temperatures. 
Zimmermann  and  Boyn  [71  applied  the  method  of  site  selective  excitation  spectroscopy 
to  Tm3+  ions  in  ZnS  crystals.  They  were  able  to  separate  three  different  sites  for  which 
4f-4f  emission,  excitation  and  absorption  spectra  were  determined  in  the  range  of  several 
4f  free-ion  transitions.  In  another  work  [81,  they  studied  the  excitation  and  thermal 
quenching  behavior  of  two  Tm3+  colters.  They  also  investigated  [9, 10]  the  excitation  and 
deexcitation  mechanisms  of  the  blue  emission  of  two  Tm3+  centers  and  demonstrated  that 
donor-acceptor  pairs,  involving  the  Tm  centers  as  the  donors,  play  an  important  role  in  the 
excitation  process. 

This  paper  reports  the  results  of  investigation  of  the  ac  electroluminescence  of 
Tm-doped  ZnS  as  a  function  of  voltage,  frequency  and  temperature.  Strong 
electroluminescence  emission  (at  all  spectral  regions)  was  observed  at  low  temperature  (8.5 
K)  as  well  as  at  room  temperature  (around  800  nm).  At  low  temperature,  the 
electroluminescence  spectra  of  ZnS:Tm  exhibit  a  fine  structure  that  consists  of  a  number 
of  groups  of  narrow  lines  with  a  halfwidth  of  about  0.5  nm.  By  investigating  the  voltage 
dependence  of  the  different  emission  lines,  it  has  been  found  that  the  direct  impact 
excitation  mechanism  is  the  dominant  mechanism  responsible  for  the  electroluminescence. 

Experimental 
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The  power  efficiency  of  the  ZnS  doped  with  Tm  polycrystalline  powder  at  478  nm 
determined  under  cathode-ray  excitation  is  0.216  WAV,  the  highest  yet  reported  for  any 
rare-earth-activated  phosphor  especially  in  view  of  the  very  narrow  emission  bandwidth 
(6  nm)  [2].  The  preparation  of  the  phosphors  can  be  found  elsewhere  [11].  To  obtain  an 
dectroluminophor  from  the  Tm-doped  ZnS  powder,  the  insulating  microcrystallites  were 
coated  with  CtijS  by  immersing  them  in  an  aqueous  solution  of  CuQ2  (0. 1  mol.  %)  at  room 
temperature  for  3  minutes  [12].  Then,  the  mixture  was  filtered  and  the  powder  was  dried 
at  room  temperature  (or  at  100°C).  Next,  the  Cu^ -coaled  ZnS:Tm  microcrystallites  were 
embedded  into  a  meta-boric  acid  (HBOj)  glass  matrix  by  heating  the  mixture  of 
ZnS:Tm:Cu2S  and  H3B03  in  weight  proportion  of  1:0.8,  [13].  The  mixture  was  heated 
to  a  temperature  in  the  range  of  170°C  -  300°C  (but  it  is  not  critical)  to  obtain  the 
meta-boric  add  glassy  medium  [13].  Boric  add  is  suitable  host  material  for  two  reasons: 
first,  it  provides  excellent  transparency  of  emission  in  the  visible  region  and  second  it  has 
a  low  mdting  point  (~169°C),  which  allows  the  powder  luminophor  to  be  embedded  in  a 
matrix,  with  a  very  small  possibility  that  chemical  reaction  or  induced  thermal  changes 
will  occur  in  the  ZnS:Tm  phosphore.  The  electroluminescent  cell,  shown  in  figure  1,  was 
prepared  by  pressing  the  glassy  mixture  at  a  temperature  above  the  melting  point,  between 
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Fig.  1  Structure  of  the  ZnS:Tm  EL  cell 
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Table  I  The  emission  intensity  ratios  of  the  main  peaks  to  the  peak  at  661.7  nm  at  10K. 
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two  transparent  conducting  glasses  (or  conducting  glass  and  stainless  steal  sheet).  Electrical 
contacts  to  the  conducting  surfaces  were  made  using  silver  paint.  The  thickness  of  the 
ZnS:Tm  layer  was  about  50  /tm  ( selected  from  a  large  number  of  samples)  and  the  area 
of  the  cell  was  4x4  mm2.  Some  investigated  samples  were  produced  using  the  teflon 
spacer  in  order  to  control  the  thickness  of  the  active  layer.  In  order  to  exclude  moisture, 
the  samples  were  kept  in  a  desiccator  with  a  powerful  drying  agent  ( the  samples  can  also 
be  protected  from  the  moisture  by  epoxy  encapsulation). 

A  sinusoidal  voltage  generator  was  used  as  the  excitation  source.  The  sample  was 
mounted  on  a  cold  finger  cooled  by  a  close-cycle  helium  cryostat  down  to  8.5  K.  The 
emission  was  dispersed  by  a  Jarrel  Ash  Model  78-490,  0.75  M,  scanning  monochrometor 
equipped  with  a  1180  grooves/mm  grating.  The  signal  was  detected  by  a  Hamamatsu 
R925  photomultiplier  with  a  spectral  response  in  the  range  of  400-930  nm.  The  total 
emission  was  detected  by  an  RCA  Model  7625  photomultiplier  .  The  signal  from  the 
photomultiplier  was  fitted  through  a  fast  ORTEC  474  preamplifier  to  a  single  photon 
counting  system,  and  then  acquired  by  the  ORTEC  7100  multichannel  analyzer.  The  data 
were  transferred  to  a  computer  for  further  manipulation  and  plotting. 

Results  and  Discussion 

The  electroluminescence  spectroscopic  investigations  of  ZnS:Tm  were  performed  in  the 
wavelength  range  between  400  nm  and  930  nm  at  different  temperatures  and  for  different 
applied  voltages  as  shown  in  figures  2  and  4,  respectively.  The  emission  spectra  show 
strong  sharp  emission  lines  of  Tm3+  in  the  visible  and  near  IR  regions  without  any 
background  emission  attributed  to  the  ZnS  host.  We  observed  five  groups  of  emission 
lines  centered  around  483  nm,  661.7  nm,  710.5  nm,  785  nm  and  805.5  nm  which 
correspond  to  the  transitions  'G4  -3H6,  ‘G4  -  3H4,  3F3  -  3H6,  *G4  -  3H5,  and  3F4  -  3H6, 
respectively  ( all  the  spectra  are  normalized  to  the  intensity  value  of  the  emission  peak  of 
661.7  nm  at  10  K,  see  Table  I  )  .  Figure  3  shows  the  energy  levels  scheme  of  Tm3+  ion 
assigned  after  [2,8].  Strong  sharp  emission  lines  have  been  observed  at  temperatures  in  the 
range  of  8.5  -  290  K.  In  figure  2  (A,  B,  C,  and  D),  the  (a)  curves  show  the  spectra  of 
ZnS:Tm  recorded  at  10  K  and  270  volts.  The  spectrum  in  the  region  475  -  500  nm  (figure 


ig.  2  AC  EL  spectra  of  Tm-doped  ZnS  for  different  temperatures  at  applied  voltage  270 
olts  and  frequency  1  kHz  .(A)  476  -  493  nm  ( transition  1G4~>  3H6),  (B)  659  -  678  nm 
transition  }G4~>  3H4 ),  (C)  700  -  715  nm  ( transition  3F3~»  3H6  ),  and  (D)  800  -  815 
m  ( transition  3F4~>  3Hg ). 


Fig.  3  Energy  levels  of  Tm3+  corresponding  to  the  emission  regions  of  (a)  476  -  493  nm 
(b)  659  -  678  nm,  (c)  700  -  715  nm,  (d)  782  -  786  nm  and  (e)  800  -  815  nm. 


Fig.  4  AC  EL  spectra  of  Tm-doped  ZnS  ( transition  1G4— »  3H6  )  for  different  applied 
voltages  at  temperature  10  K  and  frequency  1  kHz. 
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2A,  transition  'G4  -3H6 )  reveals  three  major  sharp  peaks  at  483  nm,  487  nm,  and  492  nm. 
The  emission  line  of  483  nm  is  narrower  and  about  three  times  stronger  than  the  other  two 
peaks.  The  emission  lines  at  650  -  690  nm  (figure  2B,  transition  ‘G4  -  3H4)  mainly  consists 
of  five  very  sharp  peaks  located  at  659.7  nm,  660.7  nm,  661.7  nm,  671.8  nm  and  678.2 
nm.  The  emission  line  at  661.7  nm  has  the  strongest  intensity  among  other  lines  in  this 
group  as  well  as  the  other  groups.  Two  small  peaks  were  observed  at  701.3  nm  and  710.5 
nm  (figure  2C,  transition  3F3  -  3H* ).  The  two  peaks  are  twenty  times  weaker  than  the 
emission  lines  at  661.7  nm  and  483  nm.  In  the  near  IR  region  at  800-820  nm  (  figure  2D, 
transition  3F4  -  3H6  )  the  emission  line  at  805.5  nm  was  four  times  weaker  than  the 
emission  intensity  of  the  661.7  nm  emission  line.  When  the  temperature  is  increased, 
another  emission  group  centered  at  785  nm  (transition  ‘G4  -  3H5)  appeared  at  130  K  and 
was  observed.  The  peaks  are  sharp  but  very  weak  and  disappeared  with  further  increase 
of  temperature  above  170  K.  Upon  increasing  the  temperature  up  to  room  temperature,  the 
emission  intensities  of  all  lines  in  the  483  nm  and  661.7  nm  regions.  For  the  other  two 
groups  centered  at  710.5  nm  and  805.5  nm,  the  emission  intensity  increased  with 
temperature  until  reaching  the  maximum  at  210  K.  The  emission  intensity  of  the  805.5  nm 
line  at  210  K  is  five  times  larger  than  the  emission  intensity  at  10  K.  With  further  increase 
of  temperature,  emission  intensities  of  all  lines  in  the  two  groups  decreased. 

The  Tm3+  in  ZnS  can  form  luminescence  centers  of  several  types,  depending  on 
its  surroundings.  All  these  centers  may  be  differentiated  by  the  symmetry  or  strength  of 
the  crystal  field.  On  the  other  hand,  this  field  defines  the  probabilities  of  transitions 
between  sublevels  of  Tm3+  terms.  Since  the  relative  concentration  of  various  types  of 
centers  depends  on  the  preparation  conditions,  the  different  ZnS:Tm  samples  can  differ 
strongly  in  the  relative  line  intensity  of  the  luminescence  spectrum.  It  was  report  by  [5, 
6a, b]  that  Tm3+  mainly  occupies  a  site  with  Td  symmetry. 

The  electroluminescence  intensity  of  ZnS:Tm  as  a  function  of  frequency  (  at  constant 
voltage)  was  recorded  in  the  frequency  range  1-16  kHz.  The  intensities  of  all  emission 
lines  increased  with  frequency  until  reaching  maximum  at  7  kHz  and  then  started 
decreasing  with  further  frequency  increase.  For  the  483  nm,  661.7  nm  and  805.5  nm 
emission  lines,  intensities  at  7  kHz  were  about  three  times  larger  than  at  1  and  16  kHz. 
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The  electroluminescence  spectra  were  recorded  at  different  applied  voltages  of  150  V, 
240  V,  320  V,  and  400  V.  The  intensities  of  all  the  Tm3+  emission  lines  increased  with 
increase  of  the  applied  voltage  as  shown  in  figure  4  for  the  range  475  -  495  nm  (transition 
'G4  ).  To  obtain  information  about  the  excitation  mechanism  of  electroluminescence 

in  ZnS:Tm,  we  have  investigated  the  intensities  of  different  emission  lines  as  a  function 
of  the  applied  voltage.  The  integrated  intensity  of  all  emission  lines  was  measured  as  a 
function  of  the  applied  peak  voltages  between  40  and  300  volts  and  for  different 
temperatures  in  the  range  of  10  -  290  K.  Figure  5  shows  a  plot  of  the  logarithm  of 
integrated  intensity,  I,  versus  the  applied  voltage,  V,  at  temperatures  10  K  and  290  K. 
The  emission  intensity,  I,  of  electroluminescence  due  to  the  impact  excitation 
mechanism  has  often  been  described,  following  Alfrey  and  Taylor  [14],  by 
I  =  I0exp(  -  bV',/2 ) 

where  I0  and  b  are  constants  and  V  is  the  applied  voltage.  The  plots  of  the  logarithm  of  the 
intensity,  I,  against  V'1/2  for  the  integrated  and  single-line  emissions  at  8.5  K  are  shown 
in  figure  6.  These  plots  show  that  the  emission  intensity  as  a  function  of  the  applied 
voltage  satisfies  the  above  equation  for  several  orders  of  magnitude.  This  proves  that 
ZnS:Tm  direct  impact  excitation  mechanism  by  hot  electrons  dominates  in  the 
electroluminescence  of  ZnS:Tm  crystals. 

For  more  quantitative  proof  of  the  direct  impact  excitation  mechanism  of 
electroluminescence,  we  performed  the  numerical  calculation  of  the  electroluminescence 
intensity  as  a  function  of  applied  voltage.  If  the  direct  impact  of  accelerated  electrons  with 
Tm3+  centers  is  the  excitation  mechanism  responsible  for  the  electroluminescence  in 
ZnS:Tm3+-Cu2S  system,  then  the  theoretical  model  proposed  by  [3,15,16,17]  gives  the 
expression  for  the  electroluminescence  intensity  of  the  Stark  multiplet  as: 


I=r\qnNffteJE)imo(z,zr)de 


f  (e,E)  =  £',+0  5  exp  (-be  ) 
a  =  (Ep  -  eEA.)  (2Ep  +  eEA) 1 
b'1  =  2/3  eEA.  +  l/3((eEA.)2/Ep) 


where 


Fig.  7  Comparison  of  Experimental  (star)  and  calculated  (solid  line)  Intensities  versus 

eEX.,  for  ZnS:Tm  emission  line  of  483  nm  ( transition  1G4— >  3H$). 
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Hq  is  the  quantum  efficiency,  n  is  the  density  of  hot  electrons,  N  is  the  density  of 
luminescence  centers,  f  (e,E)  is  the  Baraff  s  distribution  function  of  the  electron  energy, 
er  is  the  energy  of  the  r  multiple!  of  the  Tm3+  ion,  e  is  the  electron  charge,  E  is  the 
electric  field  in  ZnS,  and  X  is  the  electron  mean  free  path  .  Ep  is  the  optical  (LO)  phonon 
energy  (0.04  eV)  in  ZnS,  assuming  that  the  hot  electrons  are  scattered  by  the  optical 
phonons  of  the  ZnS  lattice. 

The  analytical  expression  for  the  inelastic  scattering  cross  sections  based  on  the  Bom 
approximation  considering  the  coulomb  interaction  between  hot  free  electron  and  the  4f 
electrons  of  the  rare  earth  ion  was  derived  by  Yu  and  Shen  [18].  The  impact  cross  section 
corresponding  to  the  transition  from  i  state  to  /  state  is  given  by 


18  n2me2h2c3  1 
»(« 2+2)2e3r  e 


[In 


v£+(e-er 

r~  T 

ye-(  *-*r) 


where  m  is  the  mass  of  the  electron,  c  is  the  speed  of  light  in  vacuum,  n  is  the  refractive 
index  of  ZnS,  e  is  the  energy  of  the  incident  electrons,  er  is  the  energy  corresponding  to 
the  transition  from  the  final  state  /to  the  initial  state  /,  g(/  is  the  degeneracy  of  the  initial 
and  final  states  and  Afi  is  the  electric  dipole  transition  probability. 

To  compare  with  the  experimental  results,  the  above  calculations  were  carried  out  for  the 
emission  line  of  483  nm  which  corresponds  to  transition  ‘G4  -  3H6.  The  threshold  energy 
of  this  transition,  er  is  2.57  eV,  Ep  is  0.04  eV  and  the  degeneracy  of  the  initial  and  final 
state  are  %f  =  9  and  g,  =  13  respectively.  The  refractive  index,  n,  of  ZnS  is  2.4  and  the 
electron  mass,  m,  is  the  free-electron  mass.  The  calculated  intensity  (solid  line)  with  the 
experimental  points  are  shown  in  figure  7.  The  mean  free  path  of  the  electron,  X,  was 
found  to  be  about  5.5  nm,  and  for  the  best  fit,  eEX  is  in  the  range  0.3  -  1.1  eV  (  which 
corresponds  to  the  electric  field  in  the  range  1.33  x  107  -  9.86  x  107  V/m  ).  The  computed 
intensity  increase  with  eEX  follows  the  experimental  data  well.  More  quantitative 
considerations  cannot  be  drawn  because  we  do  not  know  the  exact  value  of  the  electrical 
field  in  the  system,  the  concentration  of  the  RE3+  Centers  and  the  quantum  efficiency  nq- 
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Conclusion 

The  study  of  the  ac  electroluminescence  of  ZnS:Tm3+  embedded  in  boric  acid  matrix 
are  presented.  The  electroluminescence  spectra  at  low  and  room  temperatures  show  strong 
sharp  emission  lines  of  Tm3+,  in  the  visible  and  IR  regions  without  any  background 
emission  attributed  to  the  ZnS  host  semiconductor.  The  direct  impact  excitation 
mechanism  has  been  proved  to  be  the  dominant  mechanism  responsible  for  the 
electroluminescence  in  the  ZnS:Tm3+-Cu2S  system.  In  our  current  going  research1,  the 
powder  ZnS:Tm  luminophore  was  embedded  in  acrylic  binder  as  well  as  boric  acid.  Strong 
blue  electroluminescence  emission  was  observed  at  room  temperature  and  is  visible  in  an 
illuminated  room.  These  results  show  the  potential  for  a  blue  electroluminescence  large 
area  displays  made  from  powder  ZnS:Tm3+  luminophor. 

This  research  was  partially  supported  by  the  AFSOR  (  grant  #  90-0322  )  and  the 
Ohio  University  CMSS  program. 


1  The  results  will  be  published  elsewhere 
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In  this  work  we  have  developed  a  kinetics  model  of  energy  transfer  from  the  host  lattice  to  the  local¬ 
ized  core  excited  states  of  rare-earth  isoelectronic  structured  traps  (REI  traps).  The  presence  of  low- 
lying  empty  core  orbitals  in  rare-earth  impurities  introduces  new  excitation  and  recombination  phenom¬ 
ena.  To  adequately  describe  the  energy  transfer  to  a  REI  trap,  the  buildup  and  decay  kinetics  of  rare- 
earth  luminescence,  we  consider  six  separate  states  of  the  REI  impurity  (unoccupied,  electron  occupied, 
electron  occupied  excited,  exciton  occupied,  excited  electron  occupied,  and  excited  exciton  occupied). 
The  energy-transfer  processes  occur  through  an  Auger  mechanism  where  the  recombination  energy  of 
the  bound  electron  with  a  free  hole  is  transferred  nonradiatively  to  the  core  states,  or  energy  can  be 
transferred  from  the  bound  exciton  on  a  REI  trap  to  the  core  states.  If  the  initial  and  final  states  are  not 
resonant  (in  both  mechanisms),  the  energy  mismatch  must  be  accommodated  by  emission  or  absorption 
of  phonons.  Furthermore  we  discuss  details  of  several  quenching  processes,  which  are  incorporated  into 
the  kinetics  equations.  We  derive  two  sets  of  differential  equations  for  semi-insulating  and  n-type  semi¬ 
conductors  governing  the  kinetics  of  rare-earth  luminescence.  Equations  have  been  solved  by  a  numeri¬ 
cal  method  to  derive  the  time  dependence  of  the  rise  and  decay  kinetics  as  a  function  of  excitation  inten¬ 
sity.  The  numerically  simulated  luminescence  rise  and  decay  times  show  a  good  overall  quantitative 
agreement  with  experimental  data  obtained  for  InP:Yb,  over  a  wide  range  of  generation  rates. 


L  INTRODUCTION 

The  investigation  of  the  luminescence  properties  of 
rare-earth-doped  III-V  and  II-VI  semiconductors  is  of 
great  interest  both  from  the  scientific  and  application 
points  of  view.  The  scientific  interest  is  related  to  the 
uniqueness  of  optical  and  electrical  properties  of  rare- 
earth  impurities  in  semiconductor  hosts.  It  is  well  known 
that  rare-earth  luminescence  depends  only  slightly  on  the 
nature  of  the  host  and  on  the  temperature.  The  4 /  orbit¬ 
als  of  rare-earth  ions  incorporated  in  semiconductors  are 
so  deeply  buried  within  the  electronic  shell  that  the  ener¬ 
gy  levels  of  the  4/"  configuration  are  only  slightly  per¬ 
turbed  compared  to  free  ion  energy  levels.  The  electronic 
structure  of  the  rare-earth  luminescence  centers  and  their 
electrical  activities  as  well  as  their  indirect  photolumines¬ 
cence  and  electroluminescence  excitation  mechanisms, 
are  still  not  well  understood.  Among  the  rare-earth- 
doped  III-V  semiconductors,  InP:Yb  has  been  the  most 
extensively  studied. 1  ~9 

Ytterbium  in  InP  replaced  indium  on  a  substitutional 
site,10,11  and  acts  as  an  isoelectronic  trap.  It  was  original¬ 
ly  proposed  by  Whitney  el  al.iu)  and  confirmed  by 
others6,7,1213’  that  the  Yb  ion  creates  an  electron  trap  at 
30  meV  below  the  bottom  of  the  conduction  band.  Re¬ 
cently,  admittance  spectroscopy12'3’  was  used  to  identify 
the  electrical  activity  of  Yb  in  n-  and  p-type  InP.  It  was 
found  that  Yb  in  InP  creates  a  hole  trap  at  50  meV  above 
the  valence  band,  and  an  electron  trap  at  29  meV  below 
the  conduction  band.  The  50-meV  trap  may  be  related  tc 
other  impurities  unintentionally  incorporated  into  the 
crystal.  This  interesting  result  required  confirmation  us¬ 
ing  a  sample  grown  by  a  more  refined  crystal  growth 


technique  to  ensure  the  high  purity  of  the  crystal.  Re¬ 
cently  it  has  been  reported  that  Er  in  InP  [Ref.  12(b)]  and 
Yb  in  GaAs  [Refs.  12(c)  and  12(d)]  introduced  electron 
traps  at  60  and  63  meV,  respectively,  below  the  conduc¬ 
tion  band.  Colon  et  al.nu)  investigated  low-temperature 
photoluminescence,  selectively  exited  luminescence,  and 
deep-level  transient  spectroscopy  on  erbium-implanted 
GaAs.  Conclusions  of  these  measurements  is  that  Er  im¬ 
plantation  introduces  in  GaAs  two  hole  traps  at  84  and 
340  meV  above  the  valence  band.  Several 
authors51  *’’5<b)'6’8,9,l3(*’’B(b)  have  proposed  a  model  that 
involves  recombination  of  electron-hole  pairs  at  the  rare- 
earth  (RE)  traps  to  explain  the  excitation  of  RE  core 
states. 

In  this  paper  we  discuss  only  the  structured  isoelect¬ 
ronic  traps  in  III-V  semiconductors  introduced  by  triply 
charged  rare-earth  ions  replacing  the  element  from 
column  III  (or  another  more  complex  RE  isoelectronic 
traps).  Furthermore,  we  develop  the  luminescence  kinet¬ 
ic  models  that  describe  the  energy-transfer  and  recom¬ 
bination  processes.  The  presence  of  low-lying  empty  core 
orbitals  in  rare-earth  impurities  introduces  new  excita¬ 
tion  and  recombination  phenomena,  which  will  be  dis¬ 
cussed  in  detail.  The  RE  luminescence  rise  time  of  the 
rare-earth-doped  semiconductors,  excited  indirectly 
above  the  band  gap,  contain  information  about  the 
energy-transfer  processes  from  the  host  to  the  4/"  elec¬ 
tron  system.  It  is  shown  that  the  study  of  the  rise  time  at 
different  temperatures,  excitation  intensities,  and  excita¬ 
tion  pulse  durations  can  provide  important  information 
about  the  energy-transfer,  radiative,  and  nonradiative 
processes,  respectively.  The  numerically  simulated 
luminescence  rise  and  decay  profiles  show  a  good  quanti- 
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tative  agreement  with  experiment  over  a  wide  range  of 
generation  rates.  Finally,  the  possible  quenching  mecha¬ 
nisms  and  the  temperature  dependence  of  rare-earth 
luminescence  are  discussed. 

II.  THEORETICAL  FORMULATION 

It  is  well  known  that  isoelectronic  impurities  in  semi¬ 
conductors  produce  bound  states  in  the  forbidden  gap, 
binding  an  electron  or  a  hole. Mi  15  An  isoelectronic  center 
can  form  bound  states  because  of  a  short  range  central¬ 
cell  potential.  According  to  Thomas,15  the  primary  fac¬ 
tors  affecting  the  binding  potential  are  the  electronega¬ 
tivity  and  the  size  differences  between  the  impurity  and. 
the  host  ion  which  it  replaces.  It  is  found  experimentally 
that  only  very  large  atoms  or  very  small  atoms  produce 
isoelectronic  traps  because  they  create  large  lattice  dis¬ 
tortion  induced  by  the  substitution.  Thomas  and  co- 
workers16,15  have  pointed  out  that  to  create  a  large  bind¬ 
ing  potential,  the  substituted  atom  must  generate  a  no¬ 
ticeable  change  in  the  local  properties  of  the  lattice.  This 
is  likely  to  produce  an  unfavorable  free  energy  of  the 
solution,  and  hence  a  rather  low  solubility  can  be  ob¬ 
served.  For  instance,  the  maximum  concentration  seen 
for  bismuth  in  GaP  was  less  than  101B  cm'3. 17  Low  solu¬ 
bility  is  also  observed  for  rare  earths  in  III-V  semiconduc¬ 
tors.18 

Allen19  proposed  different  binding  mechanisms  for 
isoelectronic  traps.  According  to  Allen  the  isoelectronic 
impurity  potential  does  not  come  from  a  pseud  .' potential 
difference  of  two  isoelectronic  atoms.  Other  possible 
sources  are  the  spin-orbit  coupling,  and  the  strain  field  in 
the  close  vicinity  of  the  impurity  due  to  the  size 
difference  between  the  impurity  and  the  host  atom  which 
it  replaces.  The  main  results  of  this  theory  are  that  the 
perturbing  potential  at  an  isoelectronic  impurity  may  be 
attracted  simultaneously  to  both  the  electron  and  the 
holes,  so  bound  exciton  states  can  occur  without  bound 
single-particle  states.  Baldereschi  and  Hopfield20  have 
proposed  a  theory  of  isoelectronic  traps  assuming  that 
the  short-range  potential  arises  from  core  differences,  in¬ 
cluding  spin-orbit  interaction  between  the  dopant  atom 
and  the  host  atom  it  replaces^  The  relaxation  of  the  host 
crystal  around  the  impurity  as  well  as  the  screening  mod¬ 
el  considered  appear  to  be  important  for  the  binding  en¬ 
ergy.  However,  discrepancies  between  experimental 
binding  energies  and  those  calculated  from  the  differences 
in  the  atomic  pseudopotentials  are  observed.  Excellent 
reviews  of  existing  theories  and  experimental  data  of 
isoelectronics  impurities  were  given  by  Baldereschi21  and 
Dean22  and  Czaja.23  In  GaP  two  isoelectronic  traps  have 
been  extensively  investigated,  namely  a  nitrogen  electron 
trap,  and  a  bismuth  hole  trap  formed  by  substituting  Bi 
and  N  for  phosphorus.22  In  the  direct-band-gap  III-V 
semiconductors,  isoelectronic  impurities  have  been  inves¬ 
tigated  only  in  InP  doped  with  bismuth.24,25  The  neutral 
Bi,  in  a  P  site,  creates  a  hole  trap,  and  the  isoelectronic 
complex  (BivY)-(Bi,A')  bound-excitonic  molecules.25  In 
II-VI  semiconductors,  ZnTe  doped  with  oxygen-electron 
traps,  and  CdS  doped  with  tellurium  hole  traps,  were  in¬ 
vestigated  in  detail. 

It  is  notable  that  all  isoelectronic  impurities  discussed 


above  involve  substitutions  on  the  anion  sites.  The  previ¬ 
ous  investigation26  of  cationic  isoelectronic  substituent 
(Mg,  Ca,  Sr,  Ba)  in  ZnSe  and  ZnTe  found  no  evidence  of 
the  presence  of  isoelectronic  traps.  A  strong  photo¬ 
luminescence  (PL)  was  reported  recently  from  Mg-doped 
ZnSe  [Refs.  27(a)  and  27(b)]  and  ZnS  (Ref.  28)  which  may 
arise  from  a  new  isoelectronic  center  generated  by  mag¬ 
nesium.  We  observed27  sharp  emission  with  a  half-width 
of  16-60  meV  (for  Mg^Zn^^Se,  x  =0.04),  as  required 
for  exciton  transitions  at  low  temperatures.  Furthermore 
the  PL  characteristics  shift  dramatically  from  deeper  ex¬ 
trinsic  emission  in  ZnSe  to  be  dominated  by  narrow 
near-band-gap  emission  at  all  temperatures  in  the  range 
2-300  K  in  Mgo  ^Zn^  ^Se.  The  temperature  depen¬ 
dence  of  the  energy  position  of  this  peak  follows  the  ex¬ 
pected  behavior  for  a  free  exciton  transition  in  ZnSe. 
This  peak  remains  strong  in  the  PL  spectrum  all  the  way 
up  to  room  temperature.  This  kind  of  behavior  is  not 
usually  observed  for  bound  exciton  transitions  related  to 
neutral  donor  or  acceptor  centers.  However,  this 
behavior  is  typical  for  excitons  bound  to  isoelectronic 
traps  for  which  the  dominant  decay  mode  is  radiative. 
The  presence  of  the  third  electronic  particle  in  the  first 
two  cases,  but  not  in  the  third,  was  shown  to  introduce  a 
dominant  Auger  decay  mode  for  the  bound  exciton.  In 
the  Auger  quenching,  all  the  recombination  energy  is 
transferred  to  this  third  electronic  particle  that  is  ex¬ 
pelled  deep  into  the  conduction  (electron)  or  valence 
band  (hole),  respectively. 

The  striking  feature  of  excitons  bound  to  isoelectronic 
traps  is  a  long  luminescence  decay  time,  ranging  from  a 
few  hundred  to  a  few  thousand  nanoseconds.25,29'33  The 
lifetimes  of  a  neutral  donor  or  acceptor  bound  excitons  in 
direct-gap  semiconductors  are  in  the  range  of 
nanoseconds.  For  example  the  decay  time  of  excitons 
bound  to  neutral  donors  or  acceptors  in  InP  are  0.5  and 
1.5  ns,  respectively. 34<al  The  lifetime  of  an  exciton  bound 
to  a  neutral  donor  in  GaAs  is  1.07  ns.34<bl  In  contrast  the 
lifetime  of  an  exciton  bound  to  a  Bi  isoelectronic  trap  in 
InP  is  about  200  ns.25 

Table  I  shows  that  the  outer  electron  configurations  of 
RE3+  ions  are  the  same  (5s25p6).  If  the  rare-earth  ions 
replace  the  element  from  column  III  in  III-V  compounds 
that  are  isovaient  concerning  outer  electrons  of  RE3  + 
ions  (see  Table  I),  they  create  isoelectronic  traps  in  III-V 
semiconductors.  This  does  not  require  association  with 
other  near  distant  charge  compensating  lattice  defects  or 
impurities,  as  is  so  common  in  II-VI  semiconductors. 
The  above  conclusion  is  supported  by  the  fact  that  the 
atomic  covalent  radii  (ionic  RE3T  for  all  rare  earths  are 
bigger  than  atomic  radii  of  Ga  and  In  that  they  are  re¬ 
placing.  Pauling's  electronegativity  of  rare-earth  ele¬ 
ments  is  in  the  range  of  1.1-1.25,  and  is  smaller  than  Ga 
(1.81)  and  In  (1.78)  for  which  it  substitutes.  We  known 
from  different  investigations  that  Yb  substituted  for  In  in 
InP  behaves  according  to  the  above  experimental  rule 
and  creates  an  isoelectronic  trap.  We  have  evidence  that 
the  other  RE  ions  in  1I1-V  semiconductors  can  occupy 
different  sites  (not  only  substitutional).  The  rare-earth 
isoelectronic  trap  must  not  necessarily  be  the  “pure”  sub¬ 
stitutional  center,  if  the  rare-earth  ions  are  very  active 
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TABLE  1.  Electron  configuration  of  RE  atoms.  RE3*  ions  (and  some  elements),  ionic,  covalent  radii,  and  electronegativity. 

Ionic 


Element 

Electron 

configuration 

Electron 
configuration 
REJ+  4  others 

radius  (A) 
charge 

2+  3+ 

Covalent 

radius 

(A) 

Electronegativity 
a:  (Pauling's) 

Cerium 

4/,5rJ5p46s2 

4/'5s25p4 

1.02 

1.65 

1.12* 

Praseodymium 

4/35s2Sp‘6s2 

4/25r25p4 

1.00 

1.65 

1.13* 

Neodymium 

4/45r25p46s2 

4/35r25p4 

0.99 

1.64 

1.14* 

Promethium 

4/35r25p46r2 

4/45r25p4 

0.98 

1.63 

1.13* 

Samarium 

4/‘5i25p‘6s2 

4/35r25p4 

0.97 

1.62 

1.17* 

Europium 

4/75i25p‘6s2 

4/45s25p4 

0.97 

1.85 

1.20* 

Gadolinium 

4/75r15p45d,6s2 

4/75s25p4 

0.97 

1.61 

1.20* 

Terbium 

4/’5ii5p*6s2 

4/*5r25p‘ 

1.00 

1.59 

1.20* 

Dysprosium 

4/'°5s25pt6s2 

■4/95r25p‘ 

0.99 

1.59 

1.22* 

Holmium 

4/"5r15p‘6rJ 

4/'°5r25p4 

0.97 

1.58 

1.23* 

Erbium 

4/,25s25p46s2 

4/H5rl5p4 

0.96 

1.57 

1.24* 

Thulium 

4/'35s25p46r2 

4/,25rJ5p4 

0.95 

1.56 

1.25* 

Ytterbium 

4/'45s25p46s2 

4/l35s25p4 

0.94 

1.74 

1.10* 

Gallium 

3d,04s24p 

0.62 

1.26 

1.13b 

1.81* 

Indium 

4d,05sI5p 

0.81 

1.44 

0.99b 

1.78* 

Zinc 

3d'°4s2 

0.74 

1.23 

0.91b 

1.65* 

Cadmium 

4d,0Ss2 

0.97 

1.48 

0.83k 

1.69* 

Mercury 

Sdw6s2 

1.10 

1.49 

0.79b 

2.00* 

‘Elemental  electronegativity  in  tetrahedrally  coordinated  environments  (Ref.  35). 
‘Reference  36. 


chemically,  they  can  create  a  more  complex  center  in¬ 
volving  other  imparity  or  native  defects.  The  recent  ex¬ 
perimental  data  discussed  in  Sec.  I  shows  that  RE  ions  in¬ 
troduce  electron  or  hole  traps  in  III-V  semiconductors, 
and  we  do  not  have  evidence  that  RE  ions  act  as  donors 
or  acceptors.  The  important  roles  of  oxygen  on  RE 
luminescence  have  been  discussed  recently.  Clearly  we 
need  more  experimental  and  theoretical  investigations  de¬ 
voted  to  rare-earth  impurities  to  gain  knowledge  about 
the  electrical  activities  of  RE  ions  in  semiconductors. 

The  rare-earth  isovalent  traps  that  we  call  isoelectronic 
“structured”  impurities17  possess  unfilled  4/*  core  shells. 
The  luminescence  structure  arises  from  intraconfigura- 
tional  /-/  transitions  in  the  core  of  the  isoelectronic 
“structured”  impurities.  The  presence  of  low-lying  emp¬ 
ty  core  orbitals  in  rare-earth  impurities  introduces  new 
excitation  and  recombination  phenomena  (which  will  be 
discussed  in  detail  below).  It  distinguishes  these  impuri¬ 
ties  from  the  “simple”  impurities  (from  the  main  group  of 
elements  of  the  Periodic  Table).  The  “simple”  impurity 
typically  introduces  only  effective-mass-like  states  in  the 
forbidden  gap  of  the  host  crystals.  According  to  Robbins 
and  Dean,17  the  formation  of  a  bound  state  at  the  struc¬ 
tured  cationic  isoelectronic  impurities  is  fairly  common, 
which  is  in  contrast  to  the  situation  normally  found  for 
the  anionic  substituent  in  semiconductors  discussed 
above. 

The  isoelectronic  trap  can  be  attractive  either  for  elec¬ 
trons  or  for  holes,  and  according  to  Allen’s  theory  it  can 
bind  the  exciton  as  a  single  entity.  Since  there  is  no 
charge  involved,  the  isoelectronic  center  forms  the  bound 
states  by  a  short-range  central-cell  potential.  It  is  gen¬ 
erally  accepted  that  the  formation  of  bound  states  is  a 
specific  property  of  the  impurity  and  lattice  combinations 


discussed  above.  After  an  isoelectronic  trap  has  captured 
an  electron  or  a  hole,  the  isoelectronic  trap  is  negatively 
or  positively  charged,  and  by  Coulomb  interaction  it  will 
capture  a  carrier  of  the  opposite  charge,  creating  a  bound 
exciton. 

It  has  been  well  established  that  the  “simple”  isoelect¬ 
ronic  traps  can  act  as  very  efficient  centers  for  radiative 
recombination  in  semiconductors.  The  structured 
isoelectronic  cationic  substitutional  impurities  (trivalent 
RE3+)  in  many  phosphors  lead  to  efficient  characteristic 
luminescence.17  The  trivalent  rare-earth  ions  also  create 
structured  substitutional  isoelectronic  traps  (REI  trap)  in 
some  III-V  semiconductors. 

The  “simple”  isoelectronic  center  in  III-V  materials 
can  exist  in  three  possible  states  (to  be  identified  later),  in¬ 
stead  of  two  as  in  the  case  of  the  Shockley-Read-Hall 
(SRH)  recombination  model.32  In  the  case  of  rare-earth 
isoelectronic  traps  the  kinetics  model  is  even  more  com¬ 
plicated  because  of  energy-transfer  processes  between  the 
localized  state  in  the  forbidden  gap  of  the  host,  and  the 
localized  core  states  of  structured  isoelectronic  impuri¬ 
ties.  There  are  three  possible  mechanisms  of  energy 
transfer.  The  first  is  the  energy-transfer  process  from  ex- 
citons  bound  to  “structured”  isoelectronic  centers  to  the 
core  electrons.  It  takes  place  as  a  result  of  the  electro¬ 
static  perturbation  between  the  core  electrons  of  the 
“structured”  impurity  and  the  exciton  effective-mass-like 
particle.17  This  model  is  a  modification  of  the  Shaffer- 
Williams  model  of  intrapair  energy  transfer  to  “struc¬ 
tured”  isoelectronic  traps.37  The  second  mechanism  is 
transfer  of  energy  to  the  core  electrons  involving  the 
“structured”  isoelectronic  trap  occupied  by  electron 
(hole)  and  free  hole  (electron)  in  the  valence  (conduction) 
band.  The  third  mechanism  is  the  transfer  through  an  in- 
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elastic  scattering  process  in  which  the  energy  of  a  free  ex- 
citon  near  a  “structured”  trap  is  given  to  the  localized 
core  excited  states.17  If  the  initial  and  final  states  are  not 
resonant,  the  energy  mismatch  must  be  distributed  in 
some  way,  e.g.,  by  phonon  emission  or  absorption. 
According  to  Robbins  and  Dean,*7  if  the  atomic  core  ex¬ 
citations  are  strongly  coupled  to  the  host  phonons,  the 
energy-transfer  probability  is  likely  to  be  higher.  Strong 
phonon  coupling  may  also  be  desirable  in  ensuring  that 
relaxation  down  the  ladder  of  core  excited  states  occurs 
quickly,  thus  preventing  back  transfer.  However,  for 
efficient  radiative  recombination,  the  phonon  coupling 
should  not  be  strong,  in  order  to  prevent  core  deexcita¬ 
tion  by  nonradiative  multiphonon  emission.  In  this  re¬ 
gard  the  rare-earth  “structured”  impurity  seems  to  be 
ideal. 


III.  KINETIC  MODELS 

The  SRH  model  that  allows  for  only  two  possible  states 
fails  to  give  the  correct  general  description  of  recombina¬ 
tion  kinetics  of  isoelec tronic  traps.32  The  above  discus¬ 
sion  shows  that  isoelectronic  impurities,  both  “simple”  or 
“structured,”  can  act  as  efficient  centers  for  radiative 
recombination.  The  “simple”  traps  can  exist  in  three 
possible  states:  (1)  empty,  (2)  electron  (hole)  occupied, 
and  (3)  exciton  occupied.  In  the  case  of  RE3+  “struc¬ 
tured”  isoelectronic  centers,  the  model  is  more  compli¬ 
cated  because  the  center  can  exist  in  six  possible  states. 
Furthermore  the  energy-transfer  processes  between  the 
localized  states  in  the  forbidden  gap  and  core  states  com¬ 
plicate  the  model.  In  our  model,  we  assume  that  the 
isoelectronic  trap  is  an  electron  trap  such  as  YbJ+  in 
semi-insulating  (SI)  InP.  When  isoelectronic  traps  are 
present  in  n-  or  p-type  materials  the  model  will  be 
different,  especially  at  high  temperature,  and  must  be 
modified  separately  for  n-  and  p-type  material.  The 
asymmetry  between  n-  and  p- type  semiconductors  results 
from  the  fact  that  the  isoelectronic  center  binds  only  an 
electron  (hole).  Thus  in  an  n-type  material,  a  fraction  of 
the  isoelectronic  traps  will  be  occupied  by  electrons  even 
before  the  sample  is  excited,  while  in  p- type  material  all 
the  isoelectronic  electron  centers  will  be  empty.39’40  At 
low  temperature  HK)  there  will  be  no  difference  be¬ 
tween  the  n,  p,  or  SI  samples  because  the  electrons  or 
holes  will  be  frozen.  The  donor  and  acceptor  will  act  as 
“shunt”  recombination  centers.  Figures  1-5  show  the 
symbols’  definitions,  and  physical  models  used  in  the 
kinetic  analysis  of  the  energy-transfer  processes  and 
recombination  involving  “structured”  isoelectronic  im¬ 
purities.  The  proposed  model  accounts  for  energy 
transfer  from  the  host  to  the  core  states  through  localized 
states  in  the  forbidden  gap,  and  for  the  dependence  of  the 
rise  time  on  excitation  intensities  and  temperature,  in¬ 
cluding  nonradiative  recombination  centers.  In  cases  of 
n-  or  p-type  semiconductors,  we  incorporate  donor  or  ac¬ 
ceptor  centers  to  the  model  [see  Figs.  1(d)  and  1(e)], 

By  fitting  the  calculation  to  the  experimental  data,  we 
can  estimate  important  parameters  related  to  energy 
transfer  from  the  lattice  to  RE3  +  centers,  the  Auger  pro¬ 
cesses,  temperature  quenching  mechanisms,  and  other 


FIG.  1.  Schematic  diagram  of  recombination,  trapping 
(thermalization),  and  exciton  formation  and  dissociation  pro¬ 
cesses  involving  structured  rare-earth  isodectronic  electron 
traps  (represented  by  ellipse)  with  the  atomidike  4/"  core 
states,  (a)  Trapping  (liberation)  )  of  electron  on  RE! 
traps,  (b)  Formation  of  exciton  (shadow  ellipse)  by  hole  capture 
(r^ );  liberation  of  hole  (rv);  dissociation  of  exciton  (tJS  on  the 
REI  trap,  (c)  Shunt  electron  (hole)  trapping  time  ra  (r^).  (d) 
and  (e)  The  trapping  times  for  electron,  and  hole,  rfA,  by  a 
donor  (ionized),  and  an  acceptor  (neutral),  respectively,  and 
bound  electron  (hole)  to  free  holes  (electrons)  recombination 
rates  are  rnD  —BnD(N%)p  [rnA  —BnA(N*A  )it]. 


important  parameters  that  we  discuss  below.  We  take 
the  depth  of  the  Yb3+  REI  trap  (electron  isoelectronic 
trap)  from  the  experiment  to  be  30  meV  below  the  con¬ 
duction  band,  because  we  wish  to  apply  results  of  the  cal¬ 
culations  to  InP:Yb.  This  trap  may  exist  in  six  distinct 
states:  (1)  the  neutral  unoccupied  trap  (concentration 
N0),  (2)  the  negatively  charged  (concentration  TV_ ),  (3) 
the  exciton  occupied  (concentration  Nx),  (4)  the  neutral 
excited  (concentration  Nq),  (S)  the  excited  negatively 
charged  (concentration  TV* ),  and  (6)  the  excited  exciton 
occupied  (concentration  Nx).  The  total  concentration  of 
isoelectronic  traps  N  is  given  by 


FIG.  2.  A  schematic  representation  of  energy-transfer  pro¬ 
cesses  to  core  states  of  structured  impurity  and  radiative  and 
nonradiative  transitions,  (a)  Auger  energy-transfer  process 
where  the  recombination  energy  of  the  bound  electron  with  the 
free  hole  is  transferred  to  the  REI  impurity  core  states 
(coefficient  Br).  is  the  nonradiative  Auger  process 

coefficient  involving  bound  electron  and  two  free  holes,  BFBt  the 
coefficient  of  radiative  recombination  of  bound  electron  to  free 
hole,  and  r„  the  recombination  of  bound  electron  with  a  hole 
on  a  distant  acceptor,  (b)  The  energy-transfer  process  from  a 
bound  exciton  on  a  REI  trap  (shadow  ellipse)  to  core  states  ( BXt 
transfer  coefficient)  and  and  rx„  are  the  bound  exciton  radia¬ 
tive  and  nonradiative  recombination  times,  respectively 
( 1  /r2=  1  /rlr  +  1  /tx„  ). 
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a  b 


FIG.  3.  Energy  back  transfer  processes,  (a)  The  auto  deioni¬ 
zation  Auger  quenching  mechanism  and  its  characteristic  time 
rn  (Auger  transfer  rate  rm—BmNT.  >•  (b)  The  energy  back 
transfer  process  from  the  excited  REI  trap  to  a  bound  exciton 
with  characteristic  time  rtXT  (Auger  transfer  rate  rg xp 
=Bbxt  AT/). 


a  be  d 


FIG.  3.  Auger  quenching  processes  involving  interaction  of 
the  core  excited  REI  trap  with  (a)  and  (b)  free  electrons  or  holes 
and  their  coefficients  2?rn  and  Bn4.  Interaction  of  the  core  ex¬ 
ited  REI  trap  with  an  electron  (c)  [hole  (d)]  on  neutral  donor 
(acceptor)  and  their  coefficients  Bgn  and  2?rrt- 


N=N0+N~+Nx+NS+Nl  +Nj  . 

Many  of  the  symbols  are  the  same  as  those  defined  in 
Refs.  31-33.  Figure  1(a)  shows  how  the  neutral  isoelec- 
tronic  trap  is  transformed  into  an  JV_  state:  an  electron 
is  captured  within  a  time  defined  by  r<,=(uth<TREW)-1, 
where  is  the  thermal  velocity  of  the  free  electron,  and 
0-,^  is  the  cross  section  for  the  capture  of  electrons  by 
unoccupied  traps.  The  capture  rate  of  electrons  by  the 
isoelectronic  trap  is  (n /Ta)(NQ/N),  where  n  is  the  free 
electron  density.  The  N_  center  can  be  transformed  to 
the  neutral  rare-earth  core  excited  center  through  an 
Auger  process  where  the  recombination  energy  of  the 
bound  electron  with  free  hole  is  transferred  nonradiative- 
ly  to  the  core  states  with  a  rate 

rT—BTN-p  , 

where  BT  is  the  energy-transfer  coefficient,  and  p  is  the 
free  hole  density.  If  the  initial  and  final  states  are  not  res¬ 
onant,  the  energy  mismatch  must  be  distributed  in  some 
way  such  as  phonon  emission  or  absorption  [Figs.  2(a) 
and  2(b)],  The  N _  may  also  recombine  radiatively  or 
nonradiatively  (via  the  Auger  process)  with  a  hole  in  the 


a  b 


FIG.  4.  Auger  nonradiative  recombination  involving  the  in¬ 
teraction  of  the  core  excited  REI  trap  (Nq)  with  an  electron 
trapped  on  separate  centers  N~  (a)  and  N*  (b),  where  Bbtx  and 
BBT2  are  the  Auger  coefficients. 


valence  band  with  transition  rate  W-  or 

NL  )p  and  or  JV*  )p\  respectively  [Fig. 

2(a)],  or  recombine  radiatively  with  a  hole  trapped  on  a 
distant  acceptor  (in  the  case  of  p-type  material).  The  last 
three  processes  transform  the  N  _  center  into  the  N0  neu¬ 
tral  trap.  The  Nq  may  capture  an  electron  and  be 
transformed  into  an  AT*  center.  The  Af*  may  lose  the 
electron  through  processes  described  above,  or  by  ion  au¬ 
todeionization  Auger  process  (IV 1  -♦Af0+e-kinetic  ener¬ 
gy).  This  is  the  Auger  nonradiative  luminescence 
quenching  mechanism  of  the  excited  rare-earth  isoelect¬ 
ronic  trap,  with  characteristic  time  as  shown  in  Rg. 
3(a).  The  exciton  bound  to  a  REI  trap  with  exited  core 
states  (Nx),  can  back  transfer  energy  from  core  excita¬ 
tions  to  the  exciton,  and  dissociate  it  with  liberation  of  an 
electron  or  hole.  Once  formed,  the  AT_  and  /V*  states 
can  also  be  converted  into  N0  and  1VJ  states,  respectively, 
by  thermal  ionization  (time  tu),  or  to  Nx  and  Nx,  respec¬ 
tively,  by  capturing  a  hole  from  the  valence  band  with 
time  defined  by  =(t)tho/uAr)_1  [see  Fig.  1(b)]. 

The  second  important  energy-transfer  process  from  an 
exciton  bound  to  a  REI  trap  to  core  states  is  shown  in 
Fig.  2(b).  In  this  process  energy  is  transferred  from  the 
bound  exciton  (on  isoelectronic  trap  Nx)  to  the  core 
states  with  the  rate  rXT=BXTNx,  where  Bxj  is  the 
transfer  coefficient,  and  is  the  characteristic  time.  If 
the  initial  and  final  states  are  not  resonant,  the  energy 
mismatch  must  be  accommodated  by  emission  or  absorp¬ 
tion  of  phonons.  At  sufficiently  high  temperatures,  the 
electron  from  iV_  and  Af*  may  be  thermally  emitted  to 
the  conduction  band  at  a  rate  (n,  /ru)  (N  _  /N  or  Af  *  /AO, 
where  n,=NcP,exp[—  £RET//c7],  (1//?,)  is  the  degenera¬ 
cy  of  the  structured  isoelectronic  trap  level  at  energy 
£ret  below  the  conduction  band,  and  Nc  is  the 
conduction-band  density  of  states. 

The  exciton  bound  to  isoelectronic  traps  Nx  and  Nx 
can  thermally  dissociate  by  either  of  two  processes:  (1)  it 
can  dissociate  into  a  free  exciton  UO  and  neutral  N0,  or 
an  neutral  excited  (Nq)  REI  trap,  respectively;  or  (2)  it 
can  dissociate  with  the  liberation  of  a  hole,  with  the  time 
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1  /tv  =<1/t,x  ){pk  /N ),  where  pk  =  NeQk  *exp[  - Ek  / 
kT),  N„  is  the  effective  density  of  state  in  the  valence 
band,  and  /3*  is  the  degeneracy  of  the  hole  component  of 
the  exciton  having  binding  energy  Ek.12'*9 

The  REI  traps  NB,  N*_,  and  Nx  may  also  be  deexcited 
to  N0,  N_,  and  Nx,  respectively,  by  desirable  radiative 
transitions  between  the  4/"  core  states,  emitting  photons 
with  a  decay  time  r3. 

The  nonradiative  decay  channel  of  free  carriers 
through  an  additional  type  of  trap  Nt,  the  “shunt  path” 
with  two  states,  is  shown  in  Fig.  1(c).  The  concentrations 
of  traps  in  these  two  states  are  denoted  Nf,  the  neutral 
unoccupied,  and  N~,  the  negatively  charged  (bound  elec¬ 
tron),  respectively  (Ns  =N°+N~).  The  probability  per 
unit  time  for  an  electron  (hole)  to  be  captured  by  an  emp¬ 
ty  (electron  occupied)  shunt  is  l/r„  and  (1/r^),  respec¬ 
tively.  For  intrinsic  excitation  (band  to  band)  the  shunt 
paths  deplete  the  express  electron-hole  pairs.  Under  such 
excitation  the  shunt  paths  can  be  either  partially  or  com¬ 
pletely  saturated,  depending  on  the  concentration,  the 
time  constants,  and  die  generation  rate  of  electron-hole 
pairs.  Figures  1(d)  and  1(e)  show  the  characteristic  pa¬ 
rameters  for  donors  and  acceptors  incorporated  into  the 
model,  in  cases  of  n-  or  p-type  samples,  respectively. 

All  the  above  discussed  Auger  processes  (Figs.  4(a)  and 
4(b)  and  Figs.  S(a)-S(d)]  are  similar  to  the  two  center 
problems  analyzed  in  some  papers.32’50,41-41, 51,43,46  Fig¬ 
ure  4(a)  shows  an  example  of  nonradiative  recombination 
involving  the  interaction  of  a  core  excited  isoelectronic 
trap  (Nq)  with  an  electron  trapped  on  a  separate  center 
(N-),  or  an  electron  trapped  on  a  core  excited  REI  trap 
(NT.).  The  recombination  rates  for  these  processes  are 
proportional  to  NB,  and  the  concentrations  of  or  NT. 
traps  such  that 

rBTl=:^BT1^0^-  >  rBT2=-®BT2^0^-  > 

where  B^, .  and  are  the  Auger  coefficients. 

According  to  the  theory  developed  by  Langer, 45,48  the 
two  additional  paths  of  energy  transfer  from  the  excited 
REI  trap  are  the  first  to  free  electrons  (holes)  with  Auger 
coefficients  BBT}  (B  BT4),  and  second  to  electrons  bound  at 
the  shallow  donors  (acceptors)  [Figs.  5(a)— 5(d)],  and  are 
also  attributed  to  the  luminescence  quenching.  The  two 
processes  shown  in  Figs.  5(c),  and  5(d)  involve  interaction 
of  a  core-excited  REI  trap  with  an  electron  (hole)  on  a 
neutral  donor  (acceptor),  and  are  characterized  by  Auger 
coefficients  BBT5  and  BBT6,  respectively.  The  first  process 
is  important  during  the  exciting  pulse,  and  the  second 
processes  are  expected  to  be  important  at  low  tempera¬ 
tures  where  all  electrons  (holes)  are  frozen  out  onto  neu¬ 
tral  donors  (acceptors). 

We  can  now  complete  the  formal  description  of  the 
models  by  deriving  the  differential  equations  for  the 
REI-trap  energy-transfer  processes  and  recombination 
kinetics.  The  differential  equations  govern  the  variations 
with  the  time  of  the  concentrations  of  various  com¬ 
ponents  under  band-to-band  excitation.  By  consulting 
Figs.  1-5  (there  are  six  possible  states  in  REI  traps,  and 
two  in  shunt  paths  and  generated  free  carriers),  they  may 
be  readily  written  down: 
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Ar0=JV-AT_  -Hx~No  ~N- ~Nx  * 
Nt=Nf+Ns~  , 


p=n +JV_ +N- + 


(8) 

(9a) 

(9b) 

(10) 


Equations  (l)-(5)  govern  the  negatively  charged-rare 
earth  isoelectronic  trap  populations  (1),  the  neutral  exci- 
ton  occupied  rare-earth  isoelectronic  trap  populations  (2), 
the  neutral  core  excited  rare-earth  isoelectronic  trap  pop¬ 
ulations  (3),  the  negatively  charged  core  exited  rare-earth 
isoelectronic  trap  populations  (4),  and  the  exciton  occu¬ 
pied  core  excited  rare-earth  isoelectronic  trap  popula¬ 
tions  (5),  respectively.  Equations  (6)  and  (7)  represent  the 
changes  in  the  total  free  electron  (n)  and  hole  (p)  popula¬ 
tions,  respectively.  Equation  (8)  governs  the  density  of 
shunt  path  traps.  Equations  (9a)  and  (9b)  state  the  con¬ 
stancy  of  the  total  concentrations  of  REI  traps  and  shunt 
paths,  respectively,  and  Eq.  (10)  is  the  neutrality  condi¬ 
tion  equation.  G—gfit—T)  is  the  generation  rate  of 
free  electrons  and  holes  by  above  band-gap  laser  excita¬ 
tion.  The  electron-hole  pair  generation  rate  is  controlled 
by  a  unit  step  function  (T  is  the  pulsewidth) 

which  takes  on  the  value  of  zero  or  unity  according  to 
whether  its  argument  is  less  than  or  equal  to  zero  or 
greater  than  zero.  The  volume  generation  rate  is  approx¬ 
imated  by  g(z)stal  exp[  —az  ],  where  a  is  the  absorption 
coefficient,  and  /  is  photon  flux  in  #  photons/(cm2/sec). 

For  n-type  material  the  equations  are  even  more  com¬ 
plicated  (see  the  Appendix);  similar  equations  hold  for  p- 
type  material  Such  systems  of  coupled,  first-order  stiff 
nonlinear  differential  equations  require  specialized  nu¬ 
merical  integration  routines  designed  specifically  for  stiff 
systems.  The  kinetic  equations  were  solved  numerically 
as  a  function  of  excitation  intensities,  using  as  fitting  pa¬ 
rameters  the  time  constants  and  rate  coefficients  defined 
in  Table  II.  To  describe  the  buildup  and  decay  kinetics  of 
rare-earth  luminescence  in  a  semiconductor  host,  we  as¬ 
sumed  only  band-band  pulse  excitation.  The  above  set  of 
ten  equations  (1)— (10)  was  reduced  to  eight  dimensionless 
differential  equations  which  is  more  convenient  for  nu¬ 
merical  solutions.  To  solve  this  system  we  assumed,  for 
simplicity,  that  the  above  band-gap  excitations  take  place 
at  low  enough  temperatures  (in  our  experiment  8.6  K)  so 
that  thermal  activation  of  the  trapped  carriers  is  negligi¬ 
ble.  That  is,  terms  explicitly  dependent  on  temperature 
were  ignored,  and  only  trapping,  transferring,  and  recom¬ 
bination  transitions  were  considered.  The  radiative 
recombination  of  excitons  bound  to  “simple”  isoelectron¬ 
ic  impurities  has  a  long  lifetime  ranging  from  a  few  hun¬ 
dred  to  a  few  thousand  nanoseconds.24,28-32  That  time  is 
much  longer  than  the  energy-transfer  time  r^-  from  an 
exciton  bound  to  a  REI  trap  to  core  states.  The  measure¬ 
ments  on  n-type  InP:Yb  show  that  this  transfer  time  is 
much  less  than  10  ns.7  The  above  facts  explain  why  we 
do  not  observe  the  luminescence  of  excitons  bound  to 
REI  traps.  The  numerical  solution  of  Eqs.  (1)— (10)  was 
obtained  using  the  parameters  shown  in  Table  II.  Gor¬ 


TABLE  II.  Parameters  describing  rise  and  decay  kinetics  of 
InP:Yb. 


Symbol 

Unit 

Parameter  value 

References 

Ti 

s 

2X10'7 

24,28-32 

r3 

s 

11.6X10'4 

9,44 

rBT 

s 

5X10'7 

31,32,41-43 

Ttxr 

s 

5X 10-7 

31,32,41-43 

rM 

s 

5X10“' 

33(b),33(c) 

Tm 

s 

3X10“* 

33(b),33(c) 

rXT 

s 

1.25  X  10”* 

7 

r„ 

s 

5X10'" 

33(b),33(c) 

V 

s 

5X10“" 

33(b),33(c) 

Bt 

cmVs 

4X  10~'° 

Bfu 

ctnVs 

1X10'" 

31,32,41-43 

Bm\ 

cmVs 

1.2XHT'3 

31,32,41-43 

Barn 

cmVs 

1.2X10-13 

31,32,41-43 

Ban 

cmVs 

4X10-'3 

45-48 

■®BT« 

cmJ/s 

4X  10~u 

45-48 

N, 

cm”3 

8X10'4 

N 

cm-3 

5X10'7 

G 

#  photons/(cnr 

's)  2X10JO-1.5X1023 

don  and  Allen4**1  and  Ayling  and  Allen4*61  and  Langer, 
Van  Hong4S,4*u),48<bl  made  quantitative  measurements  of 
the  luminescence  efficiency  and  lifetime  in  ZnS:Mn, 
ZnSe:Mn,  CdF2:Mn,  and  CdF2.-Gd,  where  they  obtained 
a  value  for  the  Auger  quenching  coefficients  by  free  elec¬ 
trons  and  an  electron  bound  on  shallow  donors.48  A 
surprising  result  is  that  the  value  of  Auger  coefficients  for 
centers  in  CdF2  («5X  10-15  cm3s-1)  is  about  four  order 
of  magnitude  smaller  than  the  value  for  ZnS  and  ZnSe 
(5X10  10  cm3s  ’).  Klein,  Femeaux,  and  Henry4*11  es¬ 
timated  the  Auger  coefficient  for  the  InP:Fe  system  to  be 
=6. 7X 10-10  cm3s-1.  Excellent  reviews  of  theory  and 
experimental  data  of  Auger  processes  and  values  of 
coefficients  are  given  by  Landsberg  and  Adams.42,43  The 
values  of  the  parameters  chosen  for  the  calculations  are 
estimated  from  experimental  data  obtained  from  similar 
“simple”  isoelectronic  traps  (see  references  in  Table  II). 
We  believe  that  all  constants  are  realistic,  and  charac¬ 
teristic  for  InP:Yb.  The  volume  generation  range  G  of 
e-h  pairs  wa  chosen  in  the  same  range  as  in  the  experi¬ 
ment.44  The  kinetics  of  luminescence  of  InP:Yb  as  the 
functions  of  excitation  intensity  was  simulated  by  repeat¬ 
ing  the  numerical  calculations  for  several  different  values 
of  the  generation  rate.  We  simulated  the  kinetics  of  the 
photoluminescence  measurement  by  choosing  square  gen¬ 
erating  pulses  with  60-ps  duration,  adequate  to  establish 
a  quasiequilibrium  luminescence  intensity  during  excita¬ 
tion.  The  luminescence  intensities  are  proportional  to 
N_,  and  N£.  The  last  two  terms  can  be  ignored  be¬ 
cause  they  introduce  very  small  contributions  to  the  total 
luminescence  intensity,  mainly  during  the  excitation 
pulse.  The  numerically  simulated  luminescence  rise  and 
decay  are  shown  in  Fig.  6  as  a  function  of  generation 
rates  (G),  and  sets  of  the  parameters  from  Table  II.  The 
profiles  shown  the  buildup  of  luminescence,  allowing  the 
process  to  reach  a  steady-state  value  for  a  given  excita¬ 
tion  intensity,  and  the  decay  after  switching  off  the  exci- 
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tation  pulse.  All  profiles  shown  on  Fig.  6  are  normalized 
to  unity  at  maximum.  The  rise  times  given  in  the  figure 
decrease  with  the  increasing  generation  rates.  At  the  low 
generation  rate  of  2X102®  photon/(cm3s)  the  buildup 
curve  (which  can  be  well  approximated  by  a  single  ex¬ 
ponential  curve)  has  a  rise  time  of  1 1.98  fis,  while  at  the 
generation  rate  of  1.5  X  1023  photons/(cm3s)  the  rise  time 
is  1.7  fis.  The  decay  profiles  for  the  above  two  generation 
rates  shown  in  Fig  6  are  1 1.7  fi  and  6.9  /is,  respectively. 
The  remaining  four  decay  times  fall  in  that  window.  The 
rise  and  decay  times  for  simulated  kinetics  are  well  ap¬ 
proximated  by  single  exponential  processes,  while  the  ex¬ 
perimental  data  are  described  better  by  double  exponen¬ 
tial  functions.  Figure  7  shows  the  computed  rise  and  de¬ 
cay  times  as  a  function  of  excitation  intensities  (genera¬ 
tion  rate).  The  experimental  data  are  imposed  on  the 
computed  curve.  The  circle  shown  in  Fig.  7  represents 
experimental  rise  times  obtained  from  fitting  to  a  single 
exponential  function  (with  the  coefficient  of  determina¬ 
tion  r2«0.98).  The  squares  stand  for  the  rise  times  of  the 
dominant  component  of  the  experimental  data  fitted  to  a 
double  exponential  function.  In  Fig.  7  the  upper  solid 
line  is  the  computed  decay  time  and  the  dots  represent 
the  experimental  decay  times  of  the  dominant  component 
of  the  double  exponential  fitting.  The  second  small  ex¬ 
ponential  component  has  the  decay  times  in  the  range 
0.9- 1.8  /is.  The  double  exponential  rise  and  decay  times 
observed  in  the  experiment  are  probably  related  to  other 
energy-transfer  processes  (not  incorporated  to  our  model) 
from  Yb  to  other  accidental  impurities  where  the  energy 
is  dissipated  nonradiatively.  Iron  is  a  common  contam¬ 
inant  in  metal-organic  chemical-vapor  deposition 
(MOCVD)  grown  Yb-dopcd  InP  crystals.  The  energy- 
transfer  process  from  Yb3+  to  Fe  ions  is  attributed  to 
thermal  quenching  of  the  luminescence  intensity  and  de¬ 
cay  time  with  an  increase  in  temperature.44  The  detailed 
analysis  of  the  parameters  and  their  influence  on  the  rise 


FIG.  6.  Numerical  solutions  of  Eqs.  (1)-(10).  The  rare-earth 
luminescence  intensity  vs  time  for  different  generation  rates  and 
the  set  of  parameters  shown  in  Table  II.  The  set  of  curves  a  -f 
represents  computed  profiles  (normalized  to  unity  at  maximum), 
showing  the  buildup  of  luminescence,  the  steady-state  value, 
and  the  decay  after  switching  off  tat  60  fis)  the  excitation  pulse. 


FIG.  7.  The  numerically  computed  rise  time  (solid  line)  of 
rare-earth  luminescence  (InP:Yb)  as  a  function  of  generation 
rates.  The  circles  represent  the  experimental  values  of  the  rise¬ 
time  constants  of  ytterbium  luminescence  in  InP  at  8.4  K  (fitted 
to  single  exponential  function,  with  the  coefficient  of  determina¬ 
tion  rJ« 0.986).  The  squares  stand  for  the  rise-time  constants  of 
the  dominant  component  of  the  experimental  data  fitted  to  the 
double  exponential  function.  The  dots  stand  for  the  decay-time 
constants  of  the  dominant  component  of  the  experimental  data 
fitted  to  double  exponential  function. 


time,  and  the  efficiency  of  the  steady-state  luminescence 
and  decay  time  will  be  published  elsewhere.  The  numeri¬ 
cally  simulated  luminescence  rise  and  decay  measure¬ 
ments  show  a  good  overall  quantitative  agreement  with 
the  experiment  over  a  wide  range  of  generation  rates 
(Figs.  6  and  7).  Finally,  the  proposed  model  can  be 
refined  by  taking  into  account  the  surface  recombination 
and  carrier  diffusion  processes,  which  may  play  impor¬ 
tant  roles  in  the  reduction  of  the  overall  photolumines¬ 
cence  emission.53 

IV.  POSSIBLE  QUENCHING  MECHANISMS 
AND  TEMPERATURE  DEPENDENCE 

Discussed  below  are  several  possible  nonradiative  pro¬ 
cesses  that  quantitatively  explain  the  kinetics  (rise  and 
decay  of  luminescence  of  the  REI  trap)  and  temperature 
dependence  of  the  observed  quenching.  The  Auger 
quenching  processes  defined  in  Sec.  Ill  are  shown  in  Figs. 
3-5.  Processes  [Figs.  4(a)  and  4(b)  and  5(a)  and  5(b)]  non¬ 
radiative  recombination  involving  the  interaction  of  a 
core  excited  REI  trap,  with  an  electron  trapped  on  a 
separate  REI  trap.  The  next  two-center  recombination 
involves  the  interaction  of  the  excited  REI  trap  with  an 
electron  (hole)  on  a  separate  neutral  donor  (acceptor),  re¬ 
spectively.  The  rate  equations  describe  the  nonradiative 
recombination  involving  the  interaction  of  a  core  exited 
REI  trap,  with  a  free  electron  or  hole,  respectively.  The 
last  process  is  a  very  efficient  nonradiative  deexcitation 
mechanism  of  localized  centers  (Mn2+)  in  ZnS,  ZnSe,  and 
CdF245,47  and  Gd3+  and  Tb3+  in  CdF2.48  The  free- 
carrier  Auger  quenching  processes  are  important  at 
higher  temperatures,  and  were  shown  to  be  much  more 
efficient  than  the  Auger  process  due  to  the  shallow 
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donors.  The  last  one  is  important  only  at  low  tempera¬ 
tures.47  The  most  probable  nonradiative  mechanism 
quenching  the  rare-earth  luminescence  in  III-V  semicon¬ 
ductors  is  the  Auger  energy  back  transfer  mechanism 
shown  schematically  in  Figs.  3(a)  and  3(b).  The  rare- 
earth-excited  isoelectronic  trap  occupied  by  electron  N- 
or  exciton  N \  can  transfer  energy  to  the  trapped  electron 
or  exciton  rather  than  to  the  radiative  field;  the  electron 
is  consequently  ejected  deep  into  the  conduction  band.  It 
is  emphasized  that  at  low  temperatures  the  only  electrons 
in  the  conduction  band  come  from  either  the  exciting 
source  or  from  Auger  ionization  of  the  REI  traps  or  neu¬ 
tral  donors  (in  the  case  of  n-type  samples).  The  nonradia¬ 
tive  recombination  involving  the  interaction  of  the  core- 
exited  REI  trap  with  a  free  electron  (or  hole)  is  important 
in  seminsulating  InP  crystals  only  during  the  excitation 
pulse.  In  the  dark  the  electron  concentratio  a  at  a  tem¬ 
perature  of  300  K  is  in  the  range  1 X  10*  cm-3.  However, 
the  experimental  decay  time  of  Yb  luminescence  is  ~  12 
fis  in  SI-, n-,  and  •?  'ype  InP:Yb.3,4,6 

The  other  ir  »  rnisms  of  the  nonradiative  recombina¬ 
tion  of  the  ex.  icc  states  of  localized  centers  in  semicon¬ 
ductors  are  first,  the  mulliphonon  relaxation  processes,33 
and  second,  a  migration  of  energy  and  cross  relaxation 
processes.36,37  The  probability  of  the  multiphonon  relaxa¬ 
tion  processes  is  dependent  upon  the  type  of  coupling 
with  the  lattice  vibrations  and  the  phonon  frequency  dis¬ 
tribution.  The  ion-host  lattice  interaction  for  RE3+  4/" 
electrons  is  characteristic  of  weak  coupling,  and  the  mul¬ 
tiphonon  emission  rates  exhibit  an  approximately  ex¬ 
ponential  dependence  on  the  energy  gap  to  the  next 
lowest  level.  For  a  single  frequency  p  phonon  process  the 
nonradiative  rate  for  narrow-line  emission  is  described  by 
Kiel's  formula38 

Wr,(T)=W>f{Q)[\+nY  , 

where  W^y(0)  is  the  rate  at  low  temperature,  p  is  the 
number  of  phonons  involved,  n  is  the  occupancy  of  the 
phonon  modes: 

n  =[exp(few/k7')—  l]-1  , 

and  the  average  phonon  energy  Ho)  =  A E/p.  The  single 
frequency  model  seems  to  be  an  oversimplification,  and 
taking  a  weighted  average  over  the  phonon’s  spectrum, 
or  consideration  of  the  continuous  nature  of  the  phonon 
spectrum  would  be  more  appropriate.39 

The  results  of  many  studies  demonstrate  that  for  given 
host  crystals,  the  most  critical  parameter  affecting  the 
rate  of  multiphonon  emission  is  the  energy  gap  to  the 
nearest-  lower  level.  If  the  energy  gap  to  the  next  lower 
state  is  sufficiently  large  the  nonradiative  multiphonon 
rate  is  negligible  compared  to  the  radiative  rate,  and  this 
is  the  situation  for  most  rare-earth  ions  in  III-V  semicon¬ 
ductors.  For  example,  in  InP:  Yb3+  the  2Fi/2  ~2Fm  en_ 
ergy  gap  is  9895  cm-1,  in  GaAs:  Er3+  the  47n/2  —  4/ 15/2 
energy  gap  is  6051  cm-1,  and  in  GaAs:Nd3+  the 
*Fin—*I nn  energy  gap  is  10636  cm-1.  The  highest  en¬ 
ergy  phonons  in  InP  and  GaAs  are  345.2  and  292.7  cm  1 
LO  phonons,  respectively.  Nonradiative  decay  for  the 
three  cases  above  would  require  the  generation  of  27,  21, 


and  36  LO  phonons,  respectively.  Studies  in  different 
host  materials,351^,■35,b,  show  that  nonradiative  decay  in¬ 
volving  the  generation  of  more  than  five  phonons  is  weak¬ 
er  than  the  radiative  process. 

When  the  concentrations  of  the  rare-earth  luminescent 
ions  are  higher,  or  centers  create  pairs,  allowance  must  be 
made  of  the  possibility  of  interaction  between  centers. 
This  interaction  may  be  too  weak  to  modify  the  energy 
levels,  but  yet  be  adequately  strong  to  enable  transfer  of 
energy  from  one  to  another. 380,1,37 

The  migration  of  energy  in  InP  doped  with  Yb3+ 
occurs  through  resonant  interaction.  The  energy  level  di¬ 
agram  of  Yb3+  in  a  crystal  field  of  the  tetrahedral  sym¬ 
metry  (Td)  in  InP  indicates  that  the  only  possible  interac¬ 
tion  between  the  Yb3+  ions  at  low  temperature  is  reso¬ 
nant  interaction  between  the  lower  Stark  sublevel  of  the 
ground  state2F7/2  and  the  lower  Stark  to  the  accidental 
impurity  with  an  absorption  spectrum  overlapping  with 
that  of  the  YbJ+  ion,  emission  is  the  most  probable 
quenching  mechanism.  In  Nd3+ -doped  GaAs  the  energy 
of  a  single  excited  ion  may  be  transferred  to  the  same  lev¬ 
el  of  an  identical  ion  (resonant  energy  migration),  rnd 
finally  to  the  quenching  center.32,54,60,61,44  The  second 
possible  degradation  of  excitation  energy  is  the  cross  re¬ 
laxation  process.  The  presence  of  appropriate  intermedi¬ 
ate  levels  in  the  Nd3+  ion  (4/M/ 2  and  4/, 3/2)  between  the 
ground  4/9/2  and  excited  4/3/2  states  makes  effective  the 
cross  relaxation  process  and  increases  the  probability  of 
nonradiative  quenching  processes.38,62,63  All  the  above- 
discussed  transfer  processes  are  temperature  dependent 
and  can  be  useful  in  explaining  the  quenching  of  lumines¬ 
cence  of  rare  earths  in  semiconductors.  A  full  account  of 
our  results  on  the  temperature  dependence  of  rise,  decay, 
and  quenching  mechanisms  will  be  published  at  a  latter 
data. 

V.  SUMMARY 

In  this  work  we  have  developed  a  kinetics  model  of  en¬ 
ergy  transfer  from  the  host  lattice  to  the  localized  core 
excited  states  of  rare-earth  isoelectronic  structured  traps. 
According  to  Robbins  and  Dean17  the  formation  of 
bound  states  at  the  structured  cationic  isoelectronic  im¬ 
purities  is  fairly  common,  which  is  in  contrast  to  the  situ¬ 
ation  normally  found  for  anionic  substituents  in  semicon¬ 
ductors.  The  outer  electron  configurations  of  RE3+  ions 
are  the  same  (5s25p6).  Among  rare-earth  ions  are  those 
that,  upon  replacing  the  element  from  column  III  in  III- 
V  compounds  that  are  isovalent  concerning  outer  elec¬ 
trons  of  RE3+  ions,  create  isoelectronic  traps  in  III-V 
semiconductors.  We  have  evidence  that  the  other  RE 
ions  in  III-V  semiconductors  can  occupy  different  sites 
(not  only  substitutional).  The  rare-earth  isoelectronic 
trap  must  not  necessarily  be  the  “pure”  substitutional 
center;  if  the  rare-earth  ions  are  very  active  chemically, 
they  can  create  a  more  complex  center  involving  other 
impurity  or  native  defects.  This  conclusion  is  supported 
by  the  fact  that  the  atomic  covalent  radii  (ionic  RE3 +  )  for 
all  rare  earths  are  bigger  than  atomic  radii  of  Ga  in  In 
that  they  replace.  Pauling’s  electronegativity  of  rare- 
earth  elements  are  in  the  range  of  1.1-1.25.  and  are 
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smaller  than  Ga  (1.81)  and  In  (1.78)  for  which  they  sub¬ 
stitute.  The  rare-earth  isovalent  traps  that  we  can  call 
isoelectronic  “structured”  impurities17  posst-s  the 
unfilled  4 /"  core  shell.  The  structured  luminescence 
arises  from  intraconfiguration  /-/  transitions  in  the  core 
of  the  isoelectronic  “structured”  impurities. 

The  presence  of  low-lying  empty  core  orbitals  in  rare- 
earth  impurities  introduces  new  excitation  and  recom¬ 
bination  phenomena.  To  adequately  describe  the  energy 
transfer  to  the  REI  trap,  and  the  buildup  and  decay 
kinetics  of  rare-earth  luminescence,  we  have  considered 
six  separate  states  of  the  REI  impurity  (unoccupied,  elec¬ 
tron  occupied,  electron  occupied  excited,  neutral  exciton 
occupied,  exited  electron  occupied,  and  excited  exciton 
occupied).  The  energy-transfer  processes  occur  through 
an  Auger  mechanism,  where  the  recombination  energy  of 
the  bound  electron  with  the  free  hole  is  transferred  non- 
radiatively  to  the  core  states.  The  second  energy-transfer 
process  is  the  transfer  of  energy  from  the  bound  exciton 
on  the  REI  trap  to  the  core  states.  If  the  initial  and  final 
states  are  not  resonant  (in  both  mechanisms),  the  energy 
mismatch  must  be  accommodated  by  emission  or  absorp¬ 
tion  of  phonons. I7,38(bl- 36  Furthermore,  we  discussed  the 
details  of  several  quenching  processes  which  are  incor¬ 
porated  into  the  kinetics  equations.  We  derive  two  sets  of 
differential  equations  for  SI-,  and  n-  ip-)  type  semiconduc¬ 
tors  governing  the  kinetics  of  rare-earth  luminescence. 
The  nonradiative  pathways  present  alternative  recom¬ 
bination  possibilities  for  the  electrons  and  holes  and  are 
incorporated  in  both  models.  Equations  governing  the 
transfer  of  energy  processes  and  recombination  kinetics 
have  been  determined  and  solved  by  a  numerical  method 

_ I 


to  derive  the  time  dependence  of  the  rise  and  decay  kinet¬ 
ics  as  a  function  of  excitation  intensity.  The  parameters 
used  are  realistic,  but  only  approximately  known  (see 
Table  II).  The  numerically  simulated  luminescence  rise 
and  decay  times  show  a  good  overall  quantitative  agree¬ 
ment  with  experiment  over  a  wide  range  of  generation 
rate  values  (see  Figs.  6  and  7).  Finally,  the  proposed 
model  can  be  refined  by  taking  into  account  the  quench¬ 
ing  processes  discussed  above,  and  the  surface  recom¬ 
bination  and  carrier  diffusion  processes,  which  plays  an 
important  role  in  reduction  of  the  overall  photolumines- 
cence  emission.33 
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APPENDIX 

Consulting  Figs.  1-5,  the  equations  describing  the 
transfer  and  recombination  kinetics  in  n-type  semicon¬ 
ductor  may  be  written  as  follows: 
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Finally,  the  equation  indicating  the  neutrality  condition  is 


p=n+N~+Nt + 


(All) 


Equations  (A1)-(A5)  govern  the  negatively  charged  REI-trap  population  (Al),  the  neutral  exciton  occupied  REI-trap 
population  (A2),  the  neutral  core  excited  REI-trap  population  (A3),  the  negatively  charged  core  excited  REI-trap  popu¬ 
lation  (A4),  and  the  exciton  occupied  core  excited  (AS)  rare-earth  isoelectronic  trap  populations,  respectively.  Equa¬ 
tions  (A6)  and  (A7)  represent  the  changes  in  the  total  free-electron  and  hole  populations,  respectively.  Equation  (AS), 
governs  the  density  of  neutral  shunt  path  traps.  Equation  (A9)  governs  the  density  of  neutral  donors.  Equations  (AlO) 
state  the  constancy  of  the  total  concentrations  of  RBI  traps,  shunt  paths,  and  donors,  respectively,  and  Eq.  (Al  1)  is  the 
neutrality  condition  equation.  The  symbols  used  in  this  set  of  equations  are  the  same  as  for  SI  semiconductors. 
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KINETICS  OF  THE  LUMINESCENCE  OF  ISOELECTRONIC 
RARE-EARTH  IONS  IN  III-V  SEMICONDUCTORS 
H  J.  Lozykowski  ,  Department  of  Electrical  and  Computer  Engineering,  and  Condensed 
Matter  &  Surface  Sciences  Program  Ohio  University  Athens,  Ohio  45701 
ABSTRACT 

In  this  work  we  have  developed  a  model  for  the  ldnedcs  of  the  energy  transfer  from 
the  host  lattice  to  the  localized  core  excited  states  of  rare  earth  isoelectronic  structured  traps 
(REI-trap).  We  have  derive  a  set  of  differential  equations  for  semi-insulating  semiconductor 
governing  the  kinetics  of  rare  earth  luminescence.  The  numerically  simulated  rise  and  decay 
times  of  luminescence  show  a  good  quantitative  agreement  with  the  experimental  data 
obtained  for  InP:Yb,  over  a  wide  range  of  generation  rates. 

INTRODUCTION 

The  investigation  of  the  luminescence  properties  of  rare  earth  doped  III-V  is  of  great 
interest  both  from  scientific  and  application  points  of  view.  The  scientific  interest  is  related 
to  the  uniqueness  of  optical  and  electrical  properties  of  rare  earth  impurities  in  semiconductor 
hosts.  Among  the  rare  earth  doped  13I-V  semiconductors  InP:Yb  has  been  the  most 
extensively  studied  [1-6].  In  this  paper  we  discuss  only  the  structured  isoelectronic  traps 
(REI-trap)  in  III-V  semiconductors  introduced  by  RE+J  ions  replacing  the  element  from 
column  HI.  Furthermore,  we  develop  a  model  of  die  luminescence  kinetic  that  describes  the 
energy  transfer  from  the  host  to  the  REI-trap  core  states,  and  the  recombination  and 
quenching  processes.  Study  of  die  rise  and  decay  times  at  different  excitation  intensities, 
temperatures,  can  provide  important  information  about  the  energy  transfer  and  recombination 
(radiative,  and  non-radiative)  processes. 

THEORETICAL  FORMULATION 

An  isoelectronic  center  can  form  bound  states  because  of  a  short  range  central-cell 
potential.  According  to  Thomas  [7],  the  primary  factors  affecting  the  binding  potential  are 
the  electronegativity  and  the  size  differences  between  the  impurity  and  the  host  ion  which  it 
replaces.  It  is  found  experimentally  that  only  very  large  atoms  or  very  small  atoms  produce 
isoelectronic  traps  because  they  create  large  lattice  distortion  induced  by  the  substitution.  The 
above  conclusion  is  supported  by  the  fact  that  the  atomic  covalent  radii  (ionic  RE+3)  for  all 
rare  earths  are  bigger  than  atomic  radii  of  Ga  and  In  that  they  are  replacing.  Pauling’s 
electronegativity  of  rare  earth  elements  is  in  the  range  of  1.1-1.25,  and  is  smaller  than  Ga 
(1.81)  and  In  (1.78)  for  which  it  substitutes.  If  the  rare  earth  ions  replace  the  element  from 
column  III  in  III-V  compounds  (that  are  iso  valent  concerning  outer  electrons  of  RE+3  ions) 
,  they  create  isoelectronic  traps.  The  rare  earth  isovalent  traps  that  we  can  call  isoelectronic 
"structured"  impurities  [8]  possess  unfilled  4f  core  shells.  The  luminescence  structure 
arises  from  intra-configurational  f-f  transitions  in  the  core  of  the  isoelectronic  "structured" 
impurities. 

The  striking  feature  of  excitons  bound  to  the  isoelectronic  traps  is  a  long 
luminescence  decay  time,  ranging  from  a  few  hundred  to  few  thousand  nanoseconds  [9,10- 
12].  The  lifetimes  of  neutral  donor  or  acceptor  bound  excitons  in  direct  gap  semiconductors 
are  in  the  range  of  nanoseconds.  The  "simple"  isoelectronic  center  in  III-V  materials  can 
exist  in  three  possible  states  instead  of  two  as  in  the  case  of  the  Shockley-Read-Hall  (SRH) 
recombination  model  [11],  In  the  case  of  rare  earth  isoelectronic  traps,  the  kinetics  model 
is  even  more  complicated  because  of  an  energy  transfer  processes  between  the  localized  state 
in  the  forbidden  gap  of  the  host,  and  the  localized  core  state  of  structured  isoelectronic 
impurities.  There  are  three  possible  mechanisms  of  energy  transfer.  The  first  is  the  energy 
transfer  from  excitons  bound  to  "structured"  isoelectronic  centers  to  the  core  electrons.  It 
takes  place  as  a  result  uf  the  electrostatic  perturbation  between  the  core  electrons  and  the 
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Fig.  I(a-j).  Schematic  diagram  of  trapping, 
recombination,  (thermalization),  exciton 
formation  and  dissociation  processes 
involving  structured  REI-  traps  (represented 
by  ellipse)  with  the  atomic-like  4P  core 
states.  A  schematic  representation  of  the 
energy  transfer  processes  to  core  states  and 
radiative  and  nonradiative  transitions. 
Energy  back  transfers  processes,  and  Auger 
quenching  processes  involving  interaction  of 
core  excited  REI-trap  with  free  electrons 
(holes). 
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exciton,  (effective-mass-like  particles)  [8].  The  second  mechanism  is  the  transfer  of  energy 
to  the  core  electrons  involving  the  "structured*  isoelectronic  trap  occupied  by  electron  (hole) 
and  free  hole  (electron)  in  the  valence  (conduction)  band.  The  third  mechanism  is  the  transfer 
through  an  inelastic  scattering  process  in  which  the  energy  of  a  free  exciton  near  a 
"structured"  trap  is  given  to  the  localized  core  excited  states  [8).  If  the  initial  and  final  states 
are  not  resonant,  the  energy  mismatch  must  be  distributed  in  some  way,  e.g.  by  phonon 
emission  or  absorption  [8, 13].  If  the  atomic  core  excitations  are  strongly  coupled  to  the  host 
phonons,  the  energy  transfer  probability  is  likely  to  be  higher  [8].  Strong  phonon  coupling 
may  also  be  desirable  in  ensuring  that  relaxation  down  the  ladder  of  core  excited  states  occurs 
quickly,  thus  preventing  back  transfer.  However,  for  efficient  radiative  recombination,  the 
phonon  coupling  should  not  be  strong,  in  order  to  prevent  core  de-excitation  by  nonradiative 
multiphonon  emission.  In  this  regard  the  rare  earth  "structured"  impurity  seems  to  be  ideal. 
KINETIC  MODEL  AND  COMPUTATION  RESULTS 

The  luminescence  kinetic  model  involving  RE+S  "structured"  isoelectronic  centers, 
is  complicated  because  the  center  can  exist  in  six  possible  states.  Furthermore  the  energy 
transfer  processes  between  the  localized  state  in  the  forbidden  gap  and  core  states  complicate 
the  model.  We  assume  that  the  isoelectronic  trap  is  an  electron  trap  such  as  Yb+J  in  semi- 
insulating  (SI)  InP.  This  trap  may  exist  in  six  distinct  states:  (1)  neutral  unoccupied  trap 
(concentration  NJ,  (2)  negatively  charged  (concentration  N.),  (3)  exciton  occupied 
^(concentration  NJ,  (4)  neutral  excited  (concentration  hQ,  (S)  excited  negatively  charged 
(concentration  N*),  and  (6)  excited  exciton  occupied  (concentration  hQ.  When  isoelectronic 
traps  are  present  in  n  or  p-type  materials  the  model  will  be  different.  The  asymmetry  between 
n  and  p  type  semiconductors  results  from  the  fact  that  the  isoelectronic  center  binds  only  an 
electron  (hole).  Thus  in  a  n-type  material,  a  fraction  of  the  isoelectronic  traps  will  be 
occupied  by  electrons  even  before  the  sample  is  excited,  while  in  p  -type  material,  all  the 
isoelectronic  centers  will  be  empty  .  Figure  1  shows  the  physical  models,  symbols,  and 
definitions  used  in  the  analysis  of  excitation,  recombination  and  quenching  of  “structured" 
isoelectronic  centers.  By  fitting  the  calculation  to  the  experimental  data,  we  can  estimate 
important  parameters  related  to  energy  transfer  from  the  lattice  to  RE+3  centers,  the  Auger 
processes,  and  the  temperature  quenching  mechanisms.  The  N.  center  can  be  transformed 
to  the  neutral  rare  earth  core  excited  center  N*,  through  an  Auger  process  where  the 
recombination  energy  of  the  bound  electron  with  free  hole  is  transferred  nonradiatively  to  the 
core  states.  The  N.  may  also  recombine  radiatively  or  nonradiatively  with  a  hole  in  the 
valence  band.  The  last  three  processes  transform  N.  center  into  Nc  neutral  trap.  The  N*  may 
capture  an  electron  and  be  transformed  into  N*  center.  The  Nt  may  lose  the  electron  through 
processes  described  above,  or  by  auto  de-ionization  Auger  process.  This  Auger  nonradiative 
luminescence  quenching  mechanism  of  the  excited  REI-trap  with  characteristic  time  rBT  is 
shown  in  Fig. If.  The  exciton  bound  to  a  REI-trap  with  excited  core  states  (NQ,  can  back 
transfer  energy  from  core  excitations  to  the  exciton,  and  dissociate  it,  with  the  liberation  of 
an  electron  or  hole.  Once  formed,  the  N.,  and  N!  states  can  also  be  converted  into  Nc,  and 
N*  states,  respectively,  by  thermal  ionization  or  to  N„  and  N*  ,  respectively,  by  capturing 
a  hole  from  the  valence  band  with  time  rp,.  The  second  most  important  energy  transfer 
process  is  shown  in  Fig.  le.  In  this  process  energy  is  transferred  from  the  bound  exciton  (on 
isoelectronic  trap  NJ  to  the  core  states  with  the  rate  r,T.  The  exciton  bound  to  isoelectronic 
traps  N,  and  N'  can  thermally  dissociate  by  several  processes.  It  can  dissociate  into  a  free 
exciton  (X)  and  neutral  Nc,  or  neutral  excited  (N‘)  REI-trap  respectively.  The  REI-traps  N', 
N!  and  N*  may  also  be  de-excited  to  NOJ  N.  and  N,  respectively  by  radiative  transitions 
(time  t})  between  the  4fl,  core  states  .The  nonradiative  decay  channel  of  free  carriers  through 
an  additional  type  of  trap  N,  the  "shunt  path"  is  shown  in  Fig.lc.  Figure  l.h,i,j  show 
examples  of  nonradiative  recombination  involving  the  interaction  of  a  core  excited  REI-trap 
(N*)  with  an  electron  trapped  on  a  separate  center  (N.),  an  electron  trapped  on  a  core  excited 
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Table  I  Parameters  Describing  Rise  and  Decay  Kinetics 
of  InP  :  Yb. 


Symbol 

Unit 

Parameter  Value 

r > 

sec 

2  x  JO’ 

sec 

11.6  x  104 

sec 

1  x  10* 

r«nc 

sec 

I  x  10* 

Tm 

sec 

5  x  I04 

Tm 

sec 

3  x  to4 

sec 

1.25  X  104 

T. 

sec 

5  x  104* 

V 

sec 

5  X  10*' 

B, 

cmVsec 

4  x  I044 

Bm 

cmVsec 

1  X  I04* 

®«n 

cmVsec 

1.2  x  10° 

em’/jee 

1.2  x  10-° 

cmVsec 

4  XlO45 

cmVsec 

4  X  I041 

M. 

cm4 

8  x  10“ 

N 

cm** 

5  x  10° 

G 

#  photons  /  (cmVsec) 

2xi0“  -  1.3  x  10° 

Generation  Rale 
(10*  Photons/ctn1sec)  I 


a  1500 
b  750 
c  250 
d  200 
e  100 
f  2 
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Fig.  2.  Numerical  solutions  of  Eqs.(l)- 
(10).  The  rart  earth  luminescence 
intensity  versus  time  for  the  set  of 
parameters  shown  in  Table  I.  The  set  of 
computed  profiles  (normalized  to  unity 
at  maximum),  showing  the  buildup  of 
luminescence,  the  steady  state  value,  and 
the  decay  after  switching  off  (at  60  /ts) 
the  excitation  pulse. 


Time  (jxsec) 


Fig.  3  The  numerically  computed  rise 
and  decay  times  (solid  lines)  as  a 
function  of  generation  rates.  The 
experimental  values  of  rise  times 
(circles-single  exp.  f.),  (squares-double 
exp.  f.)  The  dots  stand  for  the  decay 
times  constant  of  the  dominant 
component  of  the  experimental  data 
fitted  to  double  exponential  function. 
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REI-trap  (N!),  or  the  free  electrons  (holes)  with  Auger  coefficients  Bbt,,(BBT4)  [14,15].  By 
consulting  Figs.  1 ,  we  can  now  complete  the  formal  description  of  the  model  by  deriving  the 
differential  equations  for  the  REI-trap  kinetics.  To  solve  this  system  ,  we  assumed, 
excitations  take  place  at  a  low  enough  temperatures  (  8.6  K)  that  thermal  activation  of  the 
trapped  carriers  is  negligible.  That  is,  terms  explicitly  dependent  on  temperature  were 
ignored,  and  only  trapping,  transferring  and  recombination  transitions  were  considered.  The 
numerical  solution  of  equations  was  obtained  using  the  parameters  shown  in  Table  I  by 
repeating  the  numerical  calculations  for  different  values  of  the  generation  rate.  The 
luminescence  intensities  are  proportional  to  N^,  N!,  and  N^.  Figure  2  shows  the  computed 
luminescence  rise  and  decay  as  a  function  of  the  generation  rate.  In  figure  3,  the  solid  lines 
are  the  computed  rise  and  decay  times  (  the  experimental  data  are  imposed  on  the  computed 
curves).  The  circles  shown  in  Fig.3  represent  the  experimental  rise  time  obtained  from  fitting 
to  a  single  exponential  function.  The  squares  stand  for  the  rise  time  of  the  dominant 
component  of  die  experimental  data  fitted  to  double  exponential  function.  The  dots  represent 
the  experimental  decay  time  of  the  dominant  component  of  the  double  exponential  fitting.  The 
numerically  simulated  luminescence  rise  and  decay  measurements  show  a  good  quantitative 
agreement  with  experiment  over  a  wide  range  of  generation  rates. 

For  n-type  material  the  equations  are  even  more  complicated,  similar  equations  hold  for  p- 
type  material.  Such  systems  of  coupled,  first  order  stiff  nonlinear  differential  equations 
require  specialized  numerical  integration  routines  designed  specifically  for  stiff  systems. 
The  detailed  analysis  of  the  parameters  and  their  influences  on  the  rise  time,  the  efficiency 
of  the  steady  state  luminescence,  and  decay  time  mil  be  published  elsewhere.  Finally,  the 
proposed  model  can  be  refined  by  taking  into  account  the  surface  recombination  and  carrier 
diffusion  processes,  which  may  play  important  roles  in  the  reduction  of  the  overall 
photoluminesoence  emission  . 

This  work  was  supported  by.AFOSR  grant  Nq.90-0322. 
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ABSTRACT 

The  photoluminescence,  time  resolved  spectra  and  kinetics  of  Yb  implanted  InP  samples  are 
studied  tinder  pulsed  and  CW  excitations  (above  and  below  band-gap)  at  different 
temperatures  and  excitation  intensity.  The  photoluminescence  intensity  and  decay  time  as  a 
function  of  temperature  is  explained  by  a  proposed  new  quenching  mechanism  involving  Be 
ion.  The  rise  and  decay  times  depend  on  excitation  intensity.  The  above  experimental  facts 
was  explained  uring  the  kinetics  model  developed  by  HJ.  Lozykowsld  (2].  The  numerically 
simulated  luminescence  rise  and  decay  times  show  a  good  quantitative  agreement  with 
experiment,  over  a  wide  range  of  generation  rates.  The  electric  field  InP:Yb 
photoluminescence  quenching  was  investigated  and  reported  for  the  first  time. 

INTRODUCTION 

There  has  been  increasing  interest  in  rare  earth  doped  semiconductor  compound  because 
of  their  potential  for  fight  emitting  derice  applications.  Betide  the  application  aspect,  the  rare 
earth  impurities  are  of  great  interest  from  the  scientific  point  of  view.  Scientific  interest  is 
related  to  the  uniqueness  of  optical  and  electrical  properties  of  rare  earth  impurities  in 
semiconductor  hosts.  Among  rare  earth  ions  are  those  that  upon  replacing  the  element  from 
column  m  in  LH-V  compounds  that  are  isovalent,  concerning  outer  electrons  of  RE*3  ions, 
create  isoelectronic  traps  in  m-V  semiconductors.  Hus  cooduskm  is  supported  by  the  fact 
that  the  atomic  covalent  radii  (ionic  RE*3)  for  all  rare  earth  are  bigger  than  atomic  radii  of 
Ga  and  In  that  they  replace  [2],  Pauling’s  electronegativity  of  rare  earth  elements  is  in  the 
range  of  1.1-1.25,  and  is  smaller  than  Ga  (1.81)  and  In  (1.78)  for  which  it  substitutes.  The 
rare  earth  isovalent  traps  (REI-trap)  that  we  can  call  isoelectronic  "structured"  impurities 
[1,2]  possess  die  unfilled  4f*  core  shell.  The  structured  luminescence  arises  from  intra- 
contigurational  f-f  transitions  in  the  core  of  the  isoelectronic  "structured"  impurities. 
Ytterbium  in  InP  replaces  indium  on  a  substitutional  site  P]  and  acts  as  an  isoelectronic  trap. 
It  was  originally  proposed  by  Whitney  at  al  [4]  and  confirmed  by  others  [5a,b]  that  the  Yb 
ion  creates  an  electron  trap  at  30  meV  below  the  bottom  of  conduction  band.  Recently, 
admittance  spectroscopy  [5b]  was  used  to  identify  the  electrical  activity  of  Yb  in  n  and  p  type 
InP.  It  was  found  that  Yb  in  InP  creates  a  hole  trap  at  50  meV  above  the  valence  band,  and 
an  electron  trap  at  29  meV  below  the  conduction  band. 

SAMPLES  AND  MEASUREMENT  RESULTS 

The  InP  material  was  produced  from  Fe-doped,  semi-insulating  iron  doped  InP. 
Samples  D  and  E  were  implanted  using  a  novel  high  current  metal  ion  implantation  facility 
in  which  a  metal  vapor  vacuum  arc  (MEW A)  ion  source  is  used.  The  Yb  was  implanted 
with  mean  ion  energy  *  140  keV  and  the  dose  (5-10)  x  10”  ion/cm2  (peak  at  —370  A  from 
the  surface,  concentration  1.8x10”  cm'3).  Sample  D  and  E  were  annealed  for  10  minutes  at 
600  °C  and  680  °C,  respectively.  Sample  C  was  implanted  using  ion  beam  energy  of  1  MeV 
and  dose  of  3  xlO13  ion/cm2(peak  at  1960  A  from  the  surface,  concentration  1.46x10**  cm' 
3).  The  sample  was  annealed  at  650  °C  for  15  minutes. 

The  photoluminescence  spectra  of  InP:  Yb  (sample  C)  shown  in  Fig.  1  was  recorded 
at  8.5  K.  The  emission  intensity  of  sample  C  is  the  highest  of  the  three  samples.  The 
spectrum  of  each  sample  consists  of  one  sharp  peak  at  1001  nra  and  broader  overlapping 
peaks  at  1005.5  nm  and  1007  nm  corresponding  to  the  2FW  transitions.  The  inset  of 
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Fig.  1  PL  Spectra  of  InP:Yb  (  sample  C  )  The  inset  shows  the  excitation  spectra. 
Fig.  2  The  PL  intensity  of  InP:  Yb  (at  1001  am)  as  a  function  of  excitation  power. 


Fig.  3  The  computed  luminescence  rise  and  decay  for  different  generation  rates. 

Kg.  4  The  calc,  and  experimental  rise  and  decay  times  as  functions  of  the  generation  rate. 

Fig.  1  shows  the  excitation  spectra  recorded  at  emission  lines  1001  nm  and  1005 .5  nm.  In 
Fig.  2,  the  intensity  of  1001  nm  peak  was  measured  as  a  function  of  the  excitation  power. 
The  PL  kinetics  measurements  were  investigated  at  different  excitation  intensities  and 
temperatures.  The  experiment  shows  that  the  rise  and  decay  times  are  function  of  excitation 
intensity.  The  above  experimental  facts  were  explained  using  the  kinetics  model  developed 
by  HJ.  Lozykowski  [2],  To  adequately  describe  the  kinetics  of  the  energy  transfer  from  the 
host  to  REI-trap,  the  buildup  and  decay  kinetics  of  rare  earth  luminescence  (as  the  function 
of  excitation  intensity),  we  have  to  consider  six  separate  states  of  the  REI-trap  (for  details  see 
[2]).  The  energy  transfer  processes  occur  through:  a)  transfer  of  energy  from  the  bound 
exdton  on  REI-trap  to  the  core  states,  and  b)  by  Auger  mechanism  where  the 
recombination  energy  of  the  bound  electron  with  free  hole  is  transferred  nonradiatively  to  the 
core  states.  If  the  initial  and  final  states  are  not  resonant  (in  both  mechanisms),  the  energy 
mismatch  must  be  accommodated  by  emission  or  absorption  of  phonons  [1,6].  In  SI  InP:Yb 
the  kinetics  are  described  by  a  set  of  differential  equations  which  were  solved  assuming  low 
temperature.  That  is,  terms  explicitly  dependent  on  temperature  were  ignored,  and  only 
trapping,  transferring  and  recombination  transitions  were  considered.  The  numerical  solution 
of  the  differential  equations  was  obtained  by  repeating  the  numerical  calculations  for 
different  generation  rates.  Fig.  3  shows  tire  computed  luminescence  rise  and  decay  as  a 
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function  of  the  generation  rate.  In  Fig.  4,  the  solid  lines  are  the  computed  rise  and  decay 
times  (the  experimental  data  are  imposed  on  the  computed  curve).  The  circles  shown  in 
Fig.4  represent  the  experimental  rise  times  obtained  from  fitting  to  single  exponential 
function.  The  squares  stand  for  the  rise  times  of  the  dominant  component  of  the  experimental 
data  fitted  to  double  exponential  function.  The  dots  represent  the  experimental  decay  times 
of  the  dominant  component  of  the  double  exponential  fitting.  The  numerically  simulated 
luminescence  rise  and  decay  measurements  show  a  good  quantitative  agreement  with 
experiment  over  a  wide  range  of  generation  rates. 

The  time  resolved  spectra  (TR,  Fig.  5)  was  recorded  at  temperatures  8.5K,  20K  and  40K 
for  different  time  delays  from  0.0  to  30/is  (for  8.S  K)  and  from  0.0  ps  to  25  fis  for 
temperatures  20K  and  40K.  The  time  resolved  emission  spectra  of  the  main  peak  at  lOOlnm 
(I,)  shows  changes  in  the  emission  line  width  and  also  slight  change  in  the  position  of  the 
peak.  The  positions  and  the  widths  of  the  emission  lines  recorded  at  long  time  after  the 
excitation  pulse  were  changed  compared  to  those  measured  just  after  the  excitation  pulse. 
The  excitation  energy  seems  to  move  within  the  regular  Yb3*  and  slightly  distorted  ion 
system,  and  also  between  Yb3*  and  other  accidental  impurities  where  the  energy  is  dissipated 
nonradiatively.  At  this  point  it  is  reasonable  to  introduce  the  hypothesis  of  energy  transfer 
from  Yb3+  to  iron  or  other  accidental  impurities.  The  energy  transfer  process  from  Yb3* 
to  Fe  ions  is  attributed  for  thermal  quenching  of  the  ytterbium  luminescence  intensity  and 
decay  time  decrease  with  temperature  increase .  Crystals  used  in  our  investigation  were  semi- 
insulating  InP  doped  with  Fe  (concentration  >  2xlOl<cm'3)  and  implanted  with  Yb.  Iron  is 
a  common  contaminant  in  MOCVD  Yb  doped  InP  semiconductor.  Uwai  etal.  [7]  presented 
the  SIMS  analysis  of  MOCVD  grown  InPtYb3*  that  showed  a  flat  ytterbium  incorporation 
profile  for  a  total  Yb3+  concentration  of  about  6x10“  cm'3.  The  SIMS  spectra  also  reveals 
significant  amounts  (I014  cm 3  or  higher)  of  Fe  and  Mn,  elements  that  were  not  present  in 
undoped  layers,  and  which  probably  appear  from  contaminated  ytterbium  [  Yb(CjHj)3]  source. 
Also  Williams[8] ,  found  an  evidence  for  impurities  such  as  Fe  and  Mn  which,  did  not  appear 
in  undoped  samples.  The  iron  contamination  of  InP:Yb  is  probably  very  common  to  most 
of  the  investigated  samples  in  literature. 

The  excellent  review  of  optical  properties  of  iron  doped  InP  (absorption, 
photoluminescence,  photoconductivity  and  photoluminescence  excitation  spectra)  was  recently 
published  by  S.  Bishop[9].  The  first  level  scheme  of  Fe*3  ion  in  InP  was  proposed  by  Stauss 
et  al  [10]  to  explain  the  photo  quenching  spectra  of  Fe5*  ESR  signal  with  onset  at  0.75  eV, 


Fig.  5  Time-resolved  spectra  of  lnP:Yb,  sampie(C),  for  different  time  delays. 
Fig.  6  The  energy  levels  schemes  for  Yb3+  and  (  Fe3+  -  Fe  2+)  in  SI  InP. 
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Fig.  8  The  intensity  ratio  of  the  emission  lineslOOl  nm  to  1005.5  nm. 

a  shoulder  near  1.13  eV,  and  maximum  between  1.2  eV  and  1.4  eV  .  These  thresholds  (at 
0.75  eV  and  1.13  eV)  were  attributed  to  direct  transfer  of  electrons  from  die  valence  band 
to  Fe?*  which  converts  it  to  the  ground  state  of  Fe2*  (*E)  or  to  excited  state  of  Fe?*  (*Tj) 
respectively.  Bishop  [11]  and  Leyral  [12]  observed  die  photoluminescence  of  Fe2*  (*E)  at 
0.35  eV  by  extrinsic  excitation  (PLE)  peaking  at  the  excitation  energy  between  1.3eV  and 
1 .35  eV.  In  conclusion  according  to  Bishop  et  al.[l  1]  the  extrinsic  excitation  of  the  Fe?*  (*Tj) 
photoluminescence  is  attributed  to  a  ligand-to-metal  charge  transfer  transition  in  which  a 
valence  electron  from  an  adjacent  phosphorus  atom  is  transferred  to  a  Fe?*  ion,  creating  an 
Fe?*  in  the  excited  *1^  state.  The  energy  levels  scheme  for  Yb3*  (4f3  in  cubic  symmetry  TJ 
and  schematic  configuration  coordinate  diagram  for  die  (Fe?*-  Fe?*)  in  SI  InP  is  shown  in 
Fig.  6.  Furthermore,  the  energy  levels  schemes  show  the  0.35  eV  luminescence  attributed  to 
the  Fe  intracenter  transition  (*T2  -»  3E),  and  energy  transfer  from  Yb3*  to  Fe?*  via  die 
state  of  Yb  ion.  Excitation  of  Fe2*  photoluminescence  at  0.35  eV  is  due  to  an  electron 
capture  by  Fe?*  giving  a  Fe2*  in  the  excited  *T2  state  which,  relaxes  radiatively  to  the  *E 
(Fe2*  ground  state).  The  structured  infrared  (0.35eV)  PL  emission  attributed  to  transition  *T2 
-*  5E  exhibits  exponential  decay.  At  10  K  the  decay  time  constant  for  semi-insulating 
sample  is  8.8  #ts  ,  and  for  n  type  sample  is  10.5  /xs  which  corresponds  to  the  lifetime  of 
the  ^Tj  excited  state  [13].  The  lifetime  of  *T2  excited  state  decreased  quickly  with  increasing 
temperature,  and  at  50  K  it  is  6.6  /ts.  The  Yb3*  -2F7/2)  emission  occurs  in  the  region 
1.229-1.242  eV  and  overlaps  the  broad  extrinsic  absorption  band  of  Fe3*  ,  a  requirement 
for  efficient  energy  transfer.  The  quenching  of  Yb*3  characteristic  PL  (integrated)  emission 
intensity  and  the  decay  time  was  investigated  as  a  function  of  temperature.  Fig.  7  shows  the 
photoluminescence  spectra  of  InP:Yb  as  a  function  of  the  temperature  in  the  range  8.5  -  80 
K.  Fig.  8  shows  the  ratio  of  the  integrated  intensity  I,  (area  under  peak  at  1001nm)to  the 
integrated  intensity  I2  (area  under  broad  peak  at  1007  nm)  as  a  function  of  temperature  for 
sample  C.  The  ratio  decreases  from  2.3  at  8.6K  to  about  1.15  at  40K  and  remains  constant 
till  80K.  Figure  9  shows  that  the  Yb3*  PL  integrated  intensity  as  a  function  of  temperature 
remains  approximately  constant  between  8.6  and  55  K,  above  which  it  decreases  rapidly  with 
further  increase  of  temperature.  The  experimental  data  (for  lOOlnm  peak)  was  fitted  to 
equation  show  in  Fig.  9  and  plotted  with  theoretical  fitting  (solid  line)  using  parameters  EA 
=  100  meV,  Eb  =7.57  meV,  a=0.9,  b  =  l. 66x10"  and  c=11.25.  The  experimental  data 
are  well  fitted  and  explained  by  the  above  simple  model  in  which  the  excited  Yb3*  ions 
lose  their  energy  via  the  2F5/2  to  the  *1^  of  Fe2*  ion  (quenching  trap).  Activation  energy  EA 
=  100  meV  obtained  from  the  experimental  data  is  in  good  agreement  with  the  gap  between 
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Fig.  9  A  fit  of  the  PL  intensity  of  InP:Yb  versus  temperature. 

Fig.  10  A  fit  of  the  experimental!  decay  times  versus  temperature. 

the  *Fi a  upper  Stark  level  of  Yb3*  and  the  excited  level  of  Fc?*.  The  energy  of  the 
excited  Ft?*  (*Tj)  ion  is  lost  via  the  *1* -• *E  photon  or  multiphonon  emissions.  A  similar  / 
simple  model  describes  the  relaxation  rate  R,  =  l/r  (decay  time)  of  die  Yb3*  luminescence 
as  a  function  of  temperature  shown  in  Fig.  10  The  relaxation  rate  (1/r)  as  a  funcuun  of 
temperature  is  given  by  die  equation  shown  in  Fig.  10  where  the  activation  energy  Eg 
=  100  mcV,  is  approximately  the  separation  between  2Fy2  (Yb3+)  and  *T2  levels,  and  a  = 
6.19x10s  sec4.  The  r,  =  11.75  /tsec  is  the  decay  time  of  Yb3*  at  very  low  temperature.  The 
lifetime  of  *Tj  (Fe3*)  excited  state  decreased  quickly  with  increasing  temperature  and  at  50 
K  it  is  6.6  fis  (13].  At  all  temperatures  investigated,  the  Ft?*  (^Tj)  lifetime  is  shorter  than 
the  Yb3+  lifetime,  therefore  Ft?*  is  appropriate  to  serve  as  an  energy  sink  for  Yb3*  excitation. 
As  mentioned  above,  the  energy  transferred  from  Yb3  to  Ft?*  ion  is  lost  via  the  *r2  ->  *E 
(0.35eV)  photon  (  or  multiphonon)  emission.  The  temperature  dependence  of  the  PL  intensity 
ofthesT1->  SE  transition  should  show  some  increase  at  the  temperature  where  the  Yb3* 

> 2Fm )  emission  intensity  drastically  decreases  (experiment  in  progress).  If  the 
outcome  of  the  experiment  will  be  positive  it  will  give  an  additional  direct  proof  of  the 
proposed  new  quenching  mechanism. 

The  width  of  the  emission  line  1001  nm  as  a  function  of  temperature  is  shown  in  Fig. 
11.  The  line  shape  in  the  temperature  range  9-80  K  is  described  by  a  convolution  of  a 
Loren tzian  and  Gaussian  functions  (the  Voigt  contour).  The  emission  line  is  homogeneously 
broadened  with  the  Lorentzian  as  the  dominant  component  (-90%)  of  the  line  shape  .  The 
width  of  the  experimental  homogeneous  component  as  a  function  of  temperature  is  plotted 
together  with  theoretical  curve  using  the  equation  shown  in  Fig.  12.  The  equation  in  Fig.  12 
is  typical  for  a  single-phonon  relaxation  process  with  absorption  of  a  phonon  of  energy  AE. 
The  fitting  parameters  are:  AE  =  5.976  meV  is  the  energy  difference  between  the  level  5  and 
4  [3],  a=0.65meV  and  b  =  11.44.  The  estimated  energy  AE  =  5.976  meV  is  dose  to  theTA 
zone  boundary  phonons  energy  TA(L)  =  6.82meV  (  TA(X)  =  8.478  meV  ),  in  InP  (14]. 

The  effect  of  the  electric  field  on  the  PL  of  InP:Yb  was  studied  for  the  first  time.To 
perform  this  experiment,  a  semi-transparent  gold  electrode  (100  A)  was  deposited  on  the 
implanted  surface  of  InP,  a  second  thick  gold  electrode  (500  A)  was  deposited  on  the 
opposite  side  of  the  sample.  Figure  12  shows  the  quenching  of  PL  of  Yb3*  emission  as  a 
function  of  the  applied  electric  field  at  both  polarities.  Assuming  uniform  absorption  ,  the 
excitation  photons  are  absorbed  within  the  surface  layer  where  a  strong  electric  field  exists. 
The  quenching  of  PL  of  Yb3*  can  be  explained  by  two  possible  mechanisms  :  (1)  exciton 
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Fig.  11  The  temperature  dependence  of  the  line  width  of  the  emission  line  1001  nm. 
Fig.  12  The  PL  intensity  of  InP:  Yb  as  a  function  of  the  applied  DC  voltage. 


bound  to  the  isoelectronic  trap  dissociation  by  hot  election  impact  ionization  or  (2)  swept  of 
electrons  or  holes  to  the  Au-InP  interface  where  they  will  recombine  nonradiatively. 
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ABSTRACT 

The  photoluminescence  (PL)  and  time  resolved  spectra  of  Nd-  and  Yb-implanted  CdS  samples 
are  studied  under  pulsed  and  CW  excitations  using  Ar*  ion  laser  for  different 
intensities  and  temperature  (9  -  300  K).  The  samples  were  annealed  under  different  conditions 
using  the  thermal-pulse  method.  For  CdS:Nd,  the  PL  spectra  was  recorded  in  the  range  890- 
930  nm  (transitions  and  for  CdStYb  it  was  recorded  in  the  range  985  -  1010  nm 

(transitions  -*  2F7/2).  The  Rise  time  and  decay  time  were  studied  for  both  CdS.'Nd  and 
CdS:Yb  at  different  emission  lines  as  a  function  of  temperature  and  excitation  power. 

INTRODUCTION 

The  photoluminescence  properties  of  rare  earth  (RE)  ions  in  II-VI  semiconductor 
compounds  have  been  of  interr  for  many  years.  The  light  emission  due  to  radiative 
transitions  in  the  partially  occupied  4f  shell  of  RE3*  ions  in  the  form  of  lines  and  narrow 
bands  is  the  basis  for  numerous  applications  such  as  solid  state  lasers.  Beside  the  application 
aspect,  the  rare  earth  impurities  are  of  great  interest  from  scientific  point  of  view.  Scientific 
interest  is  related  to  the  uniqueness  of  optical  and  electrical  properties  of  rare  earth  impurities 
in  semiconductor  hosts.  It  is  well  known  that  the  rare  earth  luminescence  depends  very  tittle 
on  the  nature  of  the  host  and  the  temperature.  The  4f  orbits  of  rare  earth  ions  incorporated 
in  semiconductors  are  so  deeply  buried  within  the  electronic  shell  that  the  energy  levels  of 
the  4P  configuration  are  only  slightly  perturbed  compared  to  free  ion  energy  levels.  The 
electronic  structure  of  the  rare  earth  luminescence  centers  and  their  electrical  activities,  as 
well  as  their  indirect  photoluminescence  excitation  mechanisms  are  still  not  well  understood. 
The  luminescence  of  rare  earth  doped  cadmium  sulphide  has  been  reported  by  many 
researchers[l-4].  The  attempt  was  to  explain  the  mechanisms  of  exciting  the  RE?*  in  CdS. 
In  [1],  Kingsley  and  et  al,  proposed  that  the  energy  is  transferred  from  luminescence  centers 
( in  his  case,  copper)  to  the  rare  earth  ions.  Anderson  [2]  proposed  that  the  transfer  of  energy 
is  due  to  the  electronic  transition  between  donor-acceptor  levels.  In  this  paper,  we  report  the 
luminescence  properties  of  Nd3*  and  Yb3*  implanted  in  CdS.  Time  resolved  spectra  was 
recorded  for  Nd3*  and  Yb3*  in  the  ranges  890-930  nm  (transitions  *F3/7-**L)a)  and  985- 
1010  nm  (transitions  2F5/2->2F7/2))  respectively.  The  term  assignments  for  these  transitions  are 
obtained  from  (5],  The  PL  spectra  was  investigated  at  different  temperatures  in  the  range  8.5 
-  300  K.  Rise  time  and  decay  time  processes  were  also  studied  as  a  function  of  temperature 
and  excitation  intensity. 

SAMPLES  AND  MEASUREMENT 

The  implantation  of  Yb  and  Nd  was  performed  with  an  unconventional  technique  using 
a  new  kind  of  high  current  metal  ion  source.  The  implantation  dosage  was  7.5  x  10‘3  cm'2 
and  the  mean  ion  energy  was  150  keV  for  Nd  and  140  keV  for  Yb.  In  the  CdS:Nd  sample, 
the  concentration  of  Nd  is  1.31  X 10'9  cm'3  and  the  peak  -to-surface  depth  is  412  A.  For  Yb, 
the  concentration  is  1.38  xlO19  cm'3  and  the  peak-to-surface  depth  is  390  A.  Both  samples 
were  annealed  at  900  °C  for  15  seconds  using  a  Thermal-Heat-Pulse  method.  The 
photoluminescence  (PL)  was  excited  by  the  488  nm  line  from  a  argon  ion  laser  (Laser  Ionics 
Model  553A).  The  pulsed  excitation  for  time  resolved  spectra  (TRS)  was  obtained  by 
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acousto-optic  light  modulator  (IntraAction  Model  AOM-125).  The  samples  were  mounted  on 
a  cold  finger  cooled  by  a  close-cycle  helium  cryostat  down  to  8.5  K.  The  emission  was 
dispersed  by  a  Jarrel  Ash  Model  78-490,  0.75  M,  scanning  monochrometer  equipped  with 
a  1180  grooves/mm  grating.  The  detecting  electronics  consists  of  a  Hamamatsu  R632-01 
photomultiplier  with  spectral  response  extended  to  1R  region. The  signal  from  the 
photomultiplier  is  amplified  by  a  fast  preamplifier  and  then  fed  to  a  dual-channel  gated 
photon  counting  system  controlled  by  a  computer,  which  also  control  the  scanning 
monochromator.  The  photoluminescence  kinetics  (rise  and  decay)  were  measured  with 
exceptional  accuracy  using  Turbo-Multi  Channel  Scaler  (Turbo-MCS,  EG&G  Ortec)  which 
offers  a  wide  range  of  channel  dwell  times  (minimum  5ns),  with  no  deadtime  between 
channels. 

RESULTS  AND  DISCUSSION 

The  Photoluminescence  spectra  of  CdS:Nd  and  CdSrYb  was  recorded  at  8.5  K  in  the 
spectral  range  895  -  930  nm  (1.385-1.335  eV)  and  985-1010  nm  (1.26-1.225  eV), 
respectively.  For  CdS:Nd,  the  spectrum  consists  of  about  nine  sharp  peaks  corresponding 
to  the  transitions  between  Stark  levels  *FW  -  4I»j  .  The  PL  spectrum  of  CdS:Yb  consists  of 
six  sharp  peaks  around  988  nm  (1.255  eV)  and  two  broader  peaks  at  995.6  nm  (  1.245  eV) 


1.34  1.36  1.38  1.23  1.24  1.25 


Photon  Energy  in  eV  Photon  Energy  in  eV 

Figure  1.  The  Time  Resolved  Spectra  f"RS)  of  (a)  CdS:Nd  and  (b)  CdS:Yb  at  T=8.6  K. 
The  gate  widths  are  100  nsec.  The  de  s  are  as  indicated  in  each  spectrum.  The  spectra 
labeled  ON  are  measured  with  gate  and  excitation  pulse  overlapped. 
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and  5)98.8  nm  (1.241  eV).  The  Yb3*  emission  peaks  correspond  to  the^Fj/j  -  2F7/2  transitions. 
The  time  resolved  spectra  were  recorded  at  8.6  K  for  different  delay  times  from  0  to  1000 
nsec  for  Nd  and  from  0  to  120  nsec  for  Yb  as  shown  in  figure  1.  For  CdS:Nd,  the  intensity 
of  all  emission  lines  decreased  with  delay  uniformly.  In  the  case  of  CdStYb  the  intensities 
of  the  emission  lines  at  986  nm  (1.257  eV),  995.6  nm  (1.245  eV)  and  998.8  nm  (1.241  eV) 
decreased  slower  than  the  other  lines. 

The  PL  spectra  were  investigated  as  a  function  of  temperature  in  the  range  8.5  -  300  K  for 
Nd  and  8.5-  150  K  for  Yb.  In  figure  2a,  the  intensity  of  all  lines  of  CdS:Nd  stayed  fairly 
constant  from  8.5  to  35  K  then  fill  down  until  no  emission  was  observed  after  100  K. 


Figure  2.  Temperature  dependence  of  the  PL  spectra  of  (a)  CdS:Nd  and  (b)  CdS:Yb 
excited  by  Ar+  laser  at  488  nm  with  80  mW  power. 

The  integrated  intensity  of  the  emission  line  903  nm  (1.372  eV)  is  plotted  versus 
temperature  as  shown  in  figure  3a.  From  figure  2b,  the  PL  intensity  of  Yb  sharp  emission 
peaks  decreases  with  the  temperature  increase  until  40  K  and  started  to  increase  with 
temperature  until  it  reached  a  maximum  at  around  55  K.  Then  started  to  decrease  with  further 
increase  in  temperature.  This  behavior  is  clear  from  figure  3b  where  the  integrated  intensity 
of  the  emission  line  995.6  nm  (1.245  eV)  is  plotted  as  a  function  of  temperature.  Another 
observation  is  that  with  the  increase  of  temperature,  the  sharp  peaks  evolved  into  one  broad 
emission  band  at  the  same  range  of  wavelength. 

In  figure  3a  the  decay  characteristic  of  the  Nd3*  4F3^  integrated  emission  intensity  is  shown. 
The  experimental  data  of  the  temperature  dependance  (  at  peak  903  nm  )  was  fitted  to  the 
relationship: 
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I(T)  =J(0)|l+A1e  *T 


The  experimental  result  is  plotted  with  theoretical  fitting  (solid  line)  using  the  parameters 
E,=72.4  meV,  1(0)  =  0.99  and  A,  =  3.5  xlO4. 


Figure  3.  PL  integrated  intensity  of  (a)  CdS:Nd  at  emission  line  903  nm  and  (b)  CdS:  Yb 
at  emission  line  995.6  nm  versus  temperature.  For  (a)  CdS:Nd  the  solid  line  is  a  fit  to  the 
experimental  data  using  the  relation  discussed  in  the  text. 

The  widths  of  the  Nd  emission  line  at  903  nm  (1.372  eV)  and  Yb  emission  line  at 
995.6  nm  (1.245  eV  )  as  functions  of  temperature  are  shown  in  figure  4.  The  line  shapes  of 
both  samples  were  described  by  a  Gaussian  (  inhomogeneously  broadened  )  function  in  the 
foremen tioned  temperature  range.  The  experimental  line  widths  for  both  Nd  and  Yb  as 
functions  of  temperature  shown  in  figure  4  (a  &  b,  respectively)  are  plotted  together  with  the 
theoretical  curves  using  equation  : 
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(t)  «  a  + 
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a  =  0.269  meV 
b  =  5.205  meV 
iE  =  30.254  meV 
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Figure  4.  The  temperature  dependence  of  line  widths  of  (a)  CdS:Nd  at  903  nm  and  (b) 

CdS:Yb  at  995.6  nm. 
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For  Nd,  AE=8.3  meV,  a =0.343  meV  and  b =0.351  meV,  and  for  Yb,  AE=30.254  meV, 
a=0.269  meV  and  b  =5.205  meV.  The  above  equation  describing  the  width  of  the  emission 
line  as  a  function  of  temperature,  is  typical  for  a  single-phonon  relaxation  process  with 
absorption  of  a  phonon  of  energy  AE=8.3  meV  (ErA=9.0meV  in  CdS)  for  Nd,  and 
AE=30.254  meV  (^=30.0  meV  in  CdS)  for  Yb. 

The  rise  time  and  decay  time  of  Nd  and  Yb  were  investigated  as  functions  of 
temperature  and  excitation  power.  Figure  5  shows  the  decay  time  of  Nd  at  emission  peak  903 
nm  (1.372  eV)  as  a  function  of  temperature. 


Figure  5.  A  fit  (solid  line)  of  the  time  decay  CdS:Nd  at  emission  line  903  nm  to  the 
experimental  data  (crosses)  as  a  function  of  temperature. 

The  decay  time  remains  constant  around  357  fiscc  from  8.6  to  50  K  and  after  it  falls  down 
below  100  nsec  at  100  K.  The  decay  time,  ,  as  a  function  of  temperature  is  given  by  the 
equation 


where  Ea  =  84.34  meV,  =  7.52  n sec  and  t,  =  357.92  //sec  where  r,  is  the  decay  time 
at  very  low  temperature.  The  activation  energy  EA  =  84.34  meV  is  close  to  the  activation 
energy  E,  =  72.4  meV  obtained  from  the  decay  characteristics  of  the  integrated  intensity  as 
a  function  of  temperature.  For  Yb,  the  decay  time  of  the  emission  line  995.6  nm  (1.245  eV) 
remained  constant  around  90  //sec  for  the  entire  range  of  temperature  of  8.6  -  100  K. 

In  figure  6,  The  rise  times  of  both  Nd  and  Yb,  were  studied  as  functions  of  excitation 
power  in  the  range  1  -  400  W/cm2.  The  rise  time  of  Nd  at  903  nm  (  1 .372  eV  )  decreased 
from  556  //sec  to  300  //sec  as  the  excitation  power  was  increased  from  1  W/cm2  to  400 
W/cm2  as  shown  in  figure  6a.  for  Yb  (figure  6b).  the  rise  time  at  988.2  nm  (  1.254  eV  ) 
decreased  from  97  //sec  to  23.5  //sec  as  the  excitation  power  was  increased  from  1  W/cm2 
to  400  W/cm2.  The  rise  time  decrease  w-ith  the  excitation  power  increase  indicates  that  the 
mechanism  of  exciting  the  RE  is  indirect. 

The  PL  intensities  of  CdStNd  (  at  emission  line  903  nm  )  and  CdSrYb  (  at  emission 
line  988.2  nm)  were  recorded  at  different  excitation  powers  in  the  range  1-400  W/cm2.  From 
figure  7,  we  can  see  that  the  intensities  increase  linearly  w-ith  the  excitation  power  until  100 
W/cm2  and  then  start  to  saturate. 

In  order  to  interpret  our  results  more  thoroughly  and  provide  explanations  foi  rgy 

transfer  and  kinetics  processes,  more  investigation  is  needed,  such  as  j  .on 

measurements,  excitation  spectra,  and  the  electroluminescence. 
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Figure  6.  The  rise  times  of  (a)  CdS:Nd  at  903  nm  and  (b)  CdS:Yb  at  988.2  nm  as 
functions  of  the  laser  excitation  power. 


Excitation  Power  (W/cm1)  Excitation  Power  (K/ctn1) 


Figure  7.  The  PL  intensities  of  (a)  CdS:Nd  at  903  nm  and  (b)  CdS:Yb  at  988.2  nm  as 
functions  of  the  laser  excitation  power 
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Digital  spectropolarimeter  for  the  measurement  of  optical  polarization 
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A  computerized  setup  for  optical  polarization  measurements  is  described.  The  heart  of  this  setup 
is  a  photoelastic  modulator  (PEM)  whose  retardation  phase  can  be  kept  constant  when  the 
wavelength  of  modulated  light  is  scanned  in  a  wide  range,  and  a  dual-channel  gated  photon 
counter  which  has  the  capability  of  compensating  the  variation  of  excitation  intensity  with  time 
and  wavelength.  Unlike  conventional  analog  setups,  where  the  difference  of  two  perpendicularly 
polarized  light  components  (/+  —  /“)  and  the  average  of  these  components  j(/+  +  /“}  are 
extracted  electronically  and  the  ratio  2(/+  —  /-)/(/+  +  /“)  is  recorded  as  degree  of 
polarization,  this  spectropolarimeter  directly  records  the  /+  and  I~  components  (for  circular 
polarization)  or  /g  and  IL  components  (for  linear  polarization)  of  the  analyzed  light  signal. 

Once  the  spectra  of  the  two  components  are  stored  as  a  numerical  file  in  the  computer,  they  can 
be  processed  to  extract  the  degree  of  polarization  as  a  function  of  wavelength,  temperature,  etc. 

The  curvature  and  asymmetry  of  the  modulator  response  are  taken  into  account  exactly  in  data 
processing.  The  spectropolarimeter  has  proved  to  be  very  accurate  due  to  the  capability  of 
compensating  the  fluctuation  in  excitation  source  and  of  electronically  tracking  the  modulation 
wavelength. 


I.  INTRODUCTION 

Optical  orientation  (or  optical  pumping)  measurement 
has  been  a  very  useful  tool  in  physical,  chemical,  and  bio¬ 
chemical  resea  ch1'3  because  of  its  sensitivity,  selectivity, 
and  accuracy.  In  this  article,  we  describe  a  computerized 
setup  for  optical  polarization  measurements.  This  setup 
consists  of  a  photoelastic  modulator  (PEM)  whose  retar¬ 
dation  phase  can  be  kept  constant  even  when  the  wave¬ 
length  of  the  modulated  light  is  scanned  in  a  wide  range, 
and  a  dual-channel  gated  photon  counter  which  can  com¬ 
pensate  the  variation  of  excitation  intensity  with  time  and 
wavelength.  Unlike  conventional  analog  setups,  where  the 
difference  of  two  perpendicularly  polarized  light  compo¬ 
nents  (/+  —  /“)  and  the  average  of  these  components 
i(I+  +  I~ )  are  extracted  electronically  and  the  ratio 
2(/+  —  /“ )/(/+  +  /“)  is  recorded  as  degree  of  polariza¬ 
tion,  this  spectropolarimeter  directly  records  the  /+  and 
1~  components  (for  circular  polarization)  or  /g  and 
components  (for  linear  polarization)  of  the  analyzed  light 
signal  in  the  same  run  of  wavelength,  temperature,  etc.  The 
spectra  of  the  two  components  are  stored  as  a  numerical 
file  in  a  computer  and  can  be  further  processed  to  extract 
the  degree  of  polarization  as  a  function  of  wavelength,  tem¬ 
perature,  etc. 

II.  PRINCIPLE 

The  circular  or  linear  polarization  of  an  optical  signal 
can  be  analyzed  by  inserting  a  A/4  retarder  and  a  linear 
polarizer  (or,  for  linear  polarization,  simply  a  linear  polar¬ 
izer)  between  the  source  and  the  detector.  To  measure  the 
degree  of  circular  polarization,  one  alternately  sets  the  lin- 
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ear  polarizer  +45*  and  —45*  with  the  fast  axis  of  the  re¬ 
tarder  so  that  <7+  and  a~  components  of  the  signal  can  be 
detected  alternately.  For  linear  polarization,  the  retarder  is 
removed  and  the  linear  polarizer  is  set  +45*  and  —45*  to 
the  horizontal  so  that  /g  and  IL  can  be  measured  alter¬ 
nately.  The  degree  of  polarization  is  defined  as 

/+_/- 


For  circular  polarization,  /+  is  the  intensity  of  right-hand 
circular  polarization  (o+ )  and  I~  is  the  intensity  of  left- 
hand  circular  polarization  (<r-).  For  linear  polarization, 
I+  and  I~  stand  for  the  beam  polarized  along  the  +45* 
and  —45*  axes,  respectively. 

The  automatic  spectropolarimeters  utilizing  the  me¬ 
chanically  rotating  polarizer  and  analog  lock-in  techniques 
are  described  in  many  publications.4,5  There  are  two  main 
drawbacks  of  this  technique.  First,  the  response  of  devices 
such  as  the  monochromator’s  grating  and  the  photomulti¬ 
plier  depends  on  the  polarization  of  the  incoming  light,  and 
this  can  induce  error  in  the  measured  degree  of  polariza¬ 
tion  [Eq.  (1)].  The  second  is  the  nonachromaticity  of  the 
retarder.  To  overcome  these  limitations,  a  fixed  analyzer 
and  a  modulated  retarder  are  put  in  series  with  their  opti¬ 
cal  axes  45°  apart.  The  retardation  is  modulated  such  that 
the  retardance  A  is  a  function  of  time  given  by 

A=A0  sin(2irvf),  (2) 

where  A0  is  proportional  to  the  driving  voltage  and  v  is  the 
modulation  frequency.  Now  that  a  fixed  analyzer  is  used, 
the  photomultiplier  as  well  as  the  grating  of  the  monochro¬ 
mator  are  exposed  to  only  one  state  of  linear  polarization 
and  hence  the  dependence  of  their  response  on  the  polar¬ 
ization  of  incoming  light  has  no  effect  on  the  accuracy  of 
the  measurements. 
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FIG.  1.  Basic  setup  for  analyzing  the  circular  polarization  of  light  /,  and 

/}  are  the  intensity  of  the  light  before  and  after  the  modulator,  respec-  p\Q  y  The  theoretical  response  of  the  P  EM -analyzer  system  to  linearly 

tively.  /  is  the  intensity  seesi  by  the  detector,  n,  and  n,  ate  the  refraction  polarized  light  The  solid  line  is  for  /(+4S*),  the  dashed  line  is  for 

indices  in  the  x  and  y  directions.  7(-45*),  and  2fi  is  the  gate  width  for  the  counters  A  and  B. 


To  assure  a  constant  retardation  for  different  wave¬ 
lengths,  an  electronic  compensation  system  was  designed. 

Figure  1  illustrates  this  setup.  The  standard  electro¬ 
magnetic  field  complex  vector  analysis  given  by  Kemp6  is 
followed  here.  The  transmitted  intensities  7+  and  I~  for 
a*  and  a~  circularly  polarized  light,  respectively,  are 
given  by 

/*  =  |a|2{l±sin(^0sin  2wvf)],  (3) 

where  ( + )  stands  for  ct+  and  ( — )  for  o~  .  To  see  7+  and 
l~  alternately  at  a  given  wavelength,  one  controls  the  driv¬ 
ing  voltage  of  the  PEM  modulator  to  obtain  A0  =  n/2  at  all 
wavelengths.  Hence  the  PEM  acts  like  a  +  A/4  and  —A/4 
plate  alternately.  Accordingly,  Eq.  (3)  becomes 

7±  =  |a|2[l±sin(ir/2sin2irvr)].  (4) 

Figure  2  shows  the  plot  of  Eq.  (4)  where  the  solid  line 
is  for  ct+  and  the  dashed  line  is  for  a~ .  It  is  clear  from  Fig. 
2  that  Eq.  (4)  has  symmetric  features  and  if  the  photon 
counters  A  and  B  are  gated  as  shown  in  the  figure,  the 
degree  of  polarization  will  be 


A-B 
P~A+B ' 


(5) 


For  an  incident  beam  linearly  polarized  along  ±45*  to 
the  vertical,  the  transmitted  intensities  are  given  by 


1.2 
i 

0.8 
0.6 
0.4 
0.2 
0 

o  */2  *  3x/2  2n 

Phase  of  Driving  Voltage 

FIG.  2.  The  theoretical  response  of  the  PEM-analyzer  system  to  circu¬ 
larly  polarized  light.  The  solid  line  is  for  a +,  the  dashed  line  is  for  a~,  and 
25  is  the  gate  width  for  the  counters  A  and  B. 


/(±45*)=M2[l±cos(4,sin2irvf)].  (6) 

To  see  7(±43*)  beams  alternately  at  a  given  wave¬ 
length,  we  control  the  driving  voltage  of  the  PEM  modu¬ 
lator  to  obtain  Aq  —  v  at  all  wavelengths  and  hence  the 
PEM  acts  alternately  as  an  achromatic  +A/2  and  —A/2 
retarder.  Accordingly,  Eq.  (6)  becomes 

7(  ±45*)=  |a|2[l  ±cos(irsin  2wvr)].  (7) 

The  plot  of  Eq.  (7)  is  shown  in  Fig.  3  where  the  solid  line 
is  for  7(  +45*)  and  the  dashed  line  is  for  7(  —45*).  It  can 
be  seen  from  the  figure  that  the  two  curves  show  asymmet¬ 
ric  features.  In  the  conventional  lock-in  technique,  addi¬ 
tional  electronic  circuits  are  needed  to  process  this  asym¬ 
metric  wave  form,  which  also  requires  that  the  PEM  be 
switched  on  and  off  (enabled  and  disabled).3  In  our  sys¬ 
tem,  we  deal  with  this  marked  asymmetry  numerically. 
Suppose  that  photon  counters  A  and  B  are  gated  as  shown 
in  Fig.  3  and  the  incident  beam  has  two  components 
7(  +45*)  and  7(  —45*).  From  Eq.  (7),  the  counts  in  A  and 
B  are  written  as 

/I  f*+s 

A=r)J-I(+45')  J  [1  +cos(ir sin  x)]  dx 

1  pr+5  \ 

+-7(-45‘)  J  [l-cos(*rsin x)]dxj  (8) 

(1  rv/2+6 

-7(+45*)  [  1  +cos(ir  sin  x)]  dx 

l  Jv/2-6 

1  r<r/2  +  S  \ 

+-7(-45')  [1— cos(irsinx)]  dx  1  , 

1  J  v/2—6  1 

(9) 

where  2 8  is  the  gate  width  and  tja  and  t)B  are  effective 
system  efficiencies  which  are  proportional  to  the  transmis¬ 
sivity  of  the  monochromator,  the  efficiency  of  the  photo¬ 
multiplier,  the  gain  of  the  preamplifier,  and  the  counter’s 
threshold.  In  the  system,  A  and  B  are  actually  using  the 
same  detecting  system  (monochromator,  photomultiplier, 
and  preamplifier),  and  if  the  detection  thresholds  of  the 
counters  are  matched: 
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HO.  4.  Block  diagram  of  the  spcctropolarimcter.  PI,  PI  are  linear  po¬ 
larizers,  it  is  an  achromatic  XJ4  retarder.  PEM  is  a  photoelastic  modu¬ 
lator,  and  PMT1  and  PMT2  are  photomultiplier  tubes. 

Va=Vb=V-  (10) 

Normally  what  we  are  interested  in  is  the  ratio  of  the 
signal;  so  in  the  following  we  are  going  to  drop  this  factor, 
■q.  Using  the  parity  of  the  function  cos(ir  sin  jc)  and  ignor¬ 
ing  17,  Eqs.  (8)  and  (9)  can  be  written  as 


B=6[I(  +45*)  +/( -45*)  ]  +£(/(  +  45*)  — /( -45*)  ]. 


where 


““  f 
-L 


cos(rr  sin  x)dx. 


cos(ir  sin  x)dx. 


Equations  (11)  and  (12)  can  be  put  in  a  matrix  form: 

\S  al  r/(+4y)+/(-45*)|_f^| 

[8  0]  [/(+45*)-/(-45*)|-[*|’  (I: 

which  gives 


/(+45*)-/(-4r)  1 8  B 

p  /(+45*)-f/(-45*)~  A  o|'  (16) 

B  0\ 

Equations  (13)  and  ( 14)  can  be  evaluated  numerically  for 
a  given  8. 

Similar  treatment  should  be  made  for  the  circularly 
polarized  case  especially  when  the  signal  is  weak  and  a 
wide  gate  is  needed.  In  that  case 

6  A  —B 

o-aTb’  (17) 


where 


i^-sinxj  dx. 


-4=5l/(+45*)+/(— 45*)  ]+<*[/( +45*)— /(— 45*)  ],  When  the  signal  is  strong  and  8  can  be  very  small,  y~8 

(11)  and  Eq.  (17)  simplifies  to  Eq.  (5). 


j 


FIG.  5.  Setup  used  to  test  the  performance  of  the  PEM  modulator-analyzer  system.  The  polarization  state  of  the  incident  beam  is  defined  by  the 
orientation  axis  of  the  polarizer  and  the  Fresnel  rhomb. 
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FIG.  6.  Measured  response  of  the  spectropolarimetcr  to  two  input  mono¬ 
chromatic  beams  (A=* 632.8  nm)  circularly  polarized  a*  and  o~.  The 
solid  line  is  the  response  to  a*  and  the  dashed  tine  to  a~  circular  polar¬ 
ization.  The  gate  width  and  step  size  are  set  to  0.2  /is. 

III.  INSTRUMENTATION 

Figure  4  shows  the  block  diagram  of  the  setup.  The 
monochromatic  light  from  the  Ar+  ion  laser  (or  Tirsap- 
phire  laser  or  dye  laser  pumped  by  Ar+  laser),  which  is 
already  linearly  polarized  by  the  Brewster  output  window, 
is  retarded  by  the  A/4  achromatic  retarder  (R)  to  get  cir¬ 
cularly  polarized  excitation  light.  For  an  experiment  that 
requires  linear  polarization,  one  can  instead  use  the  origi¬ 
nal  laser  beam  or  a  A/2  achromatic  retarder.  If  the  excita¬ 
tion  light  is  from  a  lamp,  a  monochromator  and  a  linear 
polarizer  (P, )  are  to  be  used  in  front  of  the  quarter  (half) 
wave  achromatic  retarder  (JR).  For  a  linearly  polarized 
excitation  light,  the  quarter  wave  retarder  is  removed  and 
a  A/2  achromatic  retarder  is  used  to  rotate  the  polarization 
plane.  The  samples  are  mounted  in  a  closed  cycle  cold 
finger  cryostat.  Small  angles  between  the  excitation  beam 
and  the  quantization  axis  of  the  sample,  and  between  this 
latter  and  the  direction  of  observation  (detection)  are  pre¬ 
ferred  for  circularly  polarized  optical  pumping  studies.7 
Different  arrangement  of  the  excitation  beam  and  the  sam¬ 
ple’s  orientation  have  been  used  for  linear  polarization 
studies.8,9  The  emission  (or  optical  signal  such  as  transmis¬ 
sion,  reflection,  etc.)  from  the  sample  is  modulated  by  a 


FIG.  7.  Measured  response  of  the  spectropolari meter  to  two  input  mono¬ 
chromatic  beams  (X =632.8  nm)  linearly  polarized  /(+45*)  and 
/( —45*).  The  solid  line  is  the  response  to  /( +45*)  and  the  dashed  line  to 
f(  — 45*)  linear  polarization.  The  gate  width  and  step  size  are  set  to  0.2 
MS 
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FIG.  8.  Geometry  used  to  simulate  the  partially  polarized  light.  The 
linear  polarizer  is  rotated  f>"  away  from  45*  to  the  fast  axis  of  the  retarder. 

photoelastic  modulator  (PEM)  whose  modulation  depth 
(retardance)  is  electronically  controlled  to  follow  the 
wavelength  of  the  signal  being  measured.  A  second  polar¬ 
izer  (P2)  is  used  to  analyze  the  modulated  light.  The  light 
emerging  from  P2  is  then  dispersed  by  a  monochromator 
and  detected  by  a  photomultiplier  (PMT1).  The  electric 
signal  from  PMT1  is  amplified  by  a  fast  preamplifier  and 
then  fed  into  a  dual-channel  gated  photon  counter  which  is 
capable  of  compensating  the  fluctuation  of  excitation  in¬ 
tensity.  To  sense  the  variation  of  excitation  intensity,  the 
incident  beam  is  split  and  coupled  into  a  second  photomul¬ 
tiplier  (PMT2)  through  the  optical  fiber.  The  output  of 
PMT2  is  preamplified  and  then  fed  into  a  photon  counter 
to  determine  the  counting  period. 

In  the  photon  counter  (SR-400)  system,  there  are 
three  fast  counters:  A,  B,  and  T,  all  of  which  operate  at 
rates  up  to  200  MHz.  Counter  T  is  presetable  and  deter¬ 
mines  the  counting  period  for  counters  A  and  B  which  may 
be  synchronized  to  external  events  via  two  independent 
gate  generators.  The  gate  generators  provide  gates  from  5 
ns  to  1  s  in  duration  with  a  delay  from  an  external  trigger 
ranging  from  25  ns  to  1  s.  The  counting  period  is  the  time 
it  takes  the  counter  T  to  reach  its  present  count  If  N,  is  the 
preset  count  in  counter  T,  and  and  R,  are  the  mean 

pulse  rates  of  the  inputs  to  counters  A,  B,  and  T,  then  the 
counting  period  is 


FIG.  9.  Comparison  of  the  theoretical  curve  calculated  from  Eq.  (31) 
(solid  line)  and  the  experimental  measurements  for  partially  polarized 
light.  (*)  for  linear  polarization  and  ( A )  for  circular  polarization. 
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FIG.  10.  Photolnmiiwacmct  of  GVmZno^Tc/CdTe  strained-layer  quan¬ 
tum  wdls.  (a)  Lm  «=  25  A.  (b)  Lm  «  30  A,  (c)  Lm  =  75  A. 

r=AT//i<.  (19) 

The  number  of  counts  in  counters  A  and  B  will  be 
N, 

a=**-r<  (20) 

N, 

B=Rb-g-  (21) 

Since  N,  is  a  constant  in  the  above  equations,  the  number 
of  counts  in  counters  A  and  B  are  proportional  to  the  ratio 
of  the  rates  Ra/R,  and  Rt/R„  respectively.  If  the  input  of 
counter  T  monitors  the  excitation  intensity,  the  counting 
period  will  compensate  for  the  fluctuation  in  the  excitation 
source. 

Once  counter  T  reaches  its  present  count,  all  the 
counters  are  "stopped”  and  the  computer  is  instructed  to 
read  the  data  and  adjust  the  driving  voltage  of  the  PEM  as 
well  as  the  wavelength  (depending  on  the  application, 
other  parameters  may  be  adjusted  during  this  time).  When 
everything  is  stabilized,  all  counters  are  reset  and  the  next 
counting  period  is  initiated.  All  other  components  in  the 
setup  are  standard. 


IV.  PERFORMANCE 

Figure  5  shows  the  test  setup.  The  spectropolarimeter 
was  tested  using  the  He-Ne  laser  632.8-nm  line  which  was 
100%  circularly  polarized  by  the  linear  polarizer  and  the 
achromatic  Fresnel  rhomb.  The  plots  in  Fig.  6  were  re¬ 
corded  when  the  linear  polarizer  was  set  at  +45'  and  —45' 
from  the  fast  axis  of  the  retarder  (<r+  and  a~ )  and  the  gate 
of  counter  A  is  scanned  with  0.2  /us  in  gate  width  and  delay 
increment.  The  solid  (dashed)  line  shows  the  intensity  of 
<r+  (<7~ )  circularly  polarized  beam  versus  the  phase  angle 
of  the  driving  voltage.  A  similar  test  was  also  done  for 
linearly  polarized  light  and  the  results  are  plotted  in  Fig.  7. 
It  is  clear  that  the  experimental  results  in  Figs.  6  and  7 
match  the  theoretical  results  shown  in  Figs.  2  and  3. 

As  shown  in  Fig.  8,  partially  polarized  light  can  be 
simulated  by  rotating  the  linear  polarizer  away  from  45*  to 
the  fast  axis  of  the  retarder.  The  degree  of  polarization 
with  a  rotated  angle  tp  can  be  written  as 


1.65  1.7  1.75 

Photon  Energy  (eV) 

FIG.  11.  Circularly  polarized  excitation  spectra  for 
CdanZngj'Te/CdTe  strained-layer  quantum  well  with  Lm  *»  75  A.  Curve 
(a)  a*,  curve  (b)  a~,  and  curve  (c)  the  degree  of  polarization  p. 

I+—I~  cos2  tp— sin1  <p 

^  I++I~  cos1  qj+sin2  tp 

=  cos2tp=sin  26  ( -45*<0<45*). 

(22) 

Equation  (22)  is  valid  for  both  circular  and  linear  polar¬ 
ization.  Partial  linear  polarization  was  tested  with  6=ir/5 
(and  hence  a =0.3 120, 0=  —0.6056)  in  Eq.  (16)  and  par¬ 
tial  circular  polarization  was  tested  with  6=ir/4  (and 
hence  y= 0.7684)  in  Eq.  (17).  The  results  are  compared 
with  Eq.  (22)  in  Fig.  9  where  the  solid  line  is  the  theoret¬ 
ical  curve  given  by  Eq.  (22).  The  superimposed  symbols 
(*)  and  (▲)  are  the  experimental  data  for  linear  and 
circular  polarization,  respectively. 

The  setup  of  the  apparatus  shown  in  Fig.  4  has  been 
used  for  photoluminescence  spectra  (PL),  excitation 
spectra  (PLE),  and  optical  pumping  studies  of 
Cd]  _  jZn/Te/CdTe  strained  layer  single  quantum  wells. 
Figures  10  and  11  show  the  PL  and  circularly  polarized 
excitation  spectra  of  a  Cd,  _^nxTe/CdTe  (x=0.21) 
strained-layer  single  quantum  wells  with  well  thicknesses 
of  Lw  ~  25,  50,  and  75  A  recorded  at  T =9.0  K.  The  PL 
spectrum  in  Fig.  10  was  measured  using  the  488-nm  line  of 
an  Ar+  ion  laser  as  the  excitation  source  when  the  detected 
wavelength  was  scanned.  The  circularly  polarized  excita¬ 
tion  spectra  in  Fig.  1 1  were  detected  at  the  longer  wave¬ 
length  side  of  the  peak  at  1.6  eV  of  the  PL  spectra  (Fig. 
10)  when  the  excitation  photon  (circularly  polarized 
<r+)  energy  from  the  cw  Ti-sapphire  laser  was  scanned. 
Curves  a  and  b  show  the  cr+  and  o~  components  of  the 
excitation  spectra,  respectively,  for  the  quantum  well  with 
Lw  =  75  A.  Curve  c  in  Fig.  1 1  is  calculated  by  Eq.  (5)  since 
the  signal  was  strong  enough  to  use  narrow  gates  (8=v/ 
5).  Detailed  discussion  of  the  PL,  PLE,  and  optica]  pump¬ 
ing  studies  of  Cd,  _iZnJtTe/CdTe  strained-layer  single 
quantum  wells  is  published  separately.7 
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The  electronic  states  in  a  series  of  molecular-beam-epitaxy-grown  CdTe/Cd0.wZna.2,Te  strained-layer 
'  single  quantum  wells  have  been  studied  both  theoretically  and  experimentally.  Calculations  based  on 
model  solid  theory  and  deformation-potential  theory  have  shown  that,  at  Zn  concentration  of  0.21  in  the 
barriers,  the  band  lineup  is  type  I  for  heavy  holes  and  type  II  for  light  holes.  The  confined  states  were 
calculated  using  the  three-band  envelope  function  model.  Photoluminescence  measurement  at  tempera¬ 
tures  between  9  and  100  K  indicates  that  the  broadening  mechanism  is  dominated  by  the  polar  optical 
phonon  (longitudinal  optical)  interactions.  The  circularly  polarized  excitation  spectra  have  shown 
predominant  heavy-hole  feature  confirming  that  tl  band  lineup  is  type  I  for  heavy  holes  and  type  II  for 
light  holes.  With  well  width  between  25  and  75  A,  our  study  has  shown  that  exciton  binding  energy  de¬ 
creases  with  decreasing  well  width. 


L  INTRODUCTION 

Recently,  considerable  interest  has  developed  in 
strained-layer  quantum  wells  and  superlattices  because  of 
their  potential  for  optoelectronics  device  applications. 
While  less  extensively  studied  than  III-V  compounds 
heterostructures,  II-VI  compounds  quantum  wells  and 
superlattices  are  attracting  more  and  more  attention  to¬ 
day.  They  offer  the  possibility  to  tune  the  energy  gap  to 
cover  the  spectrum  range  from  near  UV  to  the  far  in¬ 
frared.1-4  Besides,  the  strain  effect  splits  and  shifts  the 
energy  band  in  constituent  layers  and  hence  provides  a 


FIG.  1.  The  band  lineup  of  the  strained-layer 
CdTe/Cdo  7<jZn0  2,Te  quantum  well. 
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wealth  of  properties  with  which  to  design  devices  and 
study  material  physics.  In  the  case  of 
CdTe/Cdo .wZno  jiTe,  Harrison’s  linear  combination  of 
atomic  orbitais  (LCAO)  theory  indicates  that  it  has  a 
type-I  band  lineup  configuration.  When  the  strain  effect 
is  considered,  however,  it  may  end  up  to  be  simultaneous¬ 
ly  type  I  for  heavy  holes  and  type  II  for  light  holes.5 

This  paper  reports  on  the  photoluminescence  (PL), 
excitation  (PLE),  and  polarization  spectra  for 
CdTe/Cdo  79Zno  2|Te  strained-layer  single  quantum  wells 
with  the  band  lineup  configuration  of  Fig.  1.  Following 
this  introduction,  the  experiments  are  briefly  described  in 
Sec.  II.  Section  III  gives  the  calculation  of  confined  elec¬ 
tronic  states  in  the  strained-layer  single  quantum  wells 
(SL  SQW’s)  by  a  three-band  envelope  function  model  for 
strained-layer  structure  and  all  the  observation  of  spectra 
and  discussion  are  presented  in  Sec.  IV. 

H.  EXPERIMENTS 

Our  samples  were  grown  by  molecular-beam  epitaxy 
(MBE)  on  semi-insulating  <001 )  GaAs  substrates.  The 
structures  consist  of  a  4-^im-thick  Cd0  79Zn0  2i Te  buffer 
layer;  a  CdTe  single  quantum  well  with  the  nominal 
thickness  Lw  of  25,  50,  and  75  A;  and  a  0.5-/zm-thick 
Cdg  79Zn0  2,Te  cap  layer  (Fig.  1).  The  optical  spectra 
were  measured  in  the  backscattering  direction  at  ~26” 
from  the  normal  to  the  quantum-well  layer,  with  the  sam¬ 
ple  mounted  on  a  cold  finger  cooled  by  a  closed-cycle 
helium  cryostat  down  to  8.5  K.  The  PL  measurements 
were  performed  with  an  Ar+  laser  for  the  Cd0  79Zn0  21Te 
barriers  excitation  (488-nm  line)  or  by  a  tunable 
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FIG.  2.  Composition  dependence  of  Cd, -xZn,Te  alloy  band 
gap  at  6  K,  plotted  using  bowing  parameter  0.294. 


Thsapphire  laser  for  barrier  excitation  or  directly  the 
CdTe  quantum-well  excitation.  Luminescence  spectra 
were  recorded  with  a  Jarell-Ash  0.75-m  monochromator 
equipped  with  a  1180-grooves/mm  grating.  The  detec¬ 
tion  electronics  consists  of  thermoelectrically  cooled  pho¬ 
tomultiplier  with  a  fast  preamplifier  connected  to  a  dual¬ 
channel  gated  photon  counting  system  controlled  by  a 
computer,  which  also  controls  the  scanning  monochro¬ 
mator  Jarell-Ash  and  tunable  Thsapphire  laser.  The  cir¬ 
cular  polarization  of  luminescence  and  excitation  spectra 
were  measured  using  a  highly  sensitive  and  absolute  ap¬ 
paratus  involving  a  photoelastic  modulator  (PEM).  The 
system  will  be  described  in  detail  elsewhere.6 


a0  =  <  1  —  x)a(  AC)+xa(BC)  , 

A a0=a0(BC)—a0(  AC)  . 

a0's  are  lattice  constants  and  a,’s  are  the  intraband  (ab¬ 
solute)  hydrostatic  deformation  potentials  for  band  £,. 

For  variation  of  band  gap  £g(x)  with  composition  x  we 
have 


£g(x)=£c(x)-£„(x) 

=<1  -x)Eg(AC)+xEg(BC) 

+  3x(l— x)[oUC)— a(£C)]Aa0/a0  ,  (2) 

where  a  ( AC)  and  a  ( BC )  are  the  interband  (relative)  hy¬ 
drostatic  deformation  potentials:  a=ac—ae.  This  leads 
to  a  bowing  parameter 

b~l[a(AC)-a(BC)]£ia0/a0  .  (3) 


Here  aQ  chooses  the  midcomposition  lattice  constant. 
Note  that  in  (2)  and  (3)  we  ignored  the  strain  dependence 
of  spin-orbit  splitting  which  is  normally  insignificant  (see 
next  section).  For  the  ternary  alloy  CAq  21Te,  using 
the  parameters  given  in  Table  I,  Eq.  (3)  gives  b  =0. 378. 

According  to  Hill’s  theory,9  the  bowing  parameter  b 
can  be  calculated  from  the  equation 

V2 


b  = 


Ze 

8xreo 


1 


_1_ 

rB 


(r^+r^exp 


— sa- 


8 


III.  THEORY 


(4) 


A.  Interpolation  scheme  for  band  lineup 

In  ternary  alloys  A  ,_XBXC,  if  AC  and  BC  are  lattice 
matched,  linear  interpolation  is  a  valid  procedure  for 
variation  of  material  parameters  with  x.  But  in  ternary 
alloys  constituted  from  lattice  mismatched  binaries,  bew- 
ing  effects  should  be  taken  into  account.  One  material 
should  be  regarded  as  being  compressed  while  the  other 
being  expanded.  Therefore,  from  deformation  theory  we 
know  that  for  a  given  band  edge  £,  we  have7-8 

£.(x)  =  (  1  — x)£,(  AC)  +  xE,(BC) 

A  a0 

+  3x(l-x)[fl,(  AC)-a,(BC)) — -,  (1) 

ao 

where 


where  Ze  is  the  charge  of  the  substituting  ion,  a  is  the  lat¬ 
tice  constant  of  the  midcomposition  alloy,  s  is  the  screen¬ 
ing  factor  (0.25  A-1),  and  rA,  rB  are  the  Pauling’s  co¬ 
valent  radii  of  elements  A  and  B,  respectively. 

With  o(CdTe)  =  6.481  A  and  a(ZnTe)  =  6. 1037  A,4  the 
lattice  constant  of  the  midcomposition  alloy  is  calculated 
to  be  0=6.29235  A.  Using  =  1.405  A  and 
rZn  =  1.225  A,29  the  bowing  parameter  calculated  by  Eq. 
(4)  is  b  =0.294  (Fig.  2). 

The  experimental  dependence  of  the  energy  band  gap 
on  alloy  composition  (T  =  6  K)  given  in  Ref.  30  goes 

£g(Cd,_xZnxTe)=  1.606(1  — x)  +  2.394x 

—  0.26x(  1  — x)  eV  .  (5) 

Based  on  the  above  notation,  the  interpolation  scheme 


TABLE  I.  Parameter  values  used  for  calculations.  £t(Cd|_,Zn,Te)=  1.606  +  0. 494x -t-0.294x2 
(eV)  (calculated  using  Hill’s  theory),  o(Cd,-  ,ZnJITe)  =  6.481  — 0.3773x  (A),  from  Ref.  30  and  using 
m  *(Cd,  ,  Zn.Te)  =  0.099  +  0.01 7x,  y,(Cd,  1Zn,Te)  =  4. 1 1  — 0.04x,  y2(Cd, _,Zn,Te)=  1.08-0.3*, 
as  defined  in  Ref.  4. 


A 

<0 

a, 

b 

"o 

c„ 

C|  2 

(Ref  30) 

(Ref.  8) 

(Ref.  8) 

(Ref.  8! 

(Ref.  8) 

(Ref.  30) 

(Ref.  8) 

(Ref.  8) 

E  1 

*-1’  *V 

(eV) 

(eV) 

(eVi 

(eV) 

(eV) 

(A) 

C„  (Ref.  8) 

C12  (Ref  8) 

(eV) 

ZnTe 

2.364 

0  91 

079 

-5.83 

-1.26 

6.1037 

0.713 

0.407 

+  0.03 

CdTe 

1.60b 

093 

0.55 

-3.9b 

-1.10 

6.481 

0.562 

0.394 

+  0  21 

'Relative  to  the  top  of  GaAs  valence  band  given  by  HAO  theory  (Ref.  10). 
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should  be:  starting  from  Harrison  atomiclike  orbital 
(HAO)  theory’s  EUJAC)  and  £„av(£C),10  adding 
}AMC)  and  \MBC)  to  get  £„(  AC)  and  EV(BC),  and 
then  interpolating  £„(x)  by  Eq.  (1),  and  adding  the  exper¬ 
imental  Eg(x)  or  the  Eg(x)  given  by  the  bowing  equation 
onto  £„(x)  to  get  £c(x).  To  get  the  spin-orbit  splitting 
band  £so  <x),  the  linearly  interpolated  A(x)  is  subtracted 
from  £„(x).  This  gives  all  the  band  lineups  before  the 
strain  between  barrier  and  welt  layers  is  considered.  The 
valence-band  offset  used  above  is  in  good  agreement  with 
the  experimental  data  obtained  from  core-level  x-ray  pho¬ 
toemission  spectroscopy.10 

B.  Strain  effect  on  band  edges 

As  mentioned  before,  all  three  of  our  samples  have  a 
bufTer  layer  of  4  pm  which  is  much  thicker  than  its  criti¬ 
cal  thickness  so  that  the  strain  between  the  substrate  and 
the  buffer  is  relaxed,  a  cap  layer  of  0.5  /im  which  is, 
again,  much  thicker  than  the  quantum  wells  (25,  15,  and 
75  A).  Therefore,  it  is  quite  reasonable  to  assume  that 
only  quantum-well  layers  are  under  a  biaxial  compressive 
strain  which  is  equivalent  to  a  hydrostatic  strain  plus  a 
uniaxial  strain.  The  hydrostatic  component  of  the  strain 
shifts  the  energy  bands.  The  uniaxial  component  lifts  the 
J  =  -|  degeneracy  and  splits  it  into  a  pair  of  Kramers  dou¬ 
blets.11 

The  strain  Hamiltonian  Ht  for  a  p-like  multiplet  can  be 
expressed  asu 

£t=£in+£12)  . 

;n  =  -a,  Tr(e)-36,[<£;  - j£2)e„  +c.p.  ] 

-  ^3ds[(txLy  +tyLx  )cxy  +c.p.  ]  ,  (6) 

=  —  aAL-a )  Tr( e )  —  3b2[{Lxox  —  \L-a)txx  +  c.p.  ] 
-V3d2HLxoy+Lyoxkxy+c.p.]  , 


where  //j11  is  the  orbital-strain  Hamiltonian  and  is 
the  strain-dependent  spin-orbit  Hamiltonian,  e  and  ztJ 
denote  the  strain  tensor  and  its  components,  L  is  the  or¬ 
bital  angular  momentum  operator,  a  the  Pauli  matrix 
vector,  and  c.p.  stands  for  cyclic  permutation  with 
respect  to  the  indices  x,y,z.  a,  is  the  intraband  (absolute) 
hydrostatic  deformation  potential  while  6,  and  dx  are  or¬ 
bital  uniaxial  deformation  potentials  appropriate  to 
strains  of  tetragonal  and  rhombohedral  symmetries,  re¬ 
spectively.  a2,  b 2,  and  d2  are  additional  deformation  po¬ 
tentials  describing  the  effects  of  a  strain  on  the  spin-orbit 
interaction. 

The  effect  of  a  strain  on  the  conduction-band  edge  Ec 
is  to  produce  a  hydrostatic  shift  given  by 

/*'ef,=afTr(£)  ,  (7) 

where  ac  is  the  intraband  (absolute)  hydrostatic  deforma¬ 
tion  potential  of  the  conduction  band.  For  the  biaxial 
strain  in  (001 )  plane  we  have 


where  a ,  is  the  in-plane  lattice  constant,  a0  the  lattice 
constant  without  strain,  and  Cn,C12  are  elastic  con¬ 
stants. 

If  we  choose  \J,Mt )  basis  which  makes  spin-orbit  in¬ 
teraction  diagonal,  at  k  =0,  the  Kane  matrix  becomes13 


Ec0  +  bEH(. 
0 


£rjc+  -3  ~6£//.,  -bE, 


0 

0 


0 

0 


0 


0 


0 

0 


£|.af+y— 6£w,+6£,  '26  £,* 


V26  e: 


C  2A  VC- 

E,  „  +  -y  -&£«., 


(9) 


where 

6  E/i  t  =  at  ( 2  -  a  >r  , 

6  Efi  —  *  u  |  f  a  1 i(  2  a  )r  -  u  ( 2  —  A  k  , 

bE'lf  t  = » u  |  -  2a ,  It  2  -  a  )r  a ’( 2  ~  a )r  , 

6 £(  —  -  1 6,  4  26, )( 1  3  A.  Ir  =  6(1  +  A )r  , 

6 £’  -  -  '  6 1  -  6-i(l  +/.•!'  -  6 '( 1  +  A  )r  , 

A  3(  /<„  - 


r 


£  ,,  is  the  edge  of  the  conduction  band  without  strain  and 
E,  is  the  center  of  gravity  of  the  p-like  multiplet  given 
by  HAO  theory.1" 

Normally,  the  strain  dependence  of  spin-orbit  splitting 
is  not  very  significant  and  we  will  ignore  it  (a:^0, 
6,=0).  Therefore 


O'  = r<2,  , 

6 '  =  6 1  =6  . 

bE„  ~  6 Eji  =  u.  1 2  -  a  )r  , 
6 £  =  6 £'  =  6'  1  4  a  h  . 


(10) 
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The  values  for  a0  and  b  are  given  in  Table  I.  Under  these 
approximations,  the  strain  dependence  of  the  four 
energy-band  edges  of  interest  can  be  written 

EC=E,„+ j+Et+&EH'(  , 

Eht,=E,,MV  +  j-6EH'U-8Es  , 

•£ih  Ey  gv  j  &EHo  +  y8£, 

(11) 

+  {[A2+2A8E,+9<8£,)2],/2  , 

Es.o.  ~  Euty  —  8£Wu  +  j8£, 

-A[A2+2A8£J+9(8£,)2],/2  . 

C.  Three-band  envelope  function  modal 
for  strained-layer  quantum  wells 

To  calculate  the  confined  electronic  states  in  our  sam¬ 
ples,  the  strain  effects  discussed  above  need  to  be  incor¬ 
porated  into  the  three-band  envelope  function  model.14 
For  a  strained  system,  a  change  of  coordinates  for  Euler 
coordinates  rE  to  Lagrange  rL  is  necessary  to  preserve 
the  Hamiltonian  JTs  periodicity" 

rE=rL  +  u(rL),  (12) 

and  for  a  given  function  f  £(rE)~f L(rL )  we  have 

Pl/l^+^Pe/e  . 

(13) 

kL  =U+z)kE  . 


Following  the  procedure  in  normal  envelope  function  ap- 
proach,14, 15  we  end  up  with  an  equation 


cos  (k'wLJ+± 


1 


sin(fe^,Lu,)=0  . 


(14) 


This  has  the  exact  same  form  as  in  the  unstrained  case 
except  for  the  expression  of  tj  and  for  the  dispersion  rela¬ 
tion  in  each  material, 


( 1  +  bcw)k'wm'b(E) 
(\+btb)k'bm^E) 


(£-6£w.f  ])[(£  +  £,  +  6£H,,.-6£!) 

X(£  +  £g  +  A  +  6£a  —  2(8£s  >2| 

=  (\-4bt)P1h1k';1[E+Eg  +  4A  +  6£Wl  -5£,]  , 


(16) 


where 

8c  =  c. 


C 1 1  +  2  C .  i 


C„ 


c  , 


Ee,  A  are  energy-band  gap  and  s.o.  splitting  and 
m'hi  E),m'u\E)  are  given  by  the  relation 


£  = 


h2k ,2 

2m'(£)  ‘ 


(17) 


bEH  u,bEHc,bEs  are  defined  in  Eq.  (9).  The  results  here 
are  similar  to  that  of  Marzin  et  al.li  Due  to  the  strain 
effect,  our  samples  have  type-II  energy-gap  configuration 
as  shown  in  Fig.  1 . 

Our  calculated  optical  gap  (£,,—£,*)  and  experimen¬ 
tal  value  of  the  free-exciton  transition  1<-18  are  plotted  in 
Fig.  3  as  functions  of  well  width  Lw.  We  take  the 
difference  between  the  experimental  points  (A)  and  calcu¬ 
lated  optical  gaps  as  binding  energy  of  the  free  exciton 
(£4x  ).  As  predicted  by  theoretical  calculation,16- 19  Ebx 
first  increases  with  decreasing  Lw  because  of  the  tendency 
to  a  two-dimensional  behavior  and  then,  beyond  a  certain 
point,  Ebx  decreases  due  to  the  recovery  of  a  three- 
dimensional  behavior  through  the  increasing  penetration 
of  electron  and  hole  wave  function  in  the  barrier. 

The  three-dimensional  (3D)  exciton  binding  energies 
given  in  Ref.  20  are  12  meV  for  CdTe  and  13  meV  for 
ZnTe.  For  the  50-A-thick  CdTe  quantum-well  sample  in 
this  study,  the  value  of  Ebx  is  estimated  to  be  about  three 
times  larger  (36.5  meV)  jthan  that  of  the  3D  exciton.  In 
Ref.  5,  the  Ebx  for  a  50-A-thick  CdTe  quantum-well  sam¬ 
ple  is  also  reported  to  be  more  than  two  times  larger  (25 
meV)  than  that  of  a  3D  exciton.  The  main  difference  be¬ 
tween  the  samples  used  here  and  Ref.  5  is  in  the  barrier 
Zn  composition  x.  In  the  sample  used  here,  x  is  0.21, 
while  in  Ref.  5  x  is  only  0.08.  A  larger  x  value  can  have 
two  effects.  First,  a  larger  x  gives  a  wide  band  gap  in  the 
barrier  layer  and  hence  a  higher  potential  barrier. 
Second,  a  larger  x  value  in  the  barrier  layer  causes  larger 
strain  in  the  well,  pushing  both  electron  and  heavy-hole 
well  deeper.  Both  of  these  effects  enhance  the  two- 
dimensional  (2D)  feature  of  the  exciton  in  the  quantum 
well,  enlarging  the  binding  energy  Ebx-  The  accuracy  in 
the  deduced  Ebx  values  is  dependent  on  the  precision  of 
parameters  such  as  deformation  potentials  ac,  a „,  b,  and 
elastic  constants  Cn,  C12  used  for  calculations.  Different 
values  foi  these  parameters  have  been  used  in  litera- 


F1G.  3.  The  exciton  transition  energies  as  a  function  of  well 
width.  Solid  line  shows  the  calculated  transition  energies  ( lc- 
16!,  and  ▲  represents  the  experimental  points.  The  numbers 
beside  the  experimental  data  arc  the  difference  between  calcu¬ 
lated  and  measured  data  in  eV. 
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IV.  RESULTS  AND  DISCUSSION 

A.  PL  at  different  temperatures 

Figure  4  shows  the  PJ_  spectra  of  the  samples  with  well 
width  25,  50,  and  75  A,  respectively.  The  spectra  were 
measured  at  9.0  K  using  the  488-nm  Ar+  laser  line  for 
excitation.  The  excitation  is  indirect  (i.e.,  free  carriers 
are  photogenerated  in  Cd0  79Zno  2,Te  barriers  and  then 
transferred  into  CdTe  wells.)  As  expected,  the  sharp 
emission  peak  was  shifted  to  lower  energy  and  became 
narrower  as  the  well  width  was  increased.  The  main 
emission  peak  cannot  be  explained  by  a  single  line  (it  is 
neither  Gaussian  nor  Lorentzian).  Instead,  it  can  be 
fitted  gracefully  to  a  combination  of  two  lines,  one  Gauss¬ 
ian  and  one  Lorentzian  as  shown  in  Fig.  5.  The  least- 
square  fit  of  QW-3  (Lw  =  75  A)’s  spectrum  yields 
linewidths  of  7.553  and  6.65  meV,  centered  at  1 .6049  and 
1.6005  eV,  respectively.  In  the  bulk  binary  II-VI  semi¬ 
conductors,  the  linewidth  of  the  bound  exciton  is  smaller 
than  kT  since  the  bound  exciton  has  no  kinetic  energy. 
The  line  shape  is  described  by  the  Lorentzian  function.21 
The  free-exciton  line  shape  is  described  by  Gaussian 
function.  We  believe  that  the  Gaussian  peak  at  1.6049 
eV  originates  from  the  confined  n  =  1  electron  and  n  —  1 
heavy-hole  free-exciton  transition  and  the  Lorenzian 
peak  at  1 .6005  meV  is  attributed  to  the  excitons  bound  to 
neutral  acceptors  (£1A-A°)  in  a  CdTe  quantum  well. 
Similar  high-quality  fittings  were  also  obtained  for  the 
other  two  samples  except  that  sample  QW1  (Lw  =25  A) 
exhibits  one  additional  peak  on  the  lower  energy  side  of 
the  main  emission  peak,  which,  by  comparison  with  the 
PL  spectra  of  3-/xm  Cd0  7<>Zn0  2)Te  epitaxy  layer  grown  at 
the  same  conditions  as  the  buffer  layer  in  our  quantum- 
well  samples  [Fig  4(d)],  is  attributed  to  an  unidentified 
band  at  buffer-substrate  interface.  This  band  also  exists 
at  the  buffer-substrate  interface  of  the  other  two  samples 
(Lw  =  50  and  75  A).  However,  we  cannot  see  it  because 
the  optical  gap  (window)  in  these  two  samples  are  too 
narrow  to  let  this  emission  come  out.  Figure  6  shows  the 
PL  spectra  of  QW-1  (Lu.=  25  A)  at  different  excitation 
intensities.  The  spectrum  a  was  recorded  at  the  excita¬ 
tion  intensity  10  W/cm2  and  spectrum  b  at  40  W/cm2  us- 


Photon  Energy  (eV) 


FIG.  4  Photoluminescence  spectra  of  CdTe/Cd0  ,„Zn(,  ;,Te 
stngle-quantum-well  heterostructures  with  well  widths  (a)  25,  (b) 
50,  and  (ci  T5  A.  respective!),  under  indirect  excitation.  Spectra 
(dl  is  the  photoluminescence  of  Cd,,  7gZn„  ,,Te  epitaxial  layers 
shown  for  comparison. 


FIG.  5.  The  decomposition  of  the  PL  spectrum  for  (a)  QW-3 
single  quantum  well  to  two  components  spectra,  (b)  free-exciton 
Gaussian  shape,  and  (c)  bound  exciton  Lorentzian  shape. 


ing  an  Ar+  ion  laser  line  488  nm.  The  spectra  are  nor¬ 
malized  at  maximum.  The  additional  low-energy  band 
gets  saturated  at  a  higher  excitation  intensity.  Such  be¬ 
havior  suggests  that  the  additional  band  may  be  deter¬ 
mined  by  defects.22 

The  temperature  dependence  of  the  PL  spectra  of 
QW-3  (Lu.  =  75  A)  is  shown  in  Fig.  7.  The  photo¬ 
luminescence  is  quenched  mainly  due  to  the  thermally  ac¬ 
tivated  nonradiative  recombination  channel  in  the  well 
and  the  temperature  dependence  transfer  mechanism  of 
photogenerated  carriers  from  barrier  to  well.  The  inset 
shows  the  activation  energies  obtained  from  the  quench¬ 
ing  data  of  the  decomposed  free-exciton  emission  peak. 
Figure  8  shows  the  temperature  dependence  of  the 
linewidth  of  the  free-exciton  emission.  The  full  width  at 
half  maximum  (FWHM)  of  the  n  =  1  electron-heavy- 
hole  excitonic  emission  increases  slowly  in  the  tempera¬ 
ture  range  of  9-50  K,  but  increases  rapidly  with  temper¬ 
ature  above  50  K.  This  behavior  is  described  by  the  equa¬ 
tion23 


rlo.  fV  +  rph 


exp 


^LO 

kT 


+  r,mpexp 


—  (E) 
kT 


(18) 


where  the  first  term  Tq  is  the  linewidth  due  to  the  inho¬ 
mogeneous  fluctuation  of  the  well  thickness.  Second 
term  rph  describes  the  interactions  between  exciton  and 
longitude  optical  (LO)  phonon,  and  the  last  term  r,mp  is 
the  linewidth  due  to  the  ionized-impurity  scattering. 
(E  )  is  the  average  binding  energy  of  impurities.  Fitting 
the  above  equation  to  our  data  yields  the  solid  line  in  Fig. 


FIG.  6.  I’lioioluniinescencc  spectrum  for  single-quantum- 
well  QW-1  under  different  excitation  intensities 
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FIG.  7.  The  photoluminescence  spectra  of  a 
CdTe/Cdo.wZfio  jiTe  single-quantum-well  heterostructure  of 
well  width  Lm— 75  A,  measured  at  different  temperatures  be¬ 
tween  (a)  10  and  50  K  (step  10  K),  and  (b)  60  and  100  K,  under 
indirect  excitation  with  Ar+-ion  laser  line  488  nm.  The  inset 
shows  the  Arrhenius  plot  of  integrated  PL  intensity  versus  l/T 
which  reviles  two  activation  energies  25  and  5  meV. 


8  with  the  parameters  To  =75  meV,  rph=400  meV, 
£lo = 22  meV,  rjmp = 500  meV,  and  <  E  >  =  92  meV. 


B.  Circularly  polarized  excitation  spectra 

Since  the  calculation  has  shown  that  the  band  lineup  is 
type  I  for  heavy  holes  and  type  II  for  light  holes,  the 
main  emission  peak  is  expected  to  be  characterized 
predominantly  by  the  heavy-hole  feature,  which  is  clearly 
shown  in  our  circular  polarization  optical  pumping  spec¬ 
tra  (Fig.  9).  The  degree  of  circular  polarization  p  is 
defined  as  24 


P  = 


/++/-  ’ 


(19) 


FIG.  9.  (A)  The  polarized  excitation  spectra  excited  by  circu¬ 
larly  polarized  o~.  The  polarized  spectra  (a)  a*  and  (b)  a  are 
obtained  when  detecting  at  the  longer  wavelength  side  of  free- 
exciton  transition  le-1  h.  Curve  c  is  the  degree  of  circular  polar¬ 
ization.  (B)  Phonon-assisted  cooling  processors  for  hot  elec¬ 
trons. 


where  /+  and  /_  are  the  intensities  of  right  and  left  cir¬ 
cularly  polarized  (<r+  and  a _ )  radiation. 

In  bulk  material,  according  to  the  selection  rule  shown 
in  Fig.  10(a),  the  maximum  degree  of  circular  polariza¬ 
tion  that  can  be  expected  is  only  0.25  because  of  the  de¬ 
generacy  at  the  top  of  the  valence  band.23 

In  the  case  of  SL  SQW’s  or  SL  SL’s,  on  one  hand  the 
removal  of  the  degeneracy  by  strain  and  confinement  sug¬ 
gests  that  degree  of  circular  polarization  as  high  as  100% 
could  be  observed.  On  the  other  hand,  however,  the  mix¬ 
ing  of  the  states  due  to  the  in-plane  dispersion  lowers  the 
degree  of  polarization.  Furthermore,  the  spin  relaxation 


/=>/2 


/=3/2 


Temperature  (K) 

FIG.  8.  Temperature  dependence  of  linewidth  (FWHM)  of 
free-excitonic  emission.  The  solid  curve  is  the  best  fit  to  Eq. 
(18).  Experimental  data  are  indicated  by  ▼. 


FIG.  10.  (a)  Selection  rules  for  optical  transitions  and  rela¬ 
tive  intensities,  (b)  The  excitation  and  observation  geometry  in 
relation  to  sample  orientation. 


46 


OPTICAL  PROPERTIES  OF  CdTe/Cd,  ^,Zn,Te  STRAINED- . . . 


6967 


of  the  photogenerated  carrier  also  reduces  the  degree  of 
polarization.  And  finally,  the  orientation  of  the  sample 
relative  to  the  direction  of  the  excitation  beam  and  obser¬ 
vation  direction  also  affects  the  observed  degree  of  polar¬ 
ization.1425 

2(  2  coso  cosS 

- ! - —  =  -  zcosac°sP  ,  (20) 

l  +  (v-n)2  l+cos2a 

where  n,  n , ,  and  v  are  unit  vectors  in  the  direction  of  ex¬ 
citation  beam,  observation,  and  sample  growth  respec¬ 
tively  [see  Fig.  10(b)]. 

The  circularly  polarized  excitation  spectrum  of  QW-3 
(Lu=75  A)  is  shown  in  Fig.  9(a).  The  transition  moni¬ 
tored  is  the  long-wavelength  side  of  the  free  exciton  ( le- 
\h )  at  1.605  eV.  The  excitation  source  was  a  continuous 
wave  (cw)Ti:sapphire  tunable  laser.  The  light  from  the 
laser  was  circularly  polarized  (o~)  by  an  achromatics 
circular  polarizer.  The  detecting  system  record  simul¬ 
taneously  in  two  gated  photon  counters  the  a  4  and  o  ~ 
components  of  the  luminescence  light  which  are  curves  a 
and  b,  respectively.  Curve  c  is  the  circularly  polarized 
excitation  spectrum  calculated  using  Eq.  (19).  Similar 
spectra  were  obtained  for  all  the  SQW  samples  studied. 
The  forbidden  transition  (le-2/i,  \e-3h,  etc.)  energies,  ac¬ 
cording  to  our  calculations,  happen  to  be  in  the  valleys  of 
the  excitation  spectra  shown  in  Fig.  9(a).  Therefore,  the 
peaks  in  Fig.  9(a)  at  1.683,  1.663,  and  1.639  eV  are  related 
to  the  LO  phonon  replicas  (bulk  CdTe’s  LO-phonon  ener¬ 
gy  is  21  meV)  as  indicated  in  Fig.  9(b). 

As  one  can  see  in  Fig.  9(a),  when  the  excitation  energy 
hv  is  above  the  barrier  (but  smaller  than  E*  +A  of  the 
barrier),  most  of  the  spin  orientation  of  photogenerated 
carriers  created  by  a  +  excitation  are  relaxed  in  the  pro¬ 
cess  of  transferring  to  the  quantum  well  and  subsequent 
thermalization  (or  cooling)  in  the  quantum  well  before 
recombination,  and  hence  a  very  low  degree  of  polariza¬ 
tion  is  observed.  When  the  excitation  energy  hv  is  below 
the  barrier  such  that  [h  v— (£u +£g +£lhh  )]mhh/ 
(mhh  +m, )  is  a  multiple  of  the  LO  phonon  energy  h<oLO, 
the  excited  hot  electrons  emit  LO  phonons  and  fall 
directly  to  the  bottom  of  the  subband  retaining  most  of 
their  spin  orientation.26  The  gradual  decrease  of  the  de¬ 
gree  of  luminescence  polarization  with  increasing  number 
of  emitted  LO  phonons  indicates  noticeable  spin  relaxa¬ 
tion  at  each  stage  of  the  cascade.  When  the  condition 
(liv-(£|,  +  £|  +  £lhh)]mhh/(mhh  +  mf  )  =  nhioLO  is  not 
fulfilled,  the  cooling  process  always  involves  emission  of  a 
great  number  of  acoustic  phonons,  and  spin  orientation 
of  the  hot  electrons  is  significantly  relaxed.  When  the  ex¬ 
citation  photon  (o+)  energy  hv  is  resonant  with  the 
emission  peak,  the  degree  of  circular  polarization  as  high 
as  73%  is  observed. 


C.  Collection  of  photogenerated  carriers 
into  the  wells 

Assuming  a  uniform  absorption  coefficient  of 
5X10~Vcm  for  Cd0  -9Zn0  ;iTe  at  10  K  at  488  nm,  99% 
of  the  incident  photons  during  a  PL  experiment  are  ab¬ 
sorbed  within  the  0.5-/im-thick  cap  (barrier)  layer. 


though  thin  quantum  wells  give  very  intense  PL  Uw  is  50 
times  stronger  than  Ib).  Evidently,  most  of  the  carriers 
generated  in  the  barrier  are  transferred  into  the  well  be¬ 
fore  recombination  occurs.  The  carrier  collection  dy¬ 
namics  of  SQW’s  have  been  studied  by  many  groups.27,28 
The  simple  lumped  continuity  equations  are  written 


+  f 
Tnr  b  Tlr 


g  d  +  — = 

Sw  to  _ 


^nu  I'nrw 


(21) 


where  the  subscript  r,  nr,  and  tr  denote  radiative,  nonra- 
diative,  and  transfer;  b  and  w  denote  the  barrier  and  well, 
respectively;  nb  and  nw  are  the  two-dimensional 
minority-carrier  concentrations;  gbdb  and  gadm  are  the 
two-dimensional  generation  rate  in  the  barrier  and  the 
well;  and  r’s  are  the  associated  minority-carrier  lifetimes. 
The  tu  is  the  time  needed  to  transfer  photogenerated  car¬ 
riers  from  a  5000-A  barrier  into  the  well.  The  internal 
radiative  recombination  rate  per  unit  area  Ib  and  Iw  are 
given 


(22) 


From  Eqs.  (21)  and  (22)  one  can  obtain 
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(23) 


The  ratio  Iw  /Ib  can  be  approximated  by  the  ratio  of  in¬ 
tegrated  PL  intensities  if  we  assume  identical  conversion 
factors  linking  the  external  and  internal  radiative  recom¬ 
bination  rates  for  the  well  and  the  barrier.  It  is  quite 
reasonable  to  assume  t}w  =  l.28  Neglecting  of  carrier  gen¬ 
eration  in  the  well  simplifies  (23)  to 


(24) 


Now,  the  time  needed  to  transfer  the  photogenerated  car¬ 
riers  from  the  5000-A-thick  Cd0  79Zno  21Te  barrier  to  the 
CdTe  well  can  be  estimated  if  the  radiative  minority- 
carrier  lifetime  in  the  barrier  is  known. 


J-VwTu~ 


(25) 


Assuming  that  the  radiative  lifetime  in  the  barrier  is 
the  same  as  in  bulk  ZnTe,  about  ~2  ns,  and  the  quantum 
efficiency  in  well  is  ijw—  1,  the  estimated  transfer  time  rlr 
is  40  ps.  If  quantum  efficiency  t]u,  is  not  equal  to  unity, 
the  transfer  time  should  be  shorter.  From  a  Streack 
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Camera  measurement  under  direct  excitation,  the  emis¬ 
sion  from  the  well  reached  its  peak  about  3S0  ps  after  ex¬ 
citation  pulse.  That  time  is  the  sum  of  the  transfer  time, 
the  exciton  creation  time,  and  radiative  lifetime  of  the  ex- 
citon  in  the  well  (rlf  +  rcl+rTO ).  The  radiative  lifetime  of 
exciton  in  the  well  is  about  310  ps  if  we  ignore  the  exciton 
creation  time  rcx. 

V.  CONCLUSION 

In  conclusion,  we  have  reported  theoretical  calcula¬ 
tions  and  experimental  results  on  CdTe/Cdo ^^no  2|Te 
strained-layer  single  quantum  wells.  We  have  presented 
theoretical  calculations  for  the  confined  carrier  energies 
using  the  three-band  envelope  function  model  including 
the  strain  effect.  Based  on  Hill’s  and  model  solid 
theories,  we  have  computed  the  band  gap  of  the 
Cdo  79Zn0  2|Te  alloy  as  a  function  of  composition  x.  The 
calculated  bowing  parameters  are  0.294  and  0.303  from 
Hill’s  theory9  and  model  solid  theory,8  respectively. 

The  exciton  binding  energies  were  estimated  from  the 
difference  between  the  calculated  \e-\h  transitions  in  well 
and  experimental  free-exciton  transitions  (Fig.  3).  The 
binding  energy  is  reduced  from  36.5  to  12.9  meV  when 
well  width  is  reduced  from  75  to  25  A.  This  behavior  is 
predicted  by  theory.  In  the  range  of  investigated  quan¬ 
tum  well  widths,  the  binding  energy  Ebx  decreases  with 


reduction  of  the  well  width.  This  is  due  to  the  recovery 
of  a  three-dimensional  behavior  through  the  increasing 
penetration  of  the  electron  and  hole  wave  functions  in  the 
barrier. 

The  analysis  of  the  free-exciton  linewidth  as  a  function 
of  temperature  indicated  that  the  broadening  mechanism 
is  dominated  by  the  polar  optical  phonon  (LO)  interac¬ 
tions.  The  temperature  dependence  of  the  integrated  in¬ 
tensities  of  free  excitons  yield  activation  energies  of  25 
meV  in  high-temperature  regions  and  5  meV  in  low- 
temperature  regions.  The  activation  energy  25  meV  cor¬ 
responds  to  the  dissociation  of  the  free  excitons. 

The  polarized  excitation  spectra  have  shown  a  predom¬ 
inant  heavy-hole  feature.  No  light-hole  transitions  were 
observed.  This  confirms  that  the  lineup  configuration  is 
type  I  for  heavy  holes  and  type  II  for  light  holes. 
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